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Preface 


Impedance spectroscopy is a measurement method applied in many fields of science 
and technology, such as electrochemistry, material science, biology and medicine. The 
measurement method is based on a system identification approach and uses excita- 
tion signals of different frequencies to get more information about the system under 
test. It provides interesting possibilities to use the complex impedance of a material 
or a system over a wide frequency range. Thereby, effects dominating at different fre- 
quency ranges can be separated, non-accessible quantities become measurable, accu- 
rate measurements can be realized and multifunctional measurements can be carried 
out at low costs. 

For Impedance Spectroscopy (IS), competency in several fields of science and 
technology is indispensable. An understanding of electrochemical and physical phe- 
nomena is necessary for modeling and understanding such phenomena. Measure- 
ment procedures should be developed taking specific requirements of the considered 
application into account. Signal processing methods are very important for extract- 
ing target information by suitable mathematical methods and algorithms. New trends 
are emerging rapidly, involving special techniques for realizing fully automatic em- 
bedded solutions at low costs and requiring a deep overview of modern information 
technology. 

This book addresses impedance spectroscopy and its advanced applications in 
the fields of battery research, bioimpedance, measurement systems and sensors. It 
begins with a comprehensive review on impedance spectroscopy advances and future 
trends. Then, in three parts, it reports about new advances and different approaches 
in dealing with impedance spectroscopy, including theory, methodology and applica- 
tion aspects. 

In the first part, the focus is on electrochemical and physical phenomena in the 
fields of nanocomposites, fuel, battery cells, nanoparticles and alloys. In the second 
part, measurement systems for impedance spectroscopy are discussed, including un- 
certainty, measurement duration, simulation, system implementation and interface 
circuits. In the third part, different sensors are presented, including capacitive sensors 
for detecting aqueous solutions, eddy current sensors for corrosion measurement of 
steel, humidity sensing and nitride sensing. In the fourth part, about bioimpedance 
spectroscopy, the topics are manifold and cover deep brain stimulation, choice of the 
number frequencies for measurement, intraventricular bioimpedance spectroscopy 
and meat quality monitoring by classification of impedance spectra. 

The book scope is representative for the versatile measurement method with its 
highly interesting facets. In addition to recent advances concerning the main appli- 
cation fields of impedance spectroscopy, the book includes many methods and reflec- 
tions on theory, modeling and implementation aspects. This book is therefore very 
interesting for researchers and developers in the field of impedance spectroscopy. 


https: //doi.org/10.1515/9783110558920-201 
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I thank all contributors for the interesting contributions and the reviewer who 
supported the development of this book. We hope together to contribute actively to 
further developments in this field by sharing our scientific and practical knowledge. 


Prof. Dr.-Ing. Olfa Kanoun 
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Olfa Kanoun 
Impedance spectroscopy advances and future 
trends: A comprehensive review 


Abstract: Impedance spectroscopy is a highly performant measurement method, 
which is often used in many fields of science and technology. It refers to the measure- 
ment ofthe complex impedance over a range of frequencies and provides possibilities 
for separating relaxation effects, accurate measurements and measurements of non- 
accessible quantities. In material science, impedance spectroscopy allows the charac- 
terization ofthe complex conductivity. Impedance spectra of battery cells allow to get 
information in a non-destructive way and to separate different phenomena concern- 
ing electrolyte and electrodes.In the biomedical field, it has the decisive advantage by 
providing a lot of information non-invasively and for long time intervals. For sensors, 
decisive improvements of measurement accuracy can be reached and multifunctional 
sensors can be realized, which allow measurements of more than one measurement 
quantity at the same time and same state of a material or a device. 

In this contribution, after an introduction about the fundamentals of the measure- 
ment method, the links and differences between a wide range of application areas of 
impedance spectroscopy are shown for the first time. In each of these fields, we high- 
light the specific strengths of impedance spectroscopy and report on recent advances 
and future trends. 


Keywords: Impedance spectroscopy, electrochemical impedance spectroscopy, relax- 
ations, battery diagnosis, bioimpedance, sensors, biosensors, inductive sensors, ca- 
pacitive sensors, nanocomposite, corrosion, signal processing 


1 Introduction 


The complex impedance of a material or a system provides information about its com- 
position, interfaces and geometry and helps to follow changes of its state. Impedance 
measurement has therefore decisive benefits compared to other characterization 
methods, as it is experimentally efficient, because of the relative reduced time, hard- 
ware and software needed. It delivers more information than only resistive, capacitive 
or inductive measurements, because it considers at the same time both the real and 
the imaginary part of the complex impedance. It provides the possibility to separate 
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effects having linearlyindependent dispersions and dominating in different frequency 
ranges. 

Impedance spectroscopy has shown a high potential for many years and has so far 
demonstrated its potential in a big number of applications, such as characterization 
of energy storage units [1], biomedical and biological applications [2] , [3] , [4], mate- 
rial testing [5-7], corrosion and coatings [8] and inductive [9, 10], capacitive [11] and 
conductive sensors [12]. The range of applications is very broad and the advantages it 
provides are significant under the condition that it is applied in the right manner. In 
every application, field-specific strengths of the method are in the focus of interest, 
even if the method can provide more advantages. 

Impedance spectroscopy is non-invasive, as it uses small amplitudes of the exci- 
tation signal. It is capable to achieve real-time measurement by using wide frequency 
band signals. Even if it is generally known as a method for materials and systems 
characterizations in laboratory environments, impedance spectroscopy is increas- 
ingly used for real-time measurements in portable and embedded solutions due to the 
favorable developments in the field of microelectronics. 

In this contribution, we give an overview about different aspects of the method in 
relation to a wide range of application fields. The aim is to demonstrate the potential of 
impedance spectroscopy and to show the link between the different application areas. 
In each context, the reasons for the particular strengths of the method are discussed 
and special potentials are shown. Then, recent developments are explained by means 
of selected examples. 


2 Theoretical background 


Impedance spectroscopy belongs to the so-called system identification techniques, 
where the frequency-dependent transfer function of a system is determined. In 
impedance spectroscopy, the complex impedance of a material or a system is mea- 
sured at different frequencies or in a certain frequency range. Thereby both the real 
and the imaginary part of the complex impedance can be considered as useful sources 
of information, where certain phenomena are dominating in certain frequency ranges 
and can therefore be characterized. In general, the real and the imaginary part of the 
impedance show a different frequency dependency and therefore provide a compre- 
hensive source of information. This is why manifold information regarding electrical 
and electrochemical properties of the investigated system or material can be extracted. 

A wide frequency band helps to gain more information at a given working point 
of the i-u characteristic (see Fig. 1 (a) and (b)). If the measurements carried out corre- 
spond to the same phenomena (see Fig. 1 (c)), impedance spectroscopy helps to real- 
ize a better characterization of the phenomena. The measurements are from a phys- 
ical point of view redundant and provide possibilities to have a better measurement 
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Fig. 1: Principle of impedance spectroscopy. (a) Excitation with a single frequency at a given working 
point. (b) Example for excitation with a varied frequency with a sweep. (c) Improvement of measure- 
ment information by redundancy. (d) Improvement of measurement information by diversity trough a 
wide band excitation signal. 


accuracy, as it becomes possible to reduce the influence of random deviation, such 
as noise, and systematic deviations, e. g., temperature, and to benefit from a detailed 
modeling of the characteristic. 

If the impedance measurement is carried out in a wider frequency range (see 
Fig. 1 (d)), the information gained thereby becomes more probable to have diversity, 
as different physical/chemical phenomena can be dominating in different frequency 
ranges. If different spectral behaviors are included in the impedance spectrum refer- 
ring to different phenomena, impedance spectroscopy provides interesting opportu- 
nities to measure more than one quantity at the same time. The conditions for this 
are (i) to have sufficient measurement data in the spectral regions where the indi- 
vidual effects are dominating and (ii) for the time constants corresponding to these 
phenomena to be sufficiently distant. 


3 Extraction of information from impedance spectra 


Several aspects should be specifically addressed on the design process for an imped- 
ance measurement system, such as measurement procedures, investigations of physi- 
cal and chemical phenomena taking place, development of suitable impedance mod- 
els and extraction of target information by optimization techniques. 

In general, there are different approaches, which can be considered for informa- 
tion extraction from impedance spectra. Physical methods are interesting and allow 
profound insights in a material or a system (see Fig. 2). Mathematical methods have 
also been proposed, such as data mining, principal component analysis [13], neural 
networks and fuzzy logic [14]. 

For physical modeling, it is necessary to understand the physical and chemical 
phenomena taking place within the system and to formulate the necessary model 
equations by combining the descriptions of partial phenomena, e. g., in series, par- 
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Fig. 2: General approach of impedance spectroscopy and for information extraction from impedance 
spectra. 
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allel or cascade. Thereby, different formulations can be mathematically equivalent, 
but only some of them provide a meaningful model and are therefore physically rele- 
vant. The extraction of information is carried out, in this case, by determination of the 
model parameters, which are directly related to the corresponding physical phenom- 
ena. If the formulation of the model is suitable and the parameter extraction is suc- 
cessful [15], the model parameter provides information about the measured quantity 
and shows a physically expected behavior. In this case, extracted parameters inferring 
the measurement quantity are directly used for the calculation of the measurement 
quantity. 

Mathematical methods are also possible, but need a big data base with relevant 
experiments in order to build the necessary basis for signal processing and extraction 
of information. Suitable feature extraction and classification methods need to be de- 
veloped in this case to discriminate between the different classes corresponding to the 
state of the system [14]. For this purpose, the training data and the experimental data 
need to be representative for the real scenario and include all relevant cases that could 
influence the results, such as environmental parameters, aging and contamination ef- 
fects. 


4 Impedance spectroscopy in measurements and 
sensor solutions 
Impedance spectroscopy provides the possibility to follow the state of a system and 


to track its changes, e.g., due to aging or chemical reactions taking place. This is 
the case, e. g., for corrosion, battery diagnosis and bioimpedance. Especially for bat- 
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tery diagnosis and bioimpedance analysis, impedance spectroscopy allows to non- 
invasively get information about the state of a cell or a battery and to determine spe- 
cific information concerning parts of it or phenomena taking place within it. In the 
case of conductive, capacitive and inductive sensors, impedance spectroscopy allows 
to measure quantities in spite of experimental limitations, to improve the accuracy or 
even to measure more than one quantity at the same time. 

The applications for impedance spectroscopy are therefore manifold. In the next 
sections, we present a selection of the most important applications for impedance 
spectroscopy (see Fig. 3) and show the potential of this measurement method in the 
corresponding fields in comparison with other classical methods. 
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Fig. 3: Selection of impedance spectroscopy applications in measurement and sensor solutions. 


5 Material characterization 


For materials, impedance spectroscopy provides interesting possibilities to character- 
ize the complex electrical behavior including resistive, capacitive and inductive be- 
havior [5]. Thereby, it is sensitive to both surface phenomena and changes of bulk 
properties including conduction mechanisms. 

For example, impedance spectroscopy is a well-established method, which allows 
quantitative characterization of corrosion processes as well as the evaluation of the 
anti-corrosion performance of protective coatings on corrosive metals [16]. It can be 
used for several purposes, among which the detection of surface and material changes 
due to exposure, prediction of the lifetime of corrosion protection, evaluation of coat- 
ing systems, identification of the corrosion processes that lead to failures, measure- 
ment of water uptake by coatings and the development of models for coating/metal 
system performance. 

One of the important advantages of this method is that it is a non-destructive tech- 
nique. It can follow the evolution of coated metals as they undergo either natural or 
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artificial exposure to conditions that cause corrosive failure to provide information 
about the corrosion kinetics. Impedance measurements allow the prediction of the 
long-term performance of the coatings and track very effectively the effects of tem- 
perature, water ingress, coating plasticization and the glass transition temperature in 
coatings [17]. 

Impedance spectroscopy is interesting for electrical characterization of almost ev- 
ery material. It has played a big role for many years in the field of nanocomposites, 
which are becoming increasingly important in the field of lightweight structures and 
sensors. Flexible films of carbon nanotubes (CNTs) can be prepared for several sen- 
sor applications, such as strain gauges or pressure sensors showing a high sensitivity 
and performing under low costs. With less than 1% filler material, input conductive 
and percolation-based CNT polymer composite sensors can be fabricated. Thereby, a 
suitable composition and homogeneity are highly demanded and should be fulfilled 
by the fabrication process including choice of materials, preparation of dispersion, 
substrate treatment and deposition techniques. 

Impedance spectroscopy enables to deeply investigate the complex conductiv- 
ity of the nanocomposite materials, so that they can be better optimized for certain 
required properties [5, 18]. For example, for highly sensitive pressure sensors based 
on CNT polymer nanocomposites, it is very important to realize nanocomposites at 
the percolation ratio [19]. The percolation ratio is the nanofiller concentration at 
which rarely conductive paths can be built through the nanofillers without tunneling 
through the polymer. At this concentration, the nanocomposite becomes more sensi- 
tive to pressure, because an increase of pressure leads to a reduction of the distance 
between the CNTs and therefore increasingly leads to percolation. These concen- 
trations cannot be determined easily by experimental methods due to fluctuating 
processing parameters during the fabrication process. 

Impedance spectroscopy is suitable for characterization of the percolation behav- 
ior, as it includes the information about both resistive and capacitive behaviors of the 
nanocomposite [5]. Especially the imaginary part of the impedance gives information 
about the average gap between nanofillers and is therefore very important for investi- 
gating percolation (see Fig. 4). 

In [5] it has been shown that impedance spectroscopy measurements allow the 
experimental determination of the percolation threshold of the nanocomposite. At 
concentrations lower than the percolation ratio, the nanocomposite shows an almost 
purely capacitive behavior. Due to inter-tube distances, the share of capacitive behav- 
ior dominates the conduction behavior (see Fig. 4 (a)). Beyond the percolation ratio, 
the conductive paths dominate (see Fig. 4 (b)) and show a frequency independent be- 
havior. The percolation region is exactly the region where the transition between these 
two extreme cases takes place abruptly. 

This example shows that, by means of impedance spectroscopy, the percolation 
ratio can be determined experimentally even after film deposition. This allows tailor- 
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Fig. 4: CNT polymer composite. (a) Re- 
laxed: tunneling effect, the capacitive 
behavior is dominating. (b) Under pres- 
sure: resistive behavior is dominating 
due to increasing contacts between 
CNTs. 


ing of the piezoresistivity of strain-sensitive films based on nanocomposites for differ- 
ent nanofillers and even by fluctuating processing parameters. 


6 Battery characterization and diagnostics 


Modern battery management systems need accurate information concerning the state 
of charge (SoC) [20], the state of health (SoH) [21], the battery lifetime, the internal 
temperature and fault diagnosis within cell packs. Other methods are, e.g., voltage 
measurement and coulomb counting, which bring some challenges in the case of Li- 
ion batteries. In particular, the SoC-vs-voltage curve is very flat, so that the sensitivity 
of the voltage to the SoC is very low. The direct voltage measurement is not sufficient 
for a comprehensive diagnosis as it is dependent on SoC, temperature and current 
level at the same time. In general, the open-circuit voltage (OCV) provides an alterna- 
tive, which can be determined at offline state, where the battery is isolated from the 
external circuit for a long time. In online measurement, it can be estimated from the 
voltage measurement by considering the internal resistance multiplied by the current 
[22] Coulomb counting integrates the current over time to estimate the remaining bat- 
tery charge. But, due to the limited sampling rate, several current impulses may not 
be tracked with a high accuracy. Therefore, the coulomb counting drifts over time due 
to the cumulative deviation. 

In contrast, impedance spectroscopy is a non-destructive method, which can es- 
timate the remaining charge of the battery, also known as SoC and the state of health 
(SoH) of the battery. This is due to the fact, that battery impedance increases with ag- 
ing and changes due to temperature and also different phenomena taking place at low, 
middle and high frequencies. Therefore, impedance spectroscopy can be used as a di- 
agnostic tool and accurate analytical technique for measuring different critical battery 
parameters. It provides the possibility to get information about the inner state of the 
battery, including both the electrolyte and the electrodes, non-invasively and in a rel- 
atively short measurement time, without needing to carry out any cycling processes. 
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Impedance spectroscopy offers a solution that is independent on the leakage currents 
and the fluctuating voltage. This method is experimentally efficient and provides the 
possibility to separate effects having linear independent dispersions and dominating 
in different frequency ranges. It is capable to achieve real-time measurement in a rel- 
atively short measurement time as it can use wide frequency band excitation signals. 
As a battery ages, several phenomena occur, accompanied by a change of the over- 
all impedance of the battery [22-24]. The SoC of the battery can be also estimated by 
means of impedance spectroscopy [25-28]. 

The battery impedance spectrum includes information about the electrolytes, dif- 
fusion layer and electrode reactions, as explained in Fig. 5. At low frequencies, in the 
range of mHz, processes with a high time constant are stimulated, such as migration 
and diffusion. Diffusion effects are dominating in this frequency range. At middle fre- 
quencies, electrode reactions dominate and may show a constant phase behavior due 
to electrode porosity. At high frequencies, inductive effects in the battery wiring and 
porous structure are dominant. 
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Fig. 5: Typical impedance spectrum of a battery in an exaggerated Nyquist plot. Here, R, is the se- 
ries resistance, R, and R, are the charge transfer resistances of the electrodes, Can and Can are the 
double layer capacitances of the electrodes, ø is the Warburg coefficient, Rz is the intersection of the 
diffusion characteristic with the real impedance axis and L represents inductive effects. 


The individual phenomena in the impedance spectrum of a cell have their own mod- 
els, which can be expressed analytically or in form of equivalent circuits. Models based 
upon concentrated circuit elements such as capacitors, resistors and inductors can 
simulate some internal electrical properties of the cell. Nevertheless, especially for 
diffusion processes and due to electrode porosity, elements with distributed parame- 
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ters, such as constant phase elements and Warburg impedance, are also necessary to 
model the impedance behavior accurately. 

The impedance spectrum of a battery cell includes manifold information and the 
method is therefore very interesting for diagnosis. Nevertheless applying it on a bat- 
tery cell involves many challenges. The low-frequency measurements at mHz frequen- 
cies need a long measurement time and involve therefore a risk of charging or dis- 
charging the battery cell and changing its state thereby. The battery impedance is 
low and needs to be measured with a high resolution. A big challenge in modeling 
impedance spectra is that the individual phenomena may overlap in the impedance 
spectrum. This overlapping is strongly dependent on the cell technology and also on 
the cell state. In the case of overlap, the separation of effects and the parameter calcu- 
lation become more difficult, especially if the measured frequency is limited. The AC 
current distribution within cylindrical Li-ion cells depends strongly on the frequency 
and the temperature [29]. This may lead to limitations for deep interpretation. The cell 
design may overshadow some active material impedance by local inhomogeneity due 
to lithium plating or advanced aging. The spectrum may change its structure due to 
aging and models need to be adaptively adjusted to consider this effect. 

Yet, the latest developments allow to realize embedded solutions for impedance 
spectroscopy, which realize a comparable measurement quality like sophisticated lab- 
oratory devices [21, 30] having low cost and a significantly short measurement time. 
This development allows to use impedance for real-time measurements on site and 
to determine useful information for intelligent battery management systems, such as 
SoH [22], SoC [20], the internal cell temperature [31] and the remaining useful battery 
life [32]. 


7 Biomedical applications and food characterization 


For non-invasive tissue characterization, pathological health status investigation and 
diseases diagnosis, bioimpedance spectroscopy is one of the most important methods. 
In comparison with classical non-invasive techniques used in clinical research, such 
as optical spectroscopy, ultrasound and X-ray, bioimpedance spectroscopy is easy to 
use and safe. It can give information on the deepness of a material and helps to avoid 
interferences with substances reacting in the same wavelength. Therefore, impedance 
spectroscopy is gaining importance in the biological and food sectors and will play a 
more important role for the food industry in the future. 

Impedance spectroscopy can be used for characterization of different biological 
tissues and for medical [33-37] and food applications [11, 14, 38-41], because the com- 
plex impedance gives information about the composition and the state of a biological 
tissue or organ. Although there are many different methods, impedance spectroscopy 
is often preferred. 
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The impedance of a muscle includes three areas of dispersion: the a dispersion at 
low frequencies between some Hz and some kHz, ß dispersion at middle frequencies in 
the range from some kHz to MHz and y dispersion at higher frequencies [14, 33, 38, 39]. 
The £ dispersion is very informative about the state of the muscle as it provides in- 
formation about the cell membranes. Therefore, impedance spectroscopy is a very 
promising method for assessment of meat freshness and meet characterization in gen- 
eral. 

A lot of challenges should be overcome in the measurement of muscle impedance 
in order to guarantee reproducible measurements. For this, an electrode configuration 
has been developed to avoid influence of anisotropy and in order to realize a good 
contact to the biological material in spite of its fibrous structure [4, 42]. The devel- 
oped probe consists of one central cylindrical electrode and eight surrounding needle 
electrodes made of gold-plated steel (¢ of 1mm) in order to avoid corrosion. 

The ß dispersion, corresponding to radio frequency, is strongly dependent on the 
dielectric properties of the cell membranes and their interactions with the extra- and 
intracellular electrolytes, as the current passes through the cell membranes at these 
frequencies. It includes information referring to the behavior of the cell membrane in- 
tegrity, which changes significantly with aging (s. Fig. 6). One of the phenomena that 
contributes to the ß dispersion is the Maxwell—Wagner effect related to the interfacial 
relaxation process occurring in systems where the electric current passes at the inter- 
face between two different materials. 
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Fig. 6: Bioimpedance measurement of a beef longissimus dorsi muscle for 14 days [14]. 


In [43], a method for measuring bioimpedance non-invasively is presented. It consists 
of magneto-inductive measurement techniques, such as gradiometers. Gradiometers 
consist of a sending coil and two symmetrical receiving coils. For the detection of tis- 
sue impedance, they use the difference principle between the coil on the tissue side 
and the coil on the other side. They can therefore detect even small changes of electri- 
cal properties of a biological tissue. This high sensitivity has the drawback that they 
may have different levels of detuning, which cannot be easily adjusted. In [43], a con- 
cept for adjusting the gradiometer in a predefined initial state is presented. It uses ad- 
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ditional coils on both sides to adjust the detuning by a suitable current flow through 
the additional coils and realize an accurate measurement. 

Measuring bioimpedance of living beings in medical applications brings several 
challenges. One is the limitation of the current, which must be applied without do- 
ing any damage [46]. The medical requirements imposed by the standard IEC-60601 
claim that the current should be below the perception threshold of 0.5 mA at 5 kHz 
[44]. The perception threshold increases when the frequency increases following the 
equation [45] 


Imax = 10’f, (1) 


where the current is the root mean square value and the frequency is in Hz. 

Therefore, for in vivo applications, a current excitation is preferred to a voltage 
excitation, which can be affected by contact impedance between the electrode and the 
tissue and high currents through the tissue due to impedance changes or temperature 
dependence of the tissue conductivity. 

One further challenge is the fact that living beings change their state, e. g., due 
to breathing, moving and blood circulation. This is why one of the approaches which 
can be considered is to reduce the measurement time in order to maintain the system 
under test under quasi-static conditions. Many researchers investigated possibilities 
for using wide band signals including different frequencies in order to get a spectrum 
in a reduced measurement time [35, 47-49]. These are typically, e. g., chirp signals 
[35], multisine signals, binary chirp and binary multifrequency signals [49]. 

In [50], a novel measurement and identification approach, called the parametric- 
in-time approach, has been developed for time-varying bioimpedance systems with 
a quasi-periodic character. This allows to conduct myocardial impedance measure- 
ments in vivo and in situ. 


8 Cable fault detection and localization 


The increasing use of the wiring in vehicles, production plants, communication and 
power systems leads to a big nead for wiring faults detection and location. Experience 
in this area indicates that the majority of the problems encountered are due to cable 
errors and increased contact resistance of connectors. If the cables are long, not all 
methods are able to detect the fault position as signals may become weak and even 
concealed by noise. 

There are several methods for wire fault location and characterization, of which 
time domain reflectometry (TDR) is the most common one. In this case, a pulse signal 
is sent on the cable, which is supposed to be homogenous. If the cable has a fault, a sig- 
nal is reflected and this reflection can be measured to detect the position of the fault by 
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estimating the time of flight and the type of the reflection includes information about 
the type of change of the cable impedance and therefore about the type of the fault. 

Over the last decade, many methods, such as frequency domain reflectome- 
try (FDR), ultra-wide band-based TDR, spectrum TDR (STDR) and impedance spec- 
troscopy (IS), have been developed. They use different wide band incident signals 
and signal processing methods, but the aim was mainly to improve the accuracy of 
the fault location. 

In [6], impedance spectroscopy was introduced to identify both position and type 
of cable errors. The innovation in this approach is that the load impedance of the wire 
fault can be quantitatively extracted from its input impedance, so that the types of 
the wire faults can be exactly identified. For that, a transformation to the pseudo- 
frequency domain is applied to localize the wire fault, which corresponds to the peri- 
odicity of the spectrum represented in Fig. 7. 
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Fig. 7: Impedance spectrum of a 10-m coaxial cable with different terminating resistors. The ampli- 
tude is influenced by the terminating resistors and the periodicity of the impedance is influenced by 
the error position [6]. 


The input impedance of the wire at low frequencies is applied to identify the types of 
wire faults. The imaginary part of the cable impedance can be applied to detect the 
type of the wire fault. In the case of an open circuit, it is the capacitance of the cable. 
In the case of a short circuit, the inductance of the cable helps to detect the fault type. 


9 Conductive, capacitive and inductive sensors 


Impedance measurement has decisive benefits, as it is experimentally efficient and 
can be realized in low-cost embedded systems. Its advantages are particularly impor- 
tant in sensor technology as it delivers more information than only resistive, capaci- 
tive or inductive measurements. This is due to the fact that it uses both the real and the 
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imaginary part of the complex impedance at different frequencies. It provides the pos- 
sibility to separate effects having linearly independent dispersions and dominating in 
different frequency ranges. 

Impedance-based sensors can have a variety of designs and can serve for the ac- 
quisition of various measured variables (see Fig. 8 and Fig. 9). These include, e. g., 
position, humidity, material properties and filling level. Very often, only the real part 
is measured in the case of conductive sensors and the imaginary part in the case of 
inductive and capacitive sensors. The other part of the impedance is considered as 
losses and is in general ignored, even if it includes also information about the “device 
under test”. In some cases, the amplitude of the impedance and the measured variable 
is determined at a suitable sensitive frequency [53]. 
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Fig. 8: Electrode structures for capacitive and conductive sensors. 


The measurement of the impedance with its real and imaginary part at different fre- 
quencies expands the information considerably [54]. The information obtained from 
the sensor is representative of the measurements of a sensor array, in which each sen- 
sor is fictitiously operated at a frequency. We know from the field of multisensors that 
this increases the quality and reliability of the measurement and helps to achieve a 
better accuracy [55]. 

Impedance spectroscopy provides interesting possibilities for realizing a multi- 
functional measurement, so that one sensor measures several quantities at the same 
time. This provides interesting possibilities not only to reduce costs, but also to extend 
the information diversity from one sensor. In [13], an impedimetric sensor has been de- 
veloped for the measurement of soil moisture, having a design such as in Fig. 8 (c). It 
measures both the real and the imaginary part of the impedance at different frequen- 
cies and is thereby able to identify the type of the soil. Similar sensors measuring only 
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Fig. 9: Examples of inductive sensors for characterization of materials and non-destructive material 
testing. (a) immersion core principle, (b) Immersion core principle transformer principle, (c) inner 
wall characterization of a pipe, (d) material characterisation and testing by a reflective principle, (e) 
material characterisation and testing by a reflective difference principle, (f) material characterisation 
and testing by a transmissive principle. 
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Fig. 10: Dependency of the imaginary part of the soil impedance on soil type and soil moisture. 
(a) Investigated soil types classified by their inorganic components [13]. (b) Soil imaginary part in 
dependence of the frequency for different soil types. 


the capacitive part have massive problems with soil type influence, which leads to a 
high inaccuracy and needs therefore to be regularly calibrated to the corresponding 
soil. As can be seen in Fig. 10 (b), the changes of the imaginary part of the impedaance 
due to moisture show a different behaviour depending on the soil type. Measuring the 
capacity at a certain frequency, leads to ambiguityan is not sufficient for realizing an 
accurate measurement. In the case of the impedimetric sensor, the measurement data 
serve first to the classification of the soil type by means of data reduction via principal 
components analysis and classification via the k-nearest-neighbor method. After soil 
type recognition, the moisture measurements become easy and accurate by using a 
linear characteristic of the imaginary part of the admittance [13]. 
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Inductive sensors and especially eddy current sensors can be used in different 
manners for both distance measurement and material characterizations. The main ad- 
vantages of these sensors are a.o. their robustness [56], their process capability and the 
fact that they operate contactless and in a non-destructive manner [57]. 

Yet, if the imaginary part changes more, also the real part includes interesting 
information. A change of frequency leads to a different penetration depth of the mag- 
netic field in the target material (see Fig. 11) and the coil impedance can therefore give 
information about cracks and failures under the surface. 
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Fig. 11: Dependence of the eddy current effect on the frequency. (a) Schematic description of the 
penetration of a magnetic field in a target material in dependence of the frequency. (b) Numeri- 
cal simulation of the magnetic vector potential in a ferric target material in dependence of the fre- 
quency. 


10 Biosensors 


Impedance spectroscopy has recently gained a high popularity for biosensing of a 
wide range of substances owing to its label-free measurements. It reaches thereby a 
superb sensitivity down to the femto- and attomolar ranges [58-60]. It allows rapid 
screening of biocompatible surfaces and monitoring pathogenic bacteria, as well as 
the analysis of biological cells and tissues. 

The impedance spectroscopy transducer signal is in this case mainly based on 
the interactions between a certain biological receptor and the target species, which is 
selectively adsorbed from the solution (see Fig. 12) [61, 62]. Such interactions cause a 
change of the interfacial electron transfer kinetics between the solution species and 
the conducting electrode. The charge transfer resistance increases therefore with the 
increase in the quantity of targets bound to the receptive surface. Thereby, the selec- 
tivity, sensitivity, response time and detection limit are strongly dependent on the use 
of proper materials and catalysts [61, 63]. 
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Fig. 12: Three-electrode electrochemical sensor. (a) Sensor element. (b) Selectivity enhancement 
by using template molecules. (c) Strategies for improving sensitivity by using nanoparticles and 
catalysts. 


Because of their high affinity to biomolecules and high aspect ratio, nanomate- 
rials, such as metal nanoparticles, CNTs and graphene, are revolutionizing electro- 
chemical biosensors (see Fig. 12 (c)) [64, 65]. Nanomaterials are often used as elec- 
trode coating to increase their effective surface area and immobilize a large number of 
biomolecules such as enzymes and antibodies [65]. 


11 Future trends of impedance spectroscopy 


The impedance spectroscopy is versatile as a method and its importance is in general 
increasing, as it provides interesting possibilities for characterizing materials, systems 
and sensors. Sensors themselves today tend to the use of high-speed and low-cost elec- 
tronic circuits, allowing for increasingly more signal processing and advanced manu- 
facturing technologies [52]. Similar aspects can be observed nowadays for impedance 
spectroscopy. Measurement devices use novel technologies and provide increasingly 
efficient solutions, enabling the use of the method in a wider range of applications. 


11.1 Embedded and miniaturized solutions 


Impedance spectroscopy was regarded as an electrochemical laboratory method for 
several decades. Today, due to the excellent developments in microelectronics, a cost- 
effective electronic realization becomes more and more possible. The most available 
embedded solutions are in general in the medical field and a lot of them are in medical 
implants. This trend is now increasing towards manifold solutions for sensors and 
diagnostic devices. In most cases, the aim is not to operate a large measuring range, 
but rather to find solutions that are usually limited by the requirements of a particular 
application. 
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11.2 Fast impedance spectroscopy and tomography 


One of the important conditions for impedance spectroscopy is that the system does 
not change during measurement of the whole impedance spectrum. If the system has 
a higher dynamic, impedance measurement should be carried out within a short time 
to realize quasi-stable conditions. This is why reducing the measurement time is nec- 
essary to widen the scope of applications of impedance spectroscopy to more dynamic 
applications. 

Many examples can be seen in impedance tomography. In this case, measure- 
ments with different electrode combinations contribute to the reconstruction of one 
image across the device under test. All these measurements need to be adjusted to the 
same system state. In order to fulfil this condition, in general only one or two rela- 
tive high frequencies are measured. The development of fast impedance spectroscopy 
methods would enable tomography methods to better use the spectral information by 
involving more measurement frequencies. 


11.3 Impedimetric sensors 


Nowadays, in the sensor field there is an interesting trend towards multifunctional 
sensors, which are able to measure more than one quantity at the same time [51, 52]. 
Impedimetric sensors can provide more information as the impedance is measured 
and in a certain frequency range. The comprehensive data generated at one value of 
the measured quantity provide interesting possibilities to reduce influencing effects, 
improve the reliability of the information and even realize the measurement of more 
quantities at the same time. In this way, many classical capacitive, resistive or induc- 
tive sensors can reach a higher performance by using impedance spectroscopy. 
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Part I: Electrochemical and physical phenomena 


Abdulkadir Sanli, Abderrahmane Benchirouf, Christian Müller, and 
Olfa Kanoun 


Study of the humidity effect on the electrical 
impedance of MWCNT epoxy nanocomposites 


Abstract: Owing to the excellent physical properties of carbon nanotubes (CNTs), in- 
corporation of CNTs into polymer matrices would make the nanocomposites a poten- 
tial candidate for smart sensing applications. However, environmental factors such 
as humidity might influence the nanocomposite properties and performance. In this 
work, the effect of humidity on the electrical properties of multiwalled CNT (MWCNT) 
epoxy nanocomposites is investigated by electrochemical impedance spectroscopy in 
the frequency range from 40 Hz to 110 MHz. The MWCNT epoxy nanocomposites with 
different MWCNT concentrations up to 1% wt were fabricated by the direct mixing 
method and the composites were deposited on the glass epoxy substrate by the sten- 
cil printing technique. From the impedance results, it was observed that the overall 
impedance of the nanocomposites tends to increase with the increase of the humidity 
due to an increase in tunneling barriers induced by swelling of the polymer layer be- 
tween the CNT junctions within the nanocomposite matrix. The humidity effect tends 
to decrease with the increase of MWCNT concentration due to a decrease of the tunnel- 
ing effect at high MWCNT concentrations. Furthermore, the proposed nanocomposites 
at a 0.5% wt concentration show a humidity sensitivity of about 0.006 % relative hu- 
midity (RH) and 0.004 % RH and linear responses with linear correlations (R?) of 0.98 
and 0.93 for 100 Hz and 1kHz, respectively. Moreover, the sensors were covered with 
latex polymer to minimize the humidity effects. It was found that with encapsulation 
of the nanocomposites, the relative impedance change of the nanocomposite is de- 
creased from 3.76 % to 1.15 % at a 0.5 % wt MWCNT concentration. 


Keywords: carbon nanotubes, epoxy resin, nanocomposites, humidity sensing, en- 
capsulation, electrochemical impedance spectroscopy 


1 Introduction 


Since their discovery of by lijima in 1991 [1], carbon nanotubes (CNTs) have been used 
extensively in most fields of science and engineering owing to their excellent mechan- 
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ical, thermal and electrical properties [2-5]. In comparison to other fillers, such as 
nanoparticles or carbon black, CNTs are considered as an effective filler owing to their 
high aspect ratios (= 100-1000), which result in the formation of conductive paths 
at very low CNT concentrations [6, 7]. Thus, integration of CNTs into a polymer ma- 
trix would make the composites a potential candidate for smart sensing applications, 
including chemical, pressure, gas, flow and mechanical sensors [8-11]. The environ- 
mental factors, e. g., humidity, can influence the overall electrical conductivity as well 
as sensor performance of a polymer nanocomposite for the sensor applications such 
as strain and pressure sensing, in which the compensation of environmental effects is 
required. Therefore, it is important to quantify the influence of environmental factors 
on the electrical conductivity of CNT-based nanocomposites. To reveal the conduction 
mechanism of the CNT polymer nanocomposites under the humidity effects, several 
researchers have proposed various explanations [12-17]. For instance, Yoo et al. [13] 
suggested that resistivity changes are caused by polymer swelling below the percola- 
tion threshold, whereas above the percolation threshold the main reason in resistivity 
change is the charge transfer from water molecules to CNTs. On the other hand, Yu 
etal. [14] suggested that the change in the conductivity is due to the formation of a 
weak bond between H and C atoms from water and CNTs induced by water uptake. 
Asrafi et al. [16] have found that a change of humidity from 5 % relative humidity (RH) 
to 95% RH can cause an 80 % reduction in conductivity of CNT epoxy nanocompos- 
ites due to moisture ingress. However, to the best of our knowledge, there is no study 
addressing the effects of the CNT concentration and polymer coating on the electrical 
properties of MWCNT epoxy nanocomposites under humidity effects. 

Therefore, this work focuses on the quantitative investigation of the effect of 
humidity on the electrical conductivity of a free standing and encapsulated MWCNT 
epoxy nanocomposite by electrochemical impedance spectroscopy for different 
MWCNT concentrations. 


2 Experimental procedures 


2.1 Measurement setup for humidity investigations 


Prior to placing the samples inside the salt container, a precise humidity sensor is used 
to identify the values of humidity and temperature within the salt container. After the 
sample has been placed inside the container, impedance measurements were taken 
after ca. 45 min for the stabilization of the humidity. After that, the impedance data 
were recorded by utilizing the Agilent 4294A precision impedance analyzer in the fre- 
quency range of 40 Hz to 110 MHz. The measurement setup for humidity analysis of 
MWCNT epoxy nanocomposites is represented in Fig. 1. 
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Fig. 1: Schematic illustration of (a) humidity measurement setup for humidity investigations of 
MWCNT epoxy nanocomposites placed inside a glass container and (b) encapsulation of samples 
with latex solution. 


3 Results and discussion 


3.1 Effect of humidity on the impedance response 


It can be seen from both the Nyquist plot and the Bode plot (see Fig.2) that the 
impedance of the MWCNT epoxy nanocomposites tend to increase with the increase 
of RH. The Nyquist plot shows quasi-semi-circle behavior for all concentrations and 
humidity responses. In order to estimate the complex transport mechanism in MWCNT 
epoxy nanocomposites under the humidity effect, a 3D R-C model can be proposed to 
reveal the effect of the MWCNT concentration and the interaction within the nanocom- 
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Fig. 2: (a) Nyquist and (b) Bode plot of non-capsulated MWCNT epoxy nanocomposites under RH. 
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Fig. 3: Illustration of the working mechanism and corresponding equivalent circuit model of an 
MWCNT epoxy nanocomposite. 


posite matrix, shown in Fig. 3. Here, each nanotube (or nanotube agglomerate) has 
a resistance Rcyr; and the polymer medium can be modeled as Runne; and constant 
phase element (CPE), respectively. 


3.1.1 Humidity sensing mechanism 


From the accumulated studies [18-23], the effect of humidity on the CNT polymer 
nanocomposites can be attributed to the following working mechanisms: (i) the in- 
crease of intrinsic resistance of CNTs due to adsorption of water molecules [20], (ii) the 
charge transfer between the water molecules and CNT modules and the Fermi level 
that causes the resistance change across CNT-CNT junctions [21, 22] and (iii) the 
change in inter-tube distance induced by polymer swelling under the humidity [23]. 
For the case of CNT epoxy nanocomposites, the main reason behind the increase 
of impedance under humidity is the polymer swelling. Many researchers have sug- 
gested that the swelling of the polymer layer between CNT-CNT junctions within the 
nanocomposite matrix due to humidity causes the increase of the tunneling barriers 
[18, 23]. As the polymer expands, the CNTs in the network within the polymer matrix 
tend to move apart from each other as if the nanocomposite deformed under the ten- 
sion, which causes the alteration of the complex impedance of the nanocomposites. 
Thus, the electrical impedance increases under humidity as shown in Fig. 4. There- 
fore, the effect of intrinsic resistance change of CNTs under humidity can be neglected. 

Here, it was also found that the humidity sensitivity of the nanocomposite de- 
creases with increasing the MWCNT concentration due to a decrease in the effect of the 
tunneling resistance. When the CNT concentration is high, the MWCNTs locate very 
close within the nanocomposite matrix. Therefore, the inter-tube distances at high 
concentrations are less effected by the polymer swelling than those near the perco- 
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lation threshold [19]. While for the samples of a 0.5 % wt MWCNT concentration the 
relative resistance change is 3.76%, this value decreases to 0.72% for 1% wt MWC- 
NTs. Moreover, as illustrated in Fig.5, it was seen that the operation frequency has 
a great effect on the impedance in MWCNT epoxy nanocomposites under RH. It has 
been observed that the proposed nanocomposites of a 0.5 wt % MWCNT concentra- 
tion show humidity sensitivities of about 0.006 % RH and 0.004 % RH with linear re- 
sponses (R°) of 0.98 and 0.93 for 100 Hz and 1 kHz, respectively. The proposed MWCNT 
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Fig. 5: Impedance response of MWCNT epoxy nanocomposite at a 0.5 % wt MWCNT concentration 
under RH for different frequencies up to 100 MHz. Inset image shows the corresponding impedance 
response ofthe MWCNT epoxy nanocomposites up to 1kHz. 
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epoxy nanocomposite can be also used as humidity sensor as it was found to be more 
sensitive than many other nanocomposites [13-24]. 


3.2 Effect of the encapsulation on the impedance response 


In order to minimize the effect of the humidity on the electrical properties of MWCNT 
epoxy nanocomposites, samples were encapsulated with latex polymer (see Fig. 6b). 
After that, the samples underwent the same humidity measurements again. From 
Fig. 6, it can be seen that the resistance value of the composites slightly decreases 
after the encapsulation due to the effect of the latex polymer, which acts as an ad- 
ditional parallel resistance to the composite. By extracting the resistance values of 
the encapsulated nanocomposites, it can be seen that the latex encapsulation signif- 
icantly minimizes the effect of the humidity. While the relative resistance change of 
the nanocomposites for 0.5 % wt decreases from 3.76 % to 1.15 %, this value decreases 
from 0.72% to 0.52% for a 1% wt MWCNT concentration, as shown in Fig. 7. 
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4 Conclusions and outlook 


In this study, we investigated the effect of humidity on the impedance properties 
of MWCNT epoxy nanocomposites for different MWCNT concentrations. From the 
impedance results, we found that the impedance of the nanocomposites increases 
with the increase of the humidity due to an increase in tunneling barriers induced by 
swelling of the polymer layer. Additionally, the humidity effect was found to decrease 
with the increase of MWCNT concentration due to the decrease of the tunneling effect 
at high MWCNT concentrations. Furthermore, the proposed nanocomposites at a 0.5 
wt% MWCNT concentration show a humidity sensitivity of about 0.006 % RH and a 
linear response (R?) of 0.98. Moreover, by encapsulation of the nanocomposites the 
relative impedance change of the nanocomposite decreases from 3.76 % to 1.15% ata 
0.5 % wt MWCNT concentration. For the future work, it would be interesting to explore 
the effect of the different encapsulation methods on the electrical and piezoresistive 
properties of the MWCNT epoxy nanocomposites. 
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Lukasz Macioszek and Ryszard Rybski 

Temperature influence on impedance 

of premium summer diesel fuel measured 
with the use of impedance spectroscopy 


Abstract: In this chapter, an attempt is made to experimentally verify temperature in- 
fluence on the evaluation of diesel fuel impedance with the use of electrochemical 
impedance spectroscopy. Previous test findings showed that this method can be used 
to detect contamination, such as water, in constant temperature. Tests presented in 
this chapter show that temperature itself can change electrical properties of diesel fuel 
similarly to water. Moreover, tested premium diesel, as a function of temperature, be- 
haved oppositely from what hydrocarbons normally do. This makes temperature and 
additives present in fuel very important factors, which must be taken into considera- 
tion in the evaluation of any diesel properties using impedance spectroscopy. 


Keywords: diesel, fuel, impedance spectroscopy, temperature 


1 Introduction 


Diesel is a very common fuel used to drive engines of vehicles, locomotives, ships and 
even electricity generation plants. It is a result of crude oil distillation and consists of 
several types of hydrocarbons (about 75 % saturated, 25 % aromatic), as well as vari- 
ous kinds of additives and enhancers added during or after production [1]. A few per- 
cent of fatty acid methyl ester (FAME) is also added in many countries to each fuel 
commercially available, for example, in Poland there is 7 % (V/V) of biocomponents. 
Diesel analysis is very difficult in terms of interpreting obtained results. The en- 
countered problems can be divided into two basic kinds. The first is that only very few 
research institutes are studying electrical properties of fuel, including its impedance. 
They have to be studied at refinery laboratories but the results are not publicly avail- 
able. In sucha case, information can only be obtained via private communication with 
people working in those facilities. The publicly available publications are mainly fo- 
cused on fossil diesel rather than modern commercially sold fuel. The second kind of 
problems for a researcher is connected to the fact that the precise chemical compo- 
sition of diesel is never the same. Responsible for this is the distillation, which, like 
any process, is never performed with 100 % efficiency. Additionally, the crude oil it- 
self is not always the same. There are also a number of additives that make fuel safe to 
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transport, improve its lubricity, reduce the risk of freezing in winter, etc. [2]. Used ad- 
ditives are always producers’ secrets and may change at any time, even though there 
is an obvious need for a maintenance of the quality of diesel fuel that is being sold. 
Within the European Union the EN 590 standard is valid. It contains ranges of phys- 
ical properties’ permissible values and lists the corresponding test methods to verify 
them [3]. 

Premium-type diesel fuels are declared by producers as containing engine clean- 
ing substances [4], improving lubricity [5], having a higher cetane number [6, 7] and 
being sulphur-free [7]. Some are advertised as being capable of protecting engines 
against corrosion and increasing the driving range in comparison to regular fuel [5]. 
Winter premium diesel should enable proper functioning of engines in lower ambient 
temperatures than regular diesel. Due to this fact it has to contain additional sub- 
stances. Aforementioned characteristics make the analysis of premium-type fuel even 
more complicated. Nevertheless, premium diesel is believed to have more consistent 
content and this is the reason why it was chosen to be examined. 

The improved efficiency and performance of diesel engines imply an increase in 
their sensitivity to all kinds of contaminants in the fuel. Very high pressures, exceed- 
ing 300 MPa, are featured by modern direct injection systems (called common rail). 
Such values would not be possible to achieve without rigorous fuel parameters that en- 
sure the absence of substances aggressive towards metal alloys, of which engines are 
composed. Particularly newer engines’ fuel rails are especially sensitive to water con- 
tent in the fuel; therefore introducing additional contamination measurement meth- 
ods seems to be necessary. Only 200 mg/kg of water is allowed by the aforementioned 
standard EN 590. The coulometric Karl Fischer titration method is used for this deter- 
mination, as it allows obtaining precise results in the range of 0.003-0.100 % (m/m). 
The solubility limit of water in diesel is reported to be very low (100 ppm at 40°C [8]) 
but various kinds of additives can make it significantly higher [2]. Although amounts 
found in fuel can be considered as relatively low, water should still be treated as a dan- 
gerous contamination. Coulometric titration is a very precise method, though it cannot 
be used in situ as it needs to be performed in the laboratory and requires preparation 
of test samples. One of the various methods proposed for the evaluation of diesel fuel 
properties is impedance spectroscopy [9, 10] but it has not yet been used for water con- 
tent assessment in any device. Such device, if designed as portable, could be used for 
example at the gas station to perform a quick preliminary measurement and indicate 
whether the examined fuel may be contaminated and should be sent to the laboratory 
for further, more precise studies. 

In order to propose a complete method of measuring water content as well as any 
other contaminant in diesel fuel with the use of impedance spectroscopy, evaluation 
of the influence of the temperature on diesel impedance is required. Pure diesel fuel 
can be treated as dielectric [9] but each additive makes its dielectric properties more 
complex. Hydrocarbons, of which diesel is composed, tend to lower their resistivity 
with the raise of temperature. Since around 2006 the production of ultra-low sulphur 
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diesel (ULSD) has been growing and this kind of fuel dominates the market nowadays. 
Processes used to remove sulphur also significantly decrease fuel conductivity [2]. On 
the other hand, transportation standards require the conductivity of diesel to be above 
50 pS/m [11] to allow static charge to dissipate. This is the reason for adding substances 
that will not noticeably change all fuel properties except conductivity. What kind of 
additives are added is a secret of the fuel producer. There are many substances avail- 
able but their composition is also classified in general. It is known that additives are 
added in very small proportions, of a few tens to hundreds parts per million, and they 
can raise diesel conductivity tens of times [2]. This implies that additives may alter or 
overshadow the behavior of hydrocarbons impedance as a function of temperature, 
thus measurements of diesel impedance at different temperatures are very important 
to find their correlation. Once this correlation is found, it should be included in the 
diesel electric equivalent circuit in order to model fuel behavior at different tempera- 
tures. Such a model is necessary to assess fuel contamination with the use ofa portable 
device based on impedance spectroscopy in the future. 


2 Material and methods 


The main objective of this study is to experimentally check the temperature influ- 
ence on the impedance of premium summer diesel fuel measured with the use of 
impedance spectroscopy. Fuel samples were obtained from a gas station and were ex- 
amined at different temperatures in a temperature controlling container. Impedance 
values were the same during measurement series with increasing and decreasing tem- 
perature. Once impedance of samples was measured as a function of frequency, the 
results were used to fit values of diesel equivalent circuit elements. 


2.1 Diesel fuel samples 


Diesel samples of 500 ml capacity were examined in glass beakers (shown in Fig. 1 (a)). 
This rather high capacity was used in previous experiments to prevent unplanned, 
spontaneous changes of temperature during the experiment. Although such situa- 
tions did not take place, it was used again to preserve maximum possible repeatability 
of experiments. Temperature was monitored all the time and was approximately con- 
stant (+0.1°C) during each experiment. Measurements were performed at 6.9°C, 9.3°C, 
11.1°C, 14.5°C, 16.4°C, 18.4°C, 23°C and 26°C (+0.1°C). These values were chosen arbi- 
trarily as preliminary ones just to cover the basic range. With the equipment the used 
temperature could be more precisely measured than be set; nevertheless, the measur- 
ing conditions were good enough to perform planned studies on the behavior of diesel 
impedance at different temperatures. 
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(a) 


Fig. 1: (a) Sample with immersed electrodes. (b) Electrodes used in the experiment. 


2.2 Electrodes 


Measurements were performed using two circular parallel electrodes made of gold- 
coated copper with a surface area of 25 cm? each, that were spaced by 2mm. Dur- 
ing experiments they were immersed in diesel samples. Fig. 1 (b) illustrates both elec- 
trodes and four supports providing the planned distance. As can be seen, the actual 
electrodes were in the form of a copper layer (35 pm) on a laminate, the same that is 
usually used to build single-layer electronic circuits. Leaving gold-coated electrodes 
immersed in diesel for long periods of time (up to 14 days) did not reveal any chemical 
reactions that would take place on their surface. This allowed for the assumption that 
cleaning electrodes after each measurement is not necessary to maintain repeatabil- 
ity, which was confirmed by later obtained results. 


2.3 Impedance spectroscopy 


Diesel fuel sample impedance was measured in the frequency range of 0.1-500 Hz with 
a 700-mV RMS voltage. A laboratory EG&G/Princeton Applied Research electrochem- 
ical impedance spectroscopy system was used, consisting of the 263A Potentiostat- 
Galvanostat, 5210 Dual Phase Lock-In Amplifier and PowerSINE software. The system 
was calibrated according to the manufacturer’s recommendations and self-calibration 
before each measurement was also performed. Each final measuring result was an av- 
erage of four separate impedance measurements that were performed by the system. 
This number was chosen arbitrarily and was a compromise between measurement 
quality and time; each experiment took about twenty minutes. Measured impedance 
in the form of 


Z* (w) = Z'(w) - jZ" (w) (1) 
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was used to fit equivalent circuit element values [12]. The circuit was simple and con- 
tained a resistor with parallel constant phase element. Such circuit was already used 
elsewhere to describe diesel properties [9, 10]. Impedance of such circuit can be de- 
scribed by 


R 


2) = TERGO 


(2) 


where Q is the admittance 1/|Z| at w = 1rad/s and real n satisfies O < n < 1. 


3 Results and discussion 


Diesel fuel samples were examined at temperatures listed in Section 2.1. Stray effects 
can be omitted as several measurement repeats did not reveal values exceeding the 
measurement error of the system used. Fig. 2 presents the Nyquist plot with the mea- 
sured impedance of eight diesel samples. It can be seen that measured values form 
single semi-circles. The right ends of the semi-circles are a bit raised which suggests 
that within this low-frequency range there may be valuable information. However, fur- 
ther measurements at frequencies down to 1mHz did not reveal second semi-circles 
and hence they were not included in the chapter. Although the used frequency range 
of 0.1-500 Hz may seem to be narrow, it was observed in the evaluation of another 
dielectric, namely transformer oil, that using a wider frequency range did not reveal 
more information [13]. The raise may be associated with diffusion phenomena occur- 
ring in the dielectric when a metal electrode is used [14] and it is not considered as 
a primary and dominant phenomenon. The centers of the semi-circles are a bit de- 
pressed, which suggests that the calculated capacitance would vary with frequency. It 
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Fig. 2: Nyquist plot of examined diesel fuel samples. 
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confirms that choosing a simpler equivalent circuit consisting of a resistor and ideal 
capacitor would lead to greater fitting errors. 

Note that with the raise of the temperature, the diesel resistivity increases — semi- 
circles have larger diameters. The real part of the measured impedance at 0.1 Hz re- 
vealed values in the range of 1.35 GQ at 6.9°C to 2.25 GQ at 26°C. Increasing resistivity 
with the raise of temperature is normally encountered in metals, and hydrocarbons 
should behave in the opposite way; their resistivity decreases with the increasing tem- 
perature. This may suggest that additives present in fuel responsible for increasing its 
conductivity play a leading role in impedance measurements. Electrode material and 
the fact that we used premium-type fuel were not the reasons for these unexpected 
results, as experiments with the use of copper electrodes and regular type diesel re- 
vealed the same resistivity behavior [15]. 

Since the diameters of semi-circles from the Nyquist plot were the only visible 
changes to their shape according to different temperatures, it was decided that overall 
resistance of the diesel sample can be used to propose the calibration curve that would 
model the temperature influence on the diesel’s electrical properties. As can be seen 
in Fig. 3, the correlation of the samples’ calculated resistance and temperature seems 
to be linear, with a quite good coefficient of determination (r? = 0.9774). 
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Fig. 3: Temperature dependence of fitted R values of studied diesel samples. 


Table 1 contains fitted values of equivalent circuit elements and the corresponding 
x modeling errors. Relative errors of fitting individual values did not exceed 3.6 % in 
every examined sample. Note that the most changing value is the parallel resistivity of 
the equivalent circuit with 1.32 GO at 6.9°C and 2.229 GQ at 26°C. Even though values of 
Qare also varying, it is hard to assess its trend, and therefore their differences are not 
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Table 1: Fitted values of equivalent circuit elements and corresponding modeling errors. 


Temp. CC) R(GQ) Q(pSs”) n Fitting (x?) 


6.9 1.320 39.22 0.98216 0.00148 
9.3 1.395 39.76 0.98600 0.00074 
11.1 1.423 39.92 0.98282 0.00083 
14.5 1.606 45.21 0.97990 0.00074 
16.4 1.720 45.34 0.98268 0.00069 
18.4 1.890 47.94 0.97489 0.00114 
23 1.979 44.69 0.98356 0.00052 
26 2.229 42.37 0.98278 0.00062 


so significant as R values (up to about 20 % raise compared to 69 %) in the investigated 
temperature range. 

A very important and unexpected result of the performed experiments is that the 
examined diesel fuel had a positive temperature coefħcient of resistance. The experi- 
ment was repeated with the use of copper electrodes and regular (not premium) diesel. 
Although the fitted values of the equivalent circuit were a bit different, the examined 
fuel still had a positive temperature coefficient of resistance. Hydrocarbons do not 
have such property, thus some other substance present in the fuel caused this phe- 
nomenon. The exact chemical compositions of aforementioned additives that are used 
to increase the fuel’s conductivity are unknown, although there is a limited number 
of possibilities. It is only known that these enhancers are added to fuel in proportions 
up to one hundred parts per million [2]. A conductivity increase of raw ULSD from 
about 11 pS/m to required 50 pS/m could be achieved with the addition of conductive 
polymers during the production. Some of these polymers can have very high conduc- 
tivity, at the same level as good metal conductors. Thus, it is possible that even a very 
low proportion of polymers in fuel (below 50 ppm) would change conductivity of the 
mixture vitally. Moreover, some of the conductive polymers have a positive tempera- 
ture coefficient of resistance, just like metals do. Given properties incline us to assume 
that it is possible that some conductive polymers could be present in the examined 
premium diesel fuel. Confirmation or rebuttal of this thesis will require multiple ex- 
periments with the use of different techniques, as diesel has a very heterogeneous 
composition. Nevertheless, an explanation of what causes diesel fuel to increase its 
resistivity with the raise of temperature seems to be indispensable. 

The encountered behavior of diesel’s impedance as a function of temperature does 
not mean that further studies of its impedance are pointless. Correlation between tem- 
perature and diesel’s resistance in the examined range seems to be linear. This means 
that its influence is easy to predict and it definitely should be included in more detailed 
equivalent circuits that would model diesel’s electric properties. Evaluation of differ- 
ent diesel fuel types should answer the question whether all types of diesel behave 
the same way, such as summer, winter, premium and regular diesel. If the tempera- 
ture coefficients are, every new study of diesel’s electrical properties will be valid. 
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4 Conclusion 


In this chapter, the influence oftemperature on impedance of premium summer diesel 
fuel measured with the use of impedance spectroscopy is presented. Experimental re- 
sults are discussed and analyzed. A calibration curve as well as fitted values of the 
equivalent circuit are presented. They together show that temperature has a great im- 
pact on the impedance of diesel fuel. 

Evaluation of the influence of temperature on diesel impedance is required in or- 
der to propose a complete method for measuring diesel contamination with the use 
of impedance spectroscopy. Given the limitations discussed in Section 3, we can still 
state that impedance spectroscopy of diesel fuel is reasonable and compares favorably 
with methods that require preparation of the samples and thus can only be performed 
in the laboratory. Analysis of the obtained results may be very difficult, as the exact 
composition of diesel is unknown; nevertheless, impedance spectroscopy can be used 
as a comparative method. As mentioned before, there does not exist nor will there ever 
bea reference fuel. However, changing only one property of any diesel and measuring 
its impedance before and after that change will always be informative. The proposed 
methodology may not reach the same high level of achieved accuracy as laboratory 
methods used nowadays; nevertheless, accuracy is not the most important factor in 
each case. Future studies about diesel impedance behavior in a wider temperature 
range with and without the addition of other substances should answer the question 
whether there is a possibility to indicate the harmfulness of fuel based on its measured 
impedance. 

Advantages of using impedance spectroscopy as a measuring technique make it 
possible to build a device that could theoretically work in situ. After relatively fast mea- 
surements and proper calculations it could warn a user that contamination was found 
and the examined diesel fuel can be potentially dangerous to his car. If such a device 
would be designed to serve an inspector at a gas station, then it would indicate that 
the examined fuel is somehow different and it definitely should be sent to the labora- 
tory for the whole evaluation of its properties. Plausible relatively low measurement 
precision of the proposed device could cause sporadic fake alarms. However, we be- 
lieve that fake alarms are still far better than potentially expensive repairs of the fuel 
injection system in a modern car. 
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Anup Barai, Andrew McGorden, and Paul Jennings 
A study of the effect of short-term relaxation 
on the EIS test technique for EV battery cells 


Abstract: Electrochemical impedance spectroscopy (EIS) is a powerful tool employed 
to investigate the fundamental electrochemical parameters of a Li-ion battery cell. Its 
effectiveness has established it as a core method to study electrochemical behavior 
of batteries in laboratory environments. An industry-wide interest currently exists to 
employ EIS as part of battery management systems (BMSs) for state of charge, inter- 
nal temperature and state of health measurements. The major issue of employing EIS 
for online measurement as part of BMS is the reliability. In this work, it is shown that 
not all the parameters extracted by the EIS technique are reliable when measurement 
is performed shortly after charge/discharge. While pure Ohmic resistance was found 
to be consistent, the change transfer resistance and double layer capacitance are not 
reproducible. Therefore, only R, can be reliably employed for the fast online measure- 
ment by BMS. 


Keywords: electrochemical impedance spectroscopy, battery cells 


1 Introduction 


Intrinsic properties of lithium-ion (Li-ion) batteries, such as high energy, high power 
density, long cycle life and low self-discharge rate, have enabled their use in electric 
vehicle (EV) applications. Equivalent circuit models (ECMs) are commonly used by 
the battery management system (BMS) of an EV to estimate battery performance in 
real operating scenarios. In addition, investigation of the electrochemical processes, 
i.e., Li-ion diffusion in the electrolyte, migration through the solid electrolyte inter- 
face (SEI) layer, charge transfer through the electrode/electrolyte interface and solid 
state diffusion in the bulk of the active material can be achieved by using electrochem- 
ical impedance spectroscopy (EIS). As each process has a different time constant, time 
domain or frequency of excitation [1, 2], EIS allows for the separation of most of these 
processes [3]. Therefore, EIS always interests academic and industrial researchers as a 
diagnostic and research tool. A standard EIS test within a frequency range from 10 kHz 
to 10 mHz and 10 frequency points per decade of frequency, as commonly employed 
in the literature, requires around 10 minutes. However, if the number of frequency 
points is reduced, the test duration can be aggressively reduced to seconds. This po- 
tential reduction of test duration has allowed EIS to be incorporated as an integral part 
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ofBMS to estimate the state of charge (SoC) and state of health (SoH) in EVs [4-6]. For 
such onboard EIS measurement systems (incorporated as part of the BMS), a short 
measurement time is a key requirement. However, fast measurement implies that the 
cells will be measured very shortly after cycling. Recently the cell relaxation period 
has been identified as a source of impedance variation [7, 8]. When the charge or dis- 
charge process of the cell is stopped, the polarization of the cell ceases and as a result, 
the cell impedance and voltage of the cell change, which is referred to as relaxation of 
the cell and which is mainly caused by relaxation of ions within the double layer, elec- 
trolyte ionic concentration gradient redistribution and solid state diffusion of lithium 
atoms within the bulk of the materials. This lead to change of cell impedance in a 
couple of hours after the charge/discharge event. In [8] it has been shown that a mea- 
surable change of impedance will continue until 4 hours after a discharge event. This 
research indicates potential hurdles to get repeatable and reproducible results from 
EIS tests when employed shortly after cycling. Currently there is no research present 
in the literature showing how reliable and reproducible the EIS tests are when per- 
formed immediately after charge/discharge. This research focuses on the impedance 
variation of the cell within a short period (0-10 min) after removing an electrical load 
from the battery cell. In Section 2 the experimental procedure is outlined. Results and 
discussion of the results are presented in Section 3. Finally, the findings of this study 
are summarized in Section 4. 


2 Experimental methods 


Five commercially available Li-ion cells of three different chemical natures were used 
as part of this study. These cells were selected based on their power and energy rat- 
ing. Table 1 summarizes the chemistry, capacity and format of each cell. These cells, 
representing different groups of Li-ion batteries, i. e., chemistry, power capability and 
form factors, are currently used in different automotive applications. 

In order to isolate the effects of cell relaxation, the temperature and SoC were kept 
constant for each test; tests were carried out at 25°C and 50% SoC. To precisely ad- 
just the cell SoC to 50 %, the cells were discharged at a 1C rate to the manufacturer’s 


Table 1: Li-ion cell electrical properties. 


Cellnumber Chemistry, cathode/anode Capacity (Ah) Format Celltype 


1  LiNiMnCoO2(NCM)/LiC6 40 Pouch Power cell 
2  LiFePO4(LFP)/LiC6 20 Pouch Power cell 
3 NCM/LIC6 2.2 Cylindrical Energy cell 
4 NCM/Li4Ti4012(LTO) 13 Pouch Power cell 
5 NCM/LIC6 3.4 Cylindrical Energy cell 
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recommended cut-off voltage, subsequently allowed to rest for 3 hours before being 
fully recharged according to the manufacturer’s recommended charge protocol, using 
a commercial Li-ion cell cycler. At the end of charging, the cells were allowed to rest 
for 3 hours prior to being discharged to 50 % SoC at a 1C rate. 

To study the short-term effect of relaxation, measurements were performed in gal- 
vanostatic mode using a commercial EIS system. Impedance measurements between 
500 mHz and 10 kHz with five frequency points per decade were taken. This frequency 
range and this number of measurement points were selected as a balance of test dura- 
tion and number of measurement points. In this setup a measurement will take around 
10 seconds. Recording of the EIS data started as soon as the SoC adjustment was com- 
pleted (i. e., immediately after the 1C discharge to 50 % SoC), and every minute there- 
after. 


3 Results 


The shape of the Nyquist plots at high frequency is similar to that found in the lit- 
erature. However, at low frequency, cells 1, 2 and 4 show a spiral shape (Fig. 1 (a), 
(b) and (d)), which has not been reported in previous related literature. The spiral 
shape in the Nyquist plot starts to shrink from the second test done at 1 min after SoC 
adjustment, and finally disappears 10 min after SoC adjustment when the EIS plot is 
similar to the Nyquist plots found in the literature [5-8]. The change of shape of the 
Nyquist plots from 1 min to 10 min for cells 1, 2 and 4 happens at a different rate. Cell 3 
and 5 (Fig. 1 (c) and (e)) show this spinal shape only for the very first test, which was 
performed just after the discharge was stopped. Afterward the shape is similar to the 
Nyquist plots found in the literature. To investigate the experimental setup as a prob- 
able cause of the shape, the experiment has been repeated with different conductors, 
which connect the EIS equipment with the test cell, and even with different sets of EIS 
equipment; however, in all the instances similar results were found. The spiral shape 
at the low-frequency end could be related to the measurement time. For a stable EIS 
measurement, it is essential that there is no change of the cell during measurement 
at a particular frequency. At the lowest frequency employed (500 mHz), it takes 2 sec- 
onds to complete the measurement. After discharge was stopped, polarization of the 
double layer was ceased and as a result cell capacitance and thus total impedance 
change; this change slows down with the relaxation period. After 10 min of relaxation 
the change in the 2-second timescale is probably not significant; however, immedi- 
ately after discharge it probably is significant, which leads to the spiral shape. With the 
difference of chemistry, size and form factor, the impedance and frequency response 
of each cell is different, which could be the root cause of the difference of the spiral 
shapes. However, for all five cells measured, values of pure Ohmic resistance (R,) fall 
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Fig. 1: Nyquist plots of different cells from 0-10 min. 


(b) Cell 2 


0.0012 


(d) Cell 4) 


rer, 
> N 
NN 
RL N 
SEAN N 
ee N 
Am 
Ya “ 
Sale 
0.045 0.05 
---Omin --1min == 2min 
— 3min — 5min — 10min 
Zu (0) 
(e) 


0.0016 


-.-0 min 


— 3min 


Zu (0) 


(d) 


„wi min 
=: 


— 5min 


0.002 


—1imin 


— 5 min 


== 2 min 


— 10 min 


0.0024 


— 2 min 


— 10 min 


within experimental error bounds, and therefore they can be considered constant, re- 
gardless of relaxation time. This is because R, was found at a frequency higher than 
1 kHz for all the cells. Therefore, the measurement time was 1 ms or shorter and the 
change of the cell impedance within that short period was not significant enough to 
affect the measurement. Therefore, R, could be used as online measurement parame- 
ter. In contrast to R,, the charge transfer resistance and double layer capacitance will 
be inconsistent within 0-10 minutes after cycling, thus should not be employed for 
fast measurement. 
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4 Conclusion 


In this research, the variation of cell impedance as a function of relaxation period af- 
ter removing an electrical load has been investigated for five different commercially 
available Li-ion cells. It is concluded that the pure Ohmic resistance R, of the cell is 
independent of the relaxation processes occurring within the cell for the test duration. 
Therefore, R, could be used for direct comparisons, irrespective of the relaxation pe- 
riod used for the test. The charge transfer resistance and double layer capacitance were 
found to be inconsistent within 0-10 minutes of relaxation. Therefore, when these two 
parameters are used as a fast online measurement parameter, they will produce unre- 
liable results, which could affect SoC and SoH estimation. This new knowledge shows 
the possible inconsistency of SoC and SoH results via an EIS measurement system in- 
tegrated as part of BMS, which has been proposed by researchers previously. 
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Olfa Kanoun 


Humidity sensitivity investigation of reduced 
graphene oxide by impedance spectroscopy 


Abstract: In this chapter we report about graphene oxide (GO) and reduced GO (rGO) 
humidity sensitivity. The GO was reduced by thermal annealing to form rGO, which has 
increased its conductivity. The effects of the reduction temperature were investigated 
for both the electrical and humidity sensitivity properties. The electrical properties 
were studied by electrochemical impedance spectroscopy, which is in particular of 
high importance to reveal the conduction mechanism of humidity sensitivity. In addi- 
tion, the morphology and optical properties of rGO films were investigated by atomic 
force microscopy and UV-Vis-NIR spectroscopy. The reduction temperature plays an 
important rule for the rGO humidity sensitivity behavior not only by change in film 
conductivity but also by the transport mechanism. 


Keywords: graphene oxide, reduced graphene oxide, humidity sensor, impedance 
modeling 


1 Introduction 


Graphene and graphene oxide (GO) as carbon-based material with their well-known, 
outstanding properties have attracted a lot of attention in diverse fields, such as elec- 
tronic devices like transistors and diodes [1], energy storage [2] and sensors [3, 4]. GO 
is considered as a very promising material due to its water dispersibility enabling so- 
lution processing [5], as well as large-scale production and easy processing. GO was 
used as humidity sensor either without further processing or by introducing it into a 
polymer matrix as wellas by altering it using chemical functional groups by reduction 
to obtain reduced GO (rGO) [6, 7]. Reduction of GO can be carried out by several tech- 
niques, e. g., thermal, chemical or photonic reductions, by which some of the oxygen 
containing functional groups are eventually removed partially [8]. The importance of 
reduction is not only that it converts GO from a semi-insulating to a conductive mate- 
rial but also that it is controllable and tunable. This has a big influence on the GO and 
rGO electrical and sensing properties where they will be functions of type and degree 
of the reductions. 

GO as a promising material for humidity sensors as well as other types of sensors 
still needs to be characterized better to explore its tunability in different aspects. The 
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ability to lose their oxygen groups by reduction and to show different surface, con- 
ductivity and sensitivity properties is of high interest. In this chapter, we investigate 
the effect of the reduction temperature on the humidity sensitivity of GO films. The 
measurements were carried out by taking the impedance of rGO films. 


2 Experimental procedures 


GO-coated thin films were prepared by spin coating. The reduction process was carried 
out by thermal annealing at different temperatures to produce rGO. Electrical contacts 
were made by silver paste. The humidity measurements were taken under certain rela- 
tive humidity (RH) conditions by means of saturated slat solution at RH values of 0 %, 
11.3 %, 32.5 %, 52.3%, 75.2%, 84.3 % and 97.3 % and at room temperature (23°C). Elec- 
trical impedance of the films was measured by the Agilent impedance analyzer 4294A. 


3 Results and discussion 


In Fig. 1, the impedance relation to the humidity of rGO films reduced at different tem- 
peratures is shown. The change in impedance by thermal heating is one order of mag- 
nitude by increasing heat from 160°C to 250°C, where GO itself is semi-insulating, hav- 
ing a resistance of hundreds of MOhms. The sensitivity behavior shows dependency 
on the reduction temperature. At low humidity, the impedance increases for all rGO 
films. At high humidity, the impedance tends to decrease and shows a negative slope. 
However, for the films reduced at higher temperature, the impedance shows a positive 
slope. The relation is quite linear at lower humidity values. 

In the Nyquist plot, Fig. 2 (a) and (b), a semicircle fits the parallel resistance and 
capacitance, which represent the intrinsic impedance of the film. The behavior at low 


= 170°C 
4x10” š . . e 180°C 

, . a 200°C 
1x10 + ~ 250°C 
7x10° Ser 
4x10° G 


Zz (2) 


5x10° 
3x10° 
2x10° 


0 20 40 60 80 100 
Rel. Humidity (%) with RH at 100 Hz. 


Fig. 1: Relation of impedance of rGO films 
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humidity range (0-52%) shows an increase of the semicircle due to absorbed water 
molecules. This causes to charge transfer between water molecules and rGO, which 
reduces the charge carriers of the film [7]. This is applied for all the films. At high hu- 
midity (75-97 %), in the Nyquist plot of Fig. 2 (b), the semicircle with small tail shows 
a constant phase element (CPE) behavior. The equivalent circuit consists of two CPEs 
and a resistance that represents the interaction of water and contact and intrinsic re- 
sistance. 
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For rGO reduced at higher temperatures, by increase in the number of water molecules, 
the charge transfer is dominant and causes a further increase of the impedance. For 
low-temperature reduced films, the formation of water thin films has the role of de- 
creasing the impedance due to proton hopping [9]. In addition, by increasing the 
reduction temperature, rGO films become more hydrophobic, which influences the 
formation of a continuous water film. Thus the impedance only further increases at 
higher values of humidity. 
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4 Summary and conclusion 


GO coated and reduced by temperature leads to tuning not only the conductivity but 
also the film sensitivity to humidity. The impedance can be reduced from hundreds of 
MOhms to tens of kOhms by heating the films to 250°C. The reduction by temperature 
results in different humidity sensitivity behaviors. Tuning the reduction temperature 
provides useful selectivity of the humidity sensitivity. Where at lower temperature is 
more sensitive at higher humidities, at high reduction temperatures, the sensitivity 
to humidity becomes less due to the increase of the hydrophobicity of the rGO films. 
However, it improves the linearity. 

The realization of a humidity sensor or sensor matrix that covers the full range 
of humidity is also foreseen, e. g., by combining two films with two different anneal- 
ing temperatures and thus different transport mechanisms at low and high temper- 
ature. The electrical impedance was obtained by impedance spectroscopy, which re- 
veals the conduction mechanism. However, further tests and study of the equivalent 
circuit models will be a next step, which will provide an in-depth understanding of 
the different behaviors in humidity sensitivity of rGO. 
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Impedance study of Ni-Cr alloy in contact with 
biologic and organic solution using constant 
phase element 


Abstract: The electrochemical characterization of the Ni-Cr dental implant alloy was 
studied by electrochemical impedance spectroscopy (EIS) in contact with artificial 
saliva, milk and vinegar. For all solutions the data impedance has been recorded main- 
taining the sample at 37°C for 1 hour of immersion. In this chapter we are focused on 
the modeling of the response function of this biomaterial using an equivalent electric 
circuit (EEC). The experimental results show that the Ni-Cr alloy exhibits a quite good 
resistance to corrosion, with passive properties referring to the existence of a thick 
passive film. An EEC was proposed to describe the electrochemical behavior of this 
biomaterial in the frequency range of 0—40 kHz. The Bode diagram indicates the pres- 
ence of at least one time constant during the kinetic process. The nature of the film 
and the morphology of the interface biomaterial/film/electrolyte are presented in this 
chapter. The proposed circuit involved constant phase elements in order to explain 
the apparent heterogeneity on the interface biomaterial/film/electrolyte. Simulation 
based on this model showed a good agreement with the experimental data, with a best 
fit corresponding to a minimum standard deviation (the chi-square value x? is 10“). 


Keywords: Ni-Cr, biomaterial, electrochemical impedance spectroscopy (EIS), elec- 
trical equivalent circuit (EEC), constant phase element (CPE), interface biomate- 
rial/film/electrolyte 


1 Introduction 


Ni-Cr-based alloys are part of materials used on orthodontics as dental prostheses. 
Their study has become easier owing to their lower cost, by their biocompatibility and 
by the existence of a variety of software able to fit the data impedance. The Ni-Cr al- 
loys are commonly used for dental materials [1] and widely used in dental skeletal 
structures, implant fixtures and prosthodontic restoration [2]. Beside their low cost 
they exhibit acceptable mechanical and tribological properties and matching thermal 
expansion coefficient with the ceramics of metal-ceramic restorations [3, 4]. There 
also known as passive system with a spontaneous formation of thin passive oxide film 
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on the interface metal/electrolyte [5]. The reactivity of the oxide superficial layer for- 
mation, which is based especially on Cr,03, will be investigated by electrochemical 
impedance spectroscopy (EIS). Moreover, we study the reactivity of the superficial ox- 
ide layer while varying the solution (pH). We used three solutions, namely, artificial 
saliva (pH = 6.8), milk (pH = 6.5 to 6.7) and vinegar (pH = 3), to evidence the reac- 
tivity of the interface. In this work, we compare the Ni-based alloys to the Co-based 
alloys, with their formed protective film depending highly on the acidity of the solu- 
tion used [6]. We study by EIS the response of the biomaterials and the reactivity of its 
passive film, formed with equivalent electric circuit (EEC) involving CPE. 


2 Electrochemical characterization 


2.1 Impedance spectroscopy measurements 


The electrochemical impedance spectroscopy, or also AC impedance, is a non-destruc- 
tive technique and a useful tool for studying corrosion [7]. This experimental method 
makes it possible to distinguish the dielectric and electric properties of individual con- 
tributions of components [5]. 

In this method, the sample undergoes a small electronic excitation so its response 
is pseudo-linear. In a linear (or pseudo-linear) system, the current response (1) and the 
potential (2) are sinusoidal in the same frequency and shifted in phase. 

The sinusoidal perturbation of potential E(t) induces a sinusoidal current I(t); 
their representations in the polar coordinates is as follows: 


E(t) = Ep sin(wt), (1) 
I(t) = Ip sin(wt + p), (2) 


where w = 27f and E, and Ip are the amplitudes of the potential and the current, 
respectively. 

The impedance of the system (3) can be written following the Ohmic concept as 
follows: 
E(t) — Eosin(wt) sin(wt) 


=Z, (3) 


A: It) Ipsin(wt+p) sin(wt +) 


With Euler’s relationship or polar representation of E(t) and I(t), they can be writ- 
ten as follows: 


E(t) = Ege”, (4) 
I(t) = Ie". (5) 
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The impedance complex (6) of the biomaterials can be rewritten as follows: 


Z(t) = - = Ze” = Z,(cosp + sing). (6) 
The modulus of the complex impedance (7) and the phase angle shift (8) are de- 
duced from the following formulas: 


IZ(w)| = VZ + Zin +22, (7) 


tan(p) = om, (8) 


re 


where Z,. and Zim are the real and imaginary part of the complex impedance, respec- 
tively. 

The impedance data are profoundly interpreted by much recently developed soft- 
ware. We have used ZSimpwin version 3.10 to fit the response function of the Ni-Cr 
biomaterial. The EEC consisted of a combination of simple electronic elements (resis- 
tance, capacitor and inductance). The numerical simulation must be done with more 
accuracy to solve difficulties related to the fact that more than one EEC can be able to 
fit the data. 

Studying the response function Z(w) we return to interpret two famous diagrams 
or representations, namely, the diagrams of Bode and Nyquist. There is some liter- 
ature which talks about the basics of EIS [8-13], and the fields of its applications 
cover many domains, from biology and medicine to industry. The electrochemical 
impedance analysis has been performed on the commercial Ni-Cr alloy in the fre- 
quency range of 0-40 kHz, after it had been immersed for 1 hour in milk, artificial 
saliva or vinegar. 


2.2 Experimental results 


The Bode diagrams (Fig. 1 and Fig. 2) of the Ni-Cr dental alloys recorded in milk, artifi- 
cial saliva and vinegar show that at least one time constant phase element (CPE) can 
be distinguished. For artificial saliva and vinegar, the shape of the curves are similar. 
At high frequencies the absolute impedance curve is almost independent of the fre- 
quency, with a phase angle of 0°, and the impedance of the system is reduced to the 
electrolyte resistance. 

For medium frequencies, at 1-100 Hz, a linear frequency dependence of absolute 
complex impedance with a maximum angle reaching 70°, 65° and 55° for milk, artificial 
saliva and vinegar, respectively. This range of frequency is marked by the stable film 
formation. 

At low frequencies, the absolute impedance exhibits its high values for milk and 
artificial saliva, and these values are ten times higher than for vinegar. This differ- 
ence in the values of the impedance gives an idea of the film formation for the three 


56 — N.Leguedani etal. 


a Milk 


+ Saliva 


Vinegar 


Phase angle (degrees) 
& 
x 


10 
Fig. 1: Phase angle plot for Ni-Cr for milk, 
£ 0 1 10 100 1.000 10.000 100.000 saliva and vinegar. The fit based on the 
Frequency(Hz) model is represented by solid lines. 
100.000 
+ Milk 
*, 4 Saliva 
10.000 x t Vinegar 
* 


ZI(Q.cm?) 


Fig. 2: Measured (discrete points) and 
fitted (solid lines) values of the Z modulus 
o o 1 10 100 1.000 10.000 100.000 for Ni-Cr dental alloys for milk, saliva and 
Frequency(Hz) vinegar for 1 hour immersion at 37°C. 


solutions under study. This can be explained by the fact that the film is of a porous 
character in origin. x 

The Nyquist plots of Ni-Cr dental alloys are shown in Fig. 3 for the three solutions 
(milk, artificial saliva and vinegar) with the potential corrosion E.orr and at 37°C. It 
can be seen from the impedance data that the Ni-Cr implant biomaterial/solution in- 
terface exhibits capacitive behavior over a relatively wide frequency region, which is 
typical for a passive alloy system. The parameters for the solution bulk resistance R, 
(experimentally determined by the intersection at high frequency of the Nyquist plot 
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Fig. 3: Nyquist diagram of Ni-Cr dental alloys in the three solutions of milk, artificial saliva and vine- 
gar for 1hour immersion at 37°C (discrete points). 


with the real axis) are of about the same order of magnitude in the three solutions; 
12.63 Q.cm? for milk, 19.88 Q.cm? for artificial saliva and 30.9 Q.cm? for vinegar. 

The resistance of polarization, R, (experimentally determined by the diameter of 
the semi-circles on the arc in the Nyquist plot), related to the corrosion resistance, is 
found to have high values that increase significantly with milk, artificial saliva and 
vinegar, respectively. 

This high value of R, implies higher corrosion resistance of implant alloy for milk, 
reaching a value of 25000 Q.cm?, while this value is about 1600 Q.cm? for artificial 
saliva and 1000 Q.cm? for vinegar. 

The low R, values in vinegar show that the Ni-Cr dental alloys have a weak cor- 
rosion resistance when compared to milk and artificial saliva; then a high rate of re- 
leased metallic ions into the oral environment can take place. The Cole-Cole plot of 
the complex admittance (Fig. 4) shows the same geometry and the arc dependence 
of the semi-circle is affected by the nature of the solution used. The diameter of the 
semi-circle decreases with vinegar, artificial saliva and milk, respectively. 


2.3 Modeling the impedance data 


The proposed equivalent circuit based on CPE is described in Fig. 5, where R, is the 
electrolyte resistance, in series with Q,,; and in parallel with R,,;, which is in se- 
ries with Qin, which is in parallel with R;„. Here, Qout denotes the pseudo-capacitance 
(CPE,,;) and Rout the resistance of the Ni-Cr protection film system, whereas Qin and 
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Fig. 4: Admittance diagram of Ni-Cr dental alloys in the three solutions of milk, artificial saliva and 
vinegar for 1 hour immersion at 37° (discrete points). 


Fig. 5: EEC used in the generation of simulated data. 


R,, denote the double layer pseudo-capacitance (CPE;,,) and the transfer resistance of 
the corrosion of the Ni-Cr/solution interface. 

CPE (9) was mathematically used to solve the response of the biomaterial with a 
non-ideal capacitor interpreting the non-homogeneous behavior of the surface under 
investigation. The main reasons for adding CPE in a matrix of EEC were the distributed 
surface reactivity, surface inhomogeneity, roughness and electrode porosity [6]. This 
lack of homogeneity is modeled with a Q element, used to represent the CPE, i. e., 


Zo(w) = (9) 


1 

Yo Gw)” 
where Yọ is the admittance of an ideal capacitance and n is an empirical constant, 
ranging from O to 1. The exponent n is related to a non-uniform current distribution 
due to the surface roughness [7, 8]. It is noteworthy that when n = 1, the CPE behaves 
as a pure capacitor, while when n = 0, the CPE behaves as a pure resistor and when 
n = 0.5, the CPE is the equivalent of the so-called Warburg element. 

The impedance ofZCPE is defined by ZCPE u: = 1/(jw)™ Q; and ZCPE;, =1/(jw)"?Q>, 
where j? = -1, w is the angular frequency and the exponents n, and n, of the CPEyut 
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and CPE;, are related to a non-uniform current distribution due to the surface rough- 
ness [12]. The chi-square values 0) for all data impedance are of magnitude 10° in- 
dicating an excellent agreement between the experimental data and the model using 
CPE in the fitting program. 


3 Simulation criterion and discussion 


Taking into account that several circuits can fit the same data impedance, we aim to 
determine which of the circuits is more suitable to best describe the corrosion phe- 
nomena. We used two criteria to decide which circuits numerically fit the impedance 
data in search of an optimal configuration. First, only the simple elements introduced 
into the EEC must have a physical meaning and second, the chi-square value (y?) must 
be less than 10°*. Monitoring numerically the evolution of parameters, one can deter- 
mine which of the proposed models can be kept. Systematically, the standard devia- 
tion in the present work is found to be 1.732 x 10-™ for milk, 3.895 x 10% for artificial 
saliva and 1.693 x 10° for vinegar. 

From this model it was possible to calculate the polarization resistance, the resis- 
tance R, of the solution, the relaxation frequency of the system and the CPE parame- 
ters (the resistance and the capacity of the oxide layer). Thus a simple schematization 
of the interface (Fig. 6) can be given. 


Electrolyte 


Fig. 6: Modeling of the interface electrode/film/elec- 
Ni-Cr alloy trolyte of Ni-Cr dental alloys. 


Table 1 lists the parameters of the EEC used to simulate the impedance data of the Ni- 
Cr implant alloy in the presence of milk, artificial saliva and vinegar, maintained at 
37°C and for one hour immersion time. 

It was found that the values of the resistance of the electrolyte, R,, for the Ni-Cr 
implant alloy is increasing with milk (12.63 Q.cm2), artificial saliva (20.06 Q.cm?) and 
vinegar (31.06 Q.cm?). The resistance of polarization (Fig. 7) is calculated in this model 
by the sum of Rin and Rout. The value of R, decreased; it is 33.588 kQ.cm? for milk, 
28.579 kQ.cm? for saliva and 991.2Q.cm? for vinegar. 
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Table 1: Model parameters used to simulate the impedance spectra of the Ni-Cr implant alloy in the 
presence of milk, artificial saliva and vinegar. 


Rs Qout nı Rout Qin n2 Rin 
Q.cm?  pF.cm? Q.cm? pF.cm? Q.cm? 
Vinegar 31.06 313.6 0.82 715.5 3734 0.72 275.7 
Milk 12.63 171.9 0.82 2848 35.5 0.89 3.07x10% 
Saliva 20.06 163.5 0.8 699 380.4 0.44 2.79 x 104 
35000 
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Fig. 7: Resistance of polarization, R,, of Ni-Cr alloys in the three solutions of milk, artificial saliva 
and vinegar. 


The theoretical model predictions are in good agreement with the experimental 
data. Thus, the Qut capacitance values in all solutions were within a range of 172- 
314 uF.cm?, with n within a range of 0.8029-0.8213, which indicates a non-ideal ca- 
pacitance interface (so an oxide passive layer exists). 

Moreover, the Q;, capacitance values are highly dependent on milk, artificial 
saliva and especially of vinegar; the values of n have the same order in milk (0.8924) 
and in vinegar (0.7151), but not in saliva, which shows a value of 0.4361. 


4 Conclusion 


The corrosion interface of Ni-Cr alloys in biological (milk and artificial saliva) and or- 
ganic solutions (vinegar) was studied by EIS based on a model using CPE. Results 
show that this biomaterial exhibits a good corrosion resistance in milk and artificial 
saliva and a weak corrosion resistance in vinegar. All the n values in the three solutions 
were found to be within a range of 0.7151 to 0.8924, indicating that the corrosion in- 
terface deviated from a pure ideal capacitor. The phase angle maximum (Pmax) in this 
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biomaterial is measured at approximately -53 to -70 degrees. The EEC predictions are 
in good agreement with the experimental data. 
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Part Il: Measurement systems 


Thorsten Baumhöfer, Susanne Rothgang, and Dirk Uwe Sauer 
Influence of the settling behavior on the 
measurement uncertainty of EIS and how to 
decrease the measurement duration 


Abstract: To obtain a spectrum, different frequencies are applied stepwise. In an ideal 
system the measurement uncertainty is only influenced by the noise of the voltage and 
current measurements. But since a switched sine includes transient signal parts, the 
system’s settling behavior leads to systematic errors. Usually the first periods after the 
start of the excitation are thus not included in the measurement. A method based on 
a model of the measurement system will be shown, allowing for compensation of the 
error. Thus the measurement duration can be reduced significantly. 


Keywords: Uncertainty, settling behavior, measurement model 


1 Introduction 


The electrochemical impedance spectroscopy is widely used in laboratories as a tool 
for the parameterization of electrochemical energy storage models based on equiv- 
alent circuits [1]. It is also used as a tool to investigate changing processes in aging 
batteries [4]. A good knowledge of the measurement uncertainty is needed in both 
cases to be able to evaluate the determined parameters. 

Another growing field of application is quality assurance during the manufactur- 
ing of battery packs. In a series production line the reduction of cycle time has top 
priority. Therefore, the goal is to decrease the measurement duration of an impedance 
spectrum, without increasing the measurement uncertainty considerably. 

The most common way to measure the electrochemical impedance is to apply a 
sinusoidal excitation signal. The excitation and response signals are acquired using 
a constant sample rate. Both signals are then transferred to the frequency domain by 
using the discrete Fourier transformation and divided to obtain the impedance. To 
acquire a spectrum, different frequencies are applied stepwise, one after another. 

To reduce the measurement duration, multiple frequencies can be applied simul- 
taneously [2, 3]. Since the information density of an electrochemical system is not 
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scaled linearly with frequency, only a finite number of distinct frequencies has to be 
included in the signal. The phase angles can further be chosen in a way to minimize 
the crest factor [5]. Overall suitable signal-to-noise ratios can be reached. 


2 Measurement uncertainty 


The main influence on the uncertainty of measured impedance values, which cannot 
be reduced by calibration, is the noise of the voltage and current measurement. 

With the Fourier transformation the noise component at the examined frequency 
gets isolated and results in a Gaussian distribution around the true impedance 
value [6]. This behavior was verified by 22,000 measurements on a Sanyo/Panasonic 
UR18650E Li-ion cell, each with one period at 220 Hz. The resulting histogram (Fig. 1) 
shows the analytically proven Gaussian distribution. 


Fig. 1: Gaussian distribution of 
impedance measurement values 
due to noise in the excitation and 
Im{Z} (mQ) Re{Z} (mQ) response signals. 


The standard deviation of the distribution depends on the noise magnitude at the con- 
sidered frequency, the magnitude of the excitation and response signals and the num- 
ber of samples used for the calculation. 


3 Settling influence 


Electrochemical energy storage devices exhibit processes with different time con- 
stants. This results in a transient response, when the excitation signal is being 
switched on. It takes a specific time until the amplitude and phase shift can be consid- 
ered stable. To apply a Fourier transformation to this unstable part can be considered 
wrong, since the assumption of a periodic signal is not valid. But still it can be com- 
puted without numerical difficulties and provide information about the system. 
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To quantify the influence of the settling 20 periods at 220 Hz where applied to the 
previously introduced Li-ion cell. This procedure was repeated 100 times with breaks 
without excitation in between. Every period was individually evaluated, resulting in 
the impedance points seen in Fig. 2. For the first cycle the error due to the settling 
behavior is quite significant compared to the uncertainty due to noise. The error de- 
creases exponentially with increased cycle number, but in this case has to be consid- 
ered as a systematic error up to cycle number five. 


Overview 
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The error is not constant but depends on the time constants and resistive magnitude 
of the analyzed system and changes in magnitude and direction with the frequency 
being applied. 

The usual way to overcome this effect is to apply the first period but start the ac- 
quisition afterwards, during the second cycle. This is documented by the different in- 
strument manufactures in different detail (e. g., as “pre-wave” [7]). Beginning from the 
second cycle the following cycles are usually averaged to reduce the noise level. The 
effect of the settling is also being reduced by this averaging. 
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4 Reduction of measurement duration 


If the measurement duration should be reduced, long averaging cannot be used any- 
more. It would further be favorable to not have to apply at least two periods of the 
signal. Especially at lower frequencies, in the mHz range, this is considerably influ- 
encing the measurement duration. 

In Fig. 3 the resulting spectrum of the first excitation period and the true system 
response (from the 20th period) are shown. For these measurements a multisine signal 
with 39 frequencies, logarithmically spaced from 8 kHz to 100 mHz, was applied to 
the previously introduced cell. The length of the signal is 10 seconds (one period of 
the lowest frequency). It can be seen that the influence in the lower frequencies is 
greater than in the higher frequencies. This is due to deviating numbers of periods of 
the different frequency components and therefore a variation in the averaging. 


Im{Z} (mQ) 


30 40 50 60 70 80 90 100 excitation. The settling error can clearly 
Re{Z} (mQ) be identified within the first period. 


In the measurement model (see Fig. 4) the signal source, amplifier and input filter are 
assumed to be ideal. The multisine signal is periodic, so the gating in the current mea- 
surement path has no effect. But the gating of the voltage signal is the source of the 
settling error and thus has to be considered in the model. 


Specimen Input-Filter Gate u 
H(w) Hw) H,(w) 
Input-Filter 
Hw) 


Signal-Source 


Current- 
Amplifier 


Current-Shunt 


Fig. 4: Model of the measurement system. 
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Overall, equation (1) can be derived from the model, where H(w) denotes the true 
impedance of the system, H,w is the transfer function of an ideal rectangular window 
and V,I can be observed. A response impedance H* (w) can be calculated, including 
the settling error. The equation cannot be de-convoluted to be solved with respect to H, 
since the problem is ill defined and would lead to massive noise amplification. We have 


V(w) = (I(w) - H(w)) « H,(w), (1) 
AC) 
H” (w) = Te)’ (2) 


To reconstruct the true impedance, an iterative method has to be developed only 

computing the measurement model in forward direction. 

1. Set the system impedance input H,„-o to the measured response impedance H* 
interpolating in between the measured frequencies. 

2. Calculate the new response impedance H, with the measurement model. 

3. Change the input by the difference of the responses H,,,, = H, + (H,, -H,_,) at the 
frequency points excited. Apply an interpolation as above. 

4. Repeat from step 2 until the iteration converges. 


An interpolation can be used without introducing errors, since electrochemical en- 
ergy storage devices exhibit continuous impedance functions not including disconti- 
nuities. 

The first three iterations on the previously shown measurement results are shown 
in Fig. 5. It can be seen that the curve of the input impedance starts to smoothen and 
the response impedance converges towards the measured curve. It is therefore pos- 
sible to reconstruct the true impedance of the system only taking the measurement 
values from the first period into account. 
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5 Conclusion 


The settling behavior is a significant error source for impedance determination within 
the first to sixth excitation period. A method was given to reconstruct an impedance 
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spectrum utilizing a measurement model. It is sufficient to excite the system with a 
multisine excitation with the length of one period of the lowest frequency. The mea- 
surement duration for impedance spectra can thus be significantly reduced. 
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Achraf Lamlih, Vincent Kerzérho, Serge Bernard, Fabien Soulier, 
Mariane Comte, Michel Renovell, Tristan Rouyer, and 
Sylvain Bonhommeau 


Mixed-level simulation tool for design 
optimization of electrical impedance 
Spectroscopy systems 


Abstract: The design of electrical impedance spectroscopy systems is driven by target 
applications. Moreover, system functions are inter-dependent. Thus, optimizing the 
design of each block without input from the others is impossible. Furthermore, simu- 
lating the whole system at transistor level is time consuming. This chapter presents 
a mixed-level simulation tool for the efficient development of bioimpedance spec- 
troscopy measurement systems. We have developed Verilog-AMS models combined 
with transistor-level descriptions of several analog blocks. This allows us to run highly 
efficient system-level simulations for a whole circuit. We apply this tool to the design 
of an integrated circuit suitable for in vivo monitoring of bluefin tuna physiological pa- 
rameters. In fact, the variation of biological processes offers significant insights into 
a host of key parameters such as growth, survival and reproduction. The obtained re- 
sults show that the approach described in this chapter can reduce simulation time by 
80 % while consistently reaching high simulation performances. 


Keywords: bioimpedance spectroscopy, IC design, bluefin tuna, physiological param- 
eter monitoring 


1 Introduction 


Electrical Impedance Spectroscopy (EIS) can be used to characterize matters governed 
by electrochemical processes. From material science to biosensing [1-3], the EIS range 
of applications is wide. In some of these applications, the system design is ruled by 
strong constraints. The simple choice between various potential system architectures 
is challenging. Moreover, system functions are inter-dependent. As a result, optimiz- 
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Fig. 1: The main blocks of the integrated circuit architecture. 


ing the design of each block without input from the others is impossible. To address 
this issue, an efficient design strategy is needed. 

The overall architecture consists of the main blocks presented in Fig. 1. The three 
principal system functions are the frequency components generation, the stimulation 
waveform construction and the system-under-test response analysis. The hardware 
implementation of each of these blocks differs depending on the application. 

In general, three approaches are used for the design of integrated systems, i.e., 
the bottom-up, the top-down and the mixed-level approaches. The bottom-up ap- 
proach consists in implementing each of the system blocks at transistor level, pro- 
ceeding to the block verification, and then assembling all the blocks to build the sys- 
tem. This approach helps to achieve high performances at block level but without any 
system-level consideration. Thus, it is impossible to ensure optimal system-level per- 
formances. Moreover, a global simulation at transistor level is generally impractical 
because it is time consuming. Conversely, the top-down approach addresses archi- 
tectural level first by describing each block by its functional model. This approach 
can optimize system-level specifications but without considering each block con- 
straint. The high-level block requirements retrieved from the system-level simulation 
might be unsuitable at low level. The mixed-level approach combines high-level and 
low-level models for each block. The idea is to co-simulate some of the blocks at the 
transistor level with the others at high level. With this approach, we can address the 
inter-dependability of the system blocks to propose various optimized architectures, 
while customizing each block for better performances. Furthermore, by combining 
fast-simulated high-level models with accurate low-level models, we can drastically 
reduce the global simulation time. In this context, we develop a mixed-level simula- 
tion tool dedicated to help designers find the most suitable architecture for certain 
given system constraints. 

Since the tool is intended to help design integrated circuit systems, the Cadence 
platform was chosen. We have developed a library of high-level and low-level mod- 
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els of the mainly used blocks for impedance spectroscopy from the literature [4-6]. 
For high-level models we use Verilog-AMS, which is a derivative of the Verilog hard- 
ware description language that enables designers to create analog and mixed-signal 
modules containing high-level behavioral descriptions of the system blocks. For the 
low-level model description, we use transistor-level models simulated by the Cadence 
Virtuoso suite. 

We apply this tool to the design of an integrated circuit suitable for in vivo moni- 
toring of bluefin tuna physiological parameters. 


2 Bluefin tuna physiological parameter monitoring 


Although marine ecosystem services are critical for human well-being, the necessary 
scientific information to conserve and manage them efficiently is currently lacking. 
A way to improve our understanding of marine species is to use implantable sensors 
for monitoring physiological parameters related to key biological processes such as 
feeding and spawning. Indeed, such processes offer significant insights into a host of 
key population parameters such as growth, survival and reproduction. For instance, 
a key physiological parameter is muscle fat content. The variation of this parameter 
in space and time is indeed related to the previously mentioned behaviors as fish feed 
to accumulate reserves, as the fat is used during reproduction. 


2.1 System specifications related to the targeted application 


The choice of the different elements constituting the architecture of the integrated cir- 
cuit is related to the application. In our case, long monitoring periods in the wild are 
considered, typically 12 months. Therefore, the architecture should be characterized 
by low power consumption. As the device has to be implantable, solutions offering a 
low silicon surface area have to be preferred for the fish well-being. 

Previous measurements conducted on bluefin tuna showed that the expected 
impedance values are of the order of 100.0 at 50 KHz [7]. At higher frequencies, this 
value is expected to decrease. Using a handheld measurement device based on the 
AD5933 chip and a four-electrode configuration, bioimpedance spectroscopy ranging 
from 300 Hz to 100 KHz was performed on the bluefin tuna’s back at the second dorsal 
fin level. The obtained results (Fig. 2) confirm the expected impedance range: from 
sub-Ohm levels at higher frequencies to hundreds of Ohms at lower frequencies. 

The non-linear behavior of tissue could be linked to biological processes [8], thus 
allowing the collection of new data with substantial scientific potential. Therefore, the 
designed system must be able to measure non-linearity of the tissue. 
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Fig. 2: Bluefin tuna impedance magnitude. 


For a rigorous characterization of the biological tissue composition, the complex 
impedance needs to be analyzed over a large range of frequencies. Each type of tis- 
sue molecule has a different response to given frequencies. Three main frequency 
regions describing physiological processes have been identified [9]: the a dispersion 
(Hz-kHz), the £ dispersion (kHz-MHz) and the y dispersion (MHz-GHz). Since a and 
ß dispersions are associated with the most relevant biological aspects [4] for our ap- 
plication and for integrated circuit (IC) design considerations, the exploration was 
limited to these two frequency regions: 500 Hz-10 MHz. 


2.2 IC architecture definition 


For the general principles of the measurement, current stimulation and voltage mea- 
surement were chosen, mainly for safety reasons since in that case the current passing 
through the biological tissue can be controlled. Also, a four-electrode configuration is 
used to get rid of the impedance at the interfaces. 

To generate the a and ß dispersions we use a phase-locked loop as a frequency 
multiplier. The low-frequency signal coming from a crystal oscillator is multiplied to 
get the higher frequency. On the negative feedback path, successive divisions are then 
performed to get the frequency points. 

To analyze non-linearity of the tissue, sine waves have been chosen for the stim- 
ulus. In order to measure harmonics of the system (tissue), the generated sine wave 
should have the lowest total harmonic distortion THD possible, so that harmonic inter- 
modulation does not affect the measured result. To cancel harmonics of the sine wave, 
scaled and delayed versions of the signal at the output of the frequency sweep genera- 
tion block are added (Fig. 3). The scale coefficients and the phase shifts are computed 
so that the constructed signal matches the expression of a pure sine wave. 

The result is then fed to a voltage-controlled current source (VCCS) that generates 
the stimulation current. The VCCS is designed in such a way to have high stability cur- 
rent in the targeted frequency range. Synchronous detection has been chosen as the 
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response analysis technique for its simplicity and its ability to single out the wanted 
frequency response even for the smallest AC signals obscured with noise. Since the fre- 
quency components of interest are translated into the baseband by the synchronous 
detection, the constraints on the analog-to-digital converter (ADC) are then relieved, 
and simple low-rate ADCs could be used to feed the memory with bit words repre- 
senting the real and the imaginary parts of the impedance. A finite-state machine is 
used to control the different signals of the architecture, ensuring the sequentiality of 
the frequency sweep. Since the bluefin tuna is expected to migrate in several regions, 
the environmental parameters of the measurement (temperature, pressure, etc.) may 
vary drastically. A reference impedance is added to the overall architecture (Fig. 4) to 
compensate for the system drifts. 


Memory 
Charge 
Pump 


Reference 
Impedance 


Fig. 4: The system-on-chip architecture. 


2.3 Mixed-level simulation of the architecture 


The VCCS is a critical component in bioimpedance spectroscopy architectures. It 
should deliver a high-accuracy current to the tissue under test over the frequency 
range of use. Unfortunately, at higher frequencies, the stray capacitance shunts the 
output impedance of the current driver, which causes its performance to decrease. 
To illustrate the mixed-level simulation approach on the proposed architecture, we 
chose to implement the VCCS at the transistor level to target its performances and the 
other blocks at high level using Verilog-AMS. 
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2.3.1 Voltage controlled current source design 


The current driver is implemented using a symmetrical operational transconductance 
amplifier (OTA) (Fig. 5). It is composed of a differential pair (M1, M2) driving two loads 
(M3, M4) connected as diodes, three current mirrors (M3-M5, M4-M6, M7-M8) and a 
current source M9 [10]. The first stage (M10, M11, M12, M13) is designed to bias the OTA 
with a current I). 


Vdd 


Gnd 


Fig. 5: The current driver architecture. 


The OTA is designed with a 180-nm CMOS process, for a bias current I, of 50 pA, 
a dynamic input range of +0.3 V (Fig.6) around a common mode voltage of 0.6 V, 
a transconductance of 200 uS and a bandwidth of 10 MHz. 
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Fig. 6: The OTA differential DC input-output transfer curve. 


The achieved output impedance of the current driver is 270 KQ at frequencies below 
5 MHz; it drops to 259 kQ at 10 MHz (Fig. 7). The output impedance value could be in- 
creased by cascoding the output stage, for example, at the price of decreasing the 
output swing of the current driver. Although these values of output impedance are 
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Fig. 7: The current driver output impedance vs frequency, with a load of 1kQ/3 pF. 
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Fig. 8: The current driver transconductance vs frequency, with a load of 1 kQ/3 pF. 


not high, they are suitable for our application since the expected impedance of the 
bluefin tuna does not exceed 100 Ohms (Fig. 2). The achieved transconductance value 
is 190 uS and remains stable over the desired bandwidth (Fig. 8). 

Monte Carlo simulation was performed on the transconductance of the OTA for 
1000 runs with both process and mismatch; the results are shown in Fig.9. The 
transconductance has a mean value of 190.64 uS and a standard deviation of 3.68 uS. 
The value of the standard deviation could be decreased by increasing the dimensions 
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Fig. 9: Monte Carlo simulation of the transconductance of the current source. 


of the current mirrors of the output stage, at the price of decreasing the bandwidth 
due to larger parasitic capacitances. Since 97% of the runs are within 3.8 % of the 
deviation from the mean value, the design is considered as robust under process 
variation and device mismatch, and the transistor dimensions are kept unchanged. 


2.3.2 Simulation performances 


To illustrate the efficiency of the mixed-level approach, two transient simulations 
of the overall architecture described below were launched. The first simulation con- 
tained the current source designed at the transistor level and the other blocks at high 
level. The second simulation contained all the blocks of the architecture at high level; 
a realistic Verilog-AMS model of the current source was used. 

The Verilog-AMS model takes into account the values of the transconductance and 
the output resistance extracted from the simulation results presented in the previous 
part. The model could also be used to evaluate at high level the impact of process 
variation and device mismatch by using the results of the Monte Carlo simulation. 
In fact, the impedance error introduced by process and mismatch variation is +1.9 % 
within one standard deviation from the mean value. 

Simulations were performed on a machine with 25 cores (2 GHz) and 263 Gb of 
memory. The entire high-level simulation with the Verilog-AMS model of the current 
source took 15 minutes. On the other hand, the simulation with the current source at 
transistor level and the other blocks at high level took 75 minutes. The obtained results 
show that the mixed-level approach described in this chapter can reduce the simula- 
tion time by 80 % while consistently reaching high performances for each block. 
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3 Conclusion 


This chapter presents a mixed-level simulation tool to support the design of bio- 
impedance spectroscopy ICs. The mixed-level simulation addresses the inter-depend- 
ability of the different blocks of the architecture, while enabling the designers to test 
different architecture solutions. Furthermore, the tool permits to go further in the 
design customization, allowing designers to enhance system performances with the 
significant advantage of having reduced simulation times. Moreover, this approach 
could also be useful to estimate the influence of process variation and device mis- 
match on the overall architecture. The influence of the measurement environment 
could also be estimated. In fact, Verilog-AMS blocks modeling sources of inaccuracies 
(such as the effect of the temperature and flicker noise) could simply be added to the 
model, permitting the evaluation oftheimpact and ifneeded refining the architecture 
for robustness. 
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René Schmidt, Stephan Blokzyl, and Wolfram Hardt 
A highly scalable FPGA implementation for 
cross-correlation with up-sampling support 


Abstract: Cross-correlation is a fundamental algorithm which is used in many fields 
of digital signal processing. One important application is sound source localization. 
For this purpose, actual research shows a significant improvement of localization ac- 
curacy when using up-sampling-based algorithms. These state-of-the-art approaches 
are challenging for embedded realizations as up-sampling in combination with high 
sample rates and measurement window lengths leads to huge amounts of data to be 
processed and stored. 

This chapter presents hardware acceleration techniques for cross-correlation with 
up-sampling support. The work introduces a novel processing architecture which ad- 
dresses large-scale memory consumption compared to conventional solutions. The 
novel design introduces a highly scalable field programmable gate array realization 
with linear memory complexity and an interpolation factor-independent runtime. 


Keywords: mathematics algorithms, digital signal processing, splines, field pro- 
grammable gate arrays 


1 Introduction 


Cross-correlation algorithms are widely used in various fields of digital signal process- 
ing. Typical applications comprise image processing [1], audio signal processing [2], 
impedance spectroscopy [3] and sound source localization [2]. Especially for sound 
source localization, approaches like time difference of arrival (TDOA) are used, a tech- 
nique to calculate the position of a sound source by computing the delay between 
incoming sound signals of spatially divided microphones [4]. The cross-correlation is 
used to calculate the audio signal lag (signal time difference). There are two common 
models for cross-correlation. The general definition is 


u(t) @ v(t) = | u“ (T) v(t + ddr. (1) 
Since digital signals are discrete and real-valued, the definition (1) can be written 


u(t) @ v(t) = + u(t)v(t +0). (2) 
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Alternatively, the cross-correlation can be calculated by using convolution, given by 
u(t) * v(t) = | u(t)v(t + ddr. (3) 
An equivalent description is 
u(t) * v(t) = F \(F(u(t)) » F(v(t))). (4) 


The calculation (4) has recently been used in micro-controller-based realizations, 
since the fast Fourier transform (FFT) F can be implemented efficiently, in contrast 
to the traditional estimation (2), which has a runtime of O(NT ax) = O(n?) [5]. FFT 
calculation is also appropriate for field programmable gate array (FPGA)-based solu- 
tions [6]. Hence, predesigned FFT IP cores are resource-intensive and restricted in the 
number of samples (compare Section 3). 

Two major properties specify TDOA localization accuracy. One is the base distance 
between the microphones, which is proportional to the measurement precision. Small 
base distances cause high measurement uncertainty [7]. However, the base distance 
cannot always be adapted to a length which enables the required measurement accu- 
racy. The second feature is the sampling rate. The measurement error reduces with a 
higher sampling rate, but sample rate incrementation leads to a higher data volume to 
be processed. With a measurement interval of 20 ms and a default microphone sam- 
ple rate of 44.1kHz, 882 values have to be correlated each measurement cycle. This 
amount of data is not a problem, but with increasing sample rate the computational 
effort raises significantly, which is challenging for embedded devices. 

The sample rate can be increased by using faster ADCs or by generating addi- 
tional artificial data with up-sampling. Kan has presented an up-sampling method 
to improve localization accuracy [2]. The challenge can be illustrated with the follow- 
ing example. An up-sampling rate of 64 together with a sample rate of 44.1 kHz leads 
to an increased sample rate of approximately 2.8 MHz and 56,448 data samples to be 
computed. The restriction of, e. g., the Xilinx reference FFT IP core to 65,536 data sam- 
ples [8] allows a maximum sample rate of approximately 3.3 MHz. Higher sample rates 
cannot be realized without measurement window reduction. 

This chapter is structured as follows. Section 2 introduces related work and ap- 
plications and validates state-of-the-art cross-correlation methodologies with respect 
to embedded realizations. Section 3 gives a review of an FPGA-based reference design 
that is based on common cross-correlation as presented in standard literature. Sec- 
tion 4 presents a novel architecture for a highly scalable FPGA implementation of up- 
sampling supporting cross-correlation. Finally, Sections 5 and 6 conclude the chapter 
with the results and benefits of the proposed solution. 
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2 Related work 


There is a considerable amount of literature discussing versatile applications of cross- 
correlation. The proposed approaches apply software-based cross-correlation speci- 
fied and optimized for different measurement applications. Venable presented a work 
in 2006 using cross-correlation for the estimation of peptide molecular weights from 
tandem mass spectra [9]. With the tool SEQUEST, a cross-correlation-based template 
matching approach is used for spectrogram comparison [10]. 

In the field of image exploitation, the fast normalized cross-correlation has been 
established as novel template matching (TM) method [1, 11]. The TM-based cross- 
correlation has been used for medical applications like motion recognition in blood 
flow imaging [12] and frontotemporal dementia [13] and for quality assurance in in- 
dustrial production processes like defect detection of electronic devices [14]. Other 
applications use cross-correlation to calculate time differences between two related 
signals, like sound source localization methods [4], as mentioned above. 

The listed algorithms are realized in software and cause disadvantageous high 
computational costs [5], which is challenging in embedded applications. Therefore, 
FPGA technology is applied to accelerate computation with the help of parallel pro- 
cessing strategies. FPGA-based cross-correlation architectures, as proposed in the lit- 
erature, can be classified into two main categories: solutions using discrete cross- 
correlation and methods applying FFT. Approaches from the first category solve equa- 
tion (2) and implement a parallel or serial architecture [15] as depicted in Fig. 1. These 
architectures are used in the field of orthogonal frequency division multiplexing [16] 
and other high-performance applications [17]. The serial structure provides a high pro- 
cessing speed, since each sample X needs to be a processed one. However, resource 
consumption rises linearly with the amount of samples per measurement window. As 
each measurement sample requires a single multiplier for cross-correlation, e. g., 20 
samples consume 20 and 20,000 samples consume 20,000 multipliers. The relation 
between the measurement sample and multipliers is linear. With 220 DSP cores for 
multiplication, as provided by state-of-the-art FPGAs,’ the proposed implementation 
is not feasible for measurement scenarios with high sample rates or large-scale mea- 
surement windows. 

In contrast to the serial approach, the parallel architecture saves resources, since 
the calculation needs to be executed for each possible lag [15]. However, the paral- 
lel architecture consumes a big amount of multipliers to realize cross-correlation for 
high sample rates. As the count of multipliers is limited in FPGAs, both serial and par- 
allel solutions are not feasible for high sample rates and large measurement windows. 
Thus, state-of-the-art FPGA-based solutions [6, 18] use FFT-based cross-correlation re- 
alization with convolution (see equation (4)). Fig. 2 shows the appropriate hardware 


1 E.g., Xilinx Zynq®-7000 All Programmable SoC XC7Z020-CLG484. 
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Fig. 1: Comparison of common FPGA architectures for direct calculation of discrete cross-correlation. 
(a) The stepwise, serial realization of cross-correlation. All X-samples have to be processed for the 
final cross-correlation result. (b) and (c) The highly parallel architecture for cross-correlation which 
provides a result update every clock cycle. 
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Fig. 2: FFT-based cross-correlation architecture with linear spline interpolation-based preprocess- 
ing. 


implementation applying parallel FFT calculation. Various IP cores for efficient, high- 
performance FFT computation are available, easy to use, and optimized for specific 
FPGA devices, which makes them comfortable for modern FPGA designs. 

A literature survey and manifold application scenarios prove the cross-correlation 
is a potential method for a wide range of commercial, industrial and research fields. 
The high computational effort and significant computational resource utilization 
represent the major drawbacks of state-of-the-art methodologies, especially for em- 
bedded realizations. The following analysis discusses standard FPGA-based cross- 
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correlation in detail and derives an approach facing specific resource limitations of 
FPGA technology while providing higher computational performance compared to 
conventional cross-correlation approaches. 


3 Reference design 


The FPGA-based hardware-accelerated design for cross-correlation, depicted in Fig. 2, 
has been implemented on an evaluation board with Xilinx Zynq®-7000 SoC to com- 
pare and assess different realization alternatives and implementation techniques. Two 
parallel FFT modules simultaneously transform two time domain input signals A(t) 
and B(t) from spatially separated microphones into a frequency domain. The design 
substitutes time domain cross-correlation with multiplication in the frequency do- 
main of signal A and the complex conjugated of signal B. 

To support the up-sampling methodology by Kan, which improves the measure- 
ment accuracy (compare Section 1), both signal a(t) and signal b(t) have to be inter- 
polated to A(t) and B(t) prior to FFT. A common interpolation method is, e. g., linear 
spline interpolation (LSI). The inverse fast Fourier transform (IFFT) converts the multi- 
plication result back to the time domain, which represents the cross-correlation result. 
The most prominent peak of the cross-correlation result characterizes the time offset 
between signal A and signal B. 

The two FFTs have been implemented as dual-channel FFT, which provides the 
best resource optimization provided by the vendor Xilinx. Fig. 3 (a) shows the appro- 
priate utilization of Block RAM (BRAM) and DSP resources in relation to the measure- 
ment window length. The number of DSP slices is constant and independent of the 
measurement window size. The applied Radix-2 FFT algorithm [19] has a fixed struc- 
ture and does not change with measurement window variation. The BRAM count in- 
creases linearly and is proportional to the measurement window length. The BRAM 
components are used as data buffer only and not applied for calculation. The memory 
increases by expanding the measurement window length and vice versa. 

The Radix-2 algorithm uses a fixed amount of flip-flops (FFs) and lookup tables 
(LUTs) for FFT. Hence, measurement window length variation results in a logarithmic 
relation of logic resource consumption (see Fig. 3 (b)). A larger measurement window 
size requires a larger address space and additional glue logic, but the growth rate for 
buffer control logic consumption reduces with the increasing buffer size. 

As a result, the processing time increases linearly with the measurement window 
size (see Fig. 4). With a logic utilization less than 10 % on a Zedboard with Zynq®-7000 
All Programmable SoC XC7Z020-CLG484, the method is very resource saving. The 
amount of FFs and LUTs is logarithmically related to the number of samples per mea- 
surement window to be processed. However, the amount of used memory (BRAM) is 
very high and the reference cross-correlation architecture cannot be realized on small, 
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Fig. 3: Logic resources utilization for the reference design depicted in Fig. 2 on a Zedboard with Xil- 
inx Zynq®-7000 All Programmable SoC XC7Z020-CLG484. 
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Fig. 4: Linearly rising FPGA processing time for the 
0 20000 40000 60000 reference cross-correlation realization depicted in 
Measurement Window Size (Samples) Fig. 2 as a function of the measurement window size. 


low-budges FPGA devices with very limited BRAM resources. As the cross-correlation 
is just a single part of the overall solution, it shall not occupy the majority of BRAM 
resources of the host device. 

Especially the storage of additional data caused by linear spline interpolation con- 
sumes a lot of BRAM resources. As the interpolation data are artificial and do not repre- 
sent real measurements, a strategy to prevent the buffering of interpolation data will 
improve the BRAM utilization. While the reference FFT IP core buffers internally all 
data to process, the availability of interpolation data is essential for the convolution 
approach. Furthermore, the FFT core uses a measurement window size with a power of 
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two step width only. If fewer data are processed the remaining window share is filled 
with zeros, which causes unnecessary computational effort. The maximum lag can 
be determined for most of the TDOA use cases. Hence, it is not necessary to calculate 
the complete cross-correlation, which saves calculation time. The hereinafter highly 
scalable FPGA implementation addresses these disadvantages and faces the trade-off 
between resource limitation and measurement accuracy. The method uses fixed-point 
arithmetic and provides high up-sampling capabilities with linear spline calculation. 


4 Novel acceleration architecture 


The novel acceleration architecture for FPGA-based cross-correlation consists of 14 dif- 
ferent components (see Fig. 5). The first two parts are two True Dual Port Block-RAMs 
(BRAM X, BRAM Y) which buffer sampled input data that have to be processed. The 
first port is assigned for filling the BRAM and the second port provides access to the 
connected processing modules. The subsequent component is the data access (DA) 
block. This module controls input data readout from the BRAMs and handles inter- 
faces and address mapping. 


Fig. 5: Novel acceleration architecture for FPGA-based cross-correlation. 
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Next, the LSI block stores two subsequent values and calculates interpolation values 
depending on the interpolation factor &. The LSI block for x-values feeds a shadow 
register, which is organized as a shift register with a length of č. The interpolation 
block of the y-values feeds a second shift register with equal length č. The initialization 
value of LSI Y is for each interpolation step X„ -Tmax and generates k = 2T,,,,+1 values, 
while the next £ x-values are calculated simultaneously; T max describes the maximum 
delay (in clock cycles) of the cross-correlation given by Tmax = f (F; č) with sample 
rate F,. After each processing step of all k values, the shadow register X is copied to 
the processing register. The £ interpolation values in shift register Y are multiplied in 
parallel with all ë interpolation values in the processing register, which results in the 
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products w;(r) with {i € N |O <i < &}. All multiplication results w;(T) are finally 
summed up to Q(T) in the SUM module. 

In reference to equation (2), the interim result from the previous processing step 
for X,,_, is prepared by the BRAM control block, transferred to the SUM block, and ac- 
cumulated in the processing step of X,,. The BRAM control block stores the updated 
Xn in BRAM result and handles the storage access. The maximum detection index and 
the maximum detection value are stored by the MAX detection block after each sum- 
mation step. The procedure is repeated for all (N-£) x-values, where N is the amount of 
data per window. Since the multiplication generates big numbers and the range limit 
of integers is reached soon, the whole processing is done with fixed-point arithmetic. 
The width of the input data is defined as 32 bits, with 16 bits for the integer part and 
16 bits for the fractional part. 


5 Results 


The algorithm calculates for each x-value ë interpolation values and processes them 
in parallel. The computation is done 2Tmax + 1 times, which leads to a runtime of 
O(2NT max + 1), depending on the according sample rate. The runtime O is indepen- 
dent of the interpolation factor & since all interpolation values are calculated at the 
same point of time. 

Logic resource consumption of the proposed design depends on the value of in- 
terpolation factor &. FF, LUT, DSP slice and BRAM utilization is proportional to € with 
linear hardware complexity as shown in Fig. 6. Furthermore, both FF and LUT con- 
sumption are affected by the measurement window size and grow logarithmically with 
expanding window size and vice versa. The number of required DSP slices is inde- 
pendent of the measurement window length as all interpolation steps are executed in 
parallel. 

The BRAM utilization is linear and strongly linked to the runtime and changes 
with sample rate, measurement windows size and interpolation factor modifications 
(see Fig. 7). Window size and sample rate changes affect the two data input buffers 
(BRAM X/Y) and the BRAM for results, while a modification of ë affects the result 
BRAM only. The growth of BRAM consumption is moderate as existing memory uti- 
lization rises up to 100 % before a new memory block will be activated. The proposed 
architecture uses significantly less BRAM resources compared to the FFT IP core-based 
reference design. E. g., a hardware realization with an interpolation factor of č = 64 
requires only "60 % more FFs, "10 % more LUTs, and five times more DSP slices, but 
saves more than 97 % of BRAM resources for 65k samples (maximum of reference de- 
sign). With 18 FFs (0.3 %) and 60 LUTs (2%) more, the proposed design is able to pro- 
cess 262k samples (equal number of DSP slices and four times the BRAM consump- 
tion), which equals four times the performance compared to the reference design. The 
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proposed design does not store intermediate multiplication results, which leads to a 
significant reduction of BRAM resource consumption compared to the reference de- 
sign. Furthermore, the measurement window length and sample rate do not affect the 
overall processing part of the design. FF and LUT consumption is nearly constant and 
rises minimally for glue logic only, which is negligible compared to the available FF 
and LUT resources on modern FPGA devices. The independence of logic resource con- 
sumption from different measurement configurations leads to a beneficial scalability 
of the overall architecture. Additionally, the novel method provides a measurement 
window configuration of arbitrary size, which is only limited by the available logic 
resources of the applied FPGA. 


6 Conclusion and future work 


The proposed FPGA design provides a highly scalable FPGA architecture for hardware- 
accelerated cross-correlation with up-sampling support. The novel approach has a 
runtime of O(2NT nax) and a memory complexity of O(2N +32 +2T ax). The up-sampling 
is realized by LSI and characterized by an interpolation factor č. The modification of 
the interpolation depth affects the consumption of DSP slices and FFs only and does 
not influence the overall processing time, which is solely coupled to window size and 
sample rate. 

The novel approach enables arbitrary and unlimited measurement window 
lengths, which allows an extension of the maximum measurement window size com- 
pared to the reference design. The proposed solution requires significantly less BRAM 
and LUT resources for both high sample rates and large-scale measurement win- 
dow configurations. The presented architecture occupies slightly more FFs, which is 
negligible with respect to the available FFs on modern FPGA devices. 

Future activities focus on the analysis of time restrictions to ensure real-time capa- 
bilities of the design. Furthermore, techniques for a significant reduction of DSP slices 
have to be developed. 
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Ahmed Fendri, Racem Jribi, Olfa Kanoun, and Hamadi Ghariani 
Interface circuit for oil quality assessment 
considering dielectric losses and stray 
capacitances 


Abstract: Dielectric spectroscopy is becoming a popular technique in addition to the 
chemical and optical methods for oil quality assessment. It is based on the measure- 
ment of the variation of the dielectric properties of oil corresponding to the change 
of its chemical and physical properties. In this investigation an interface circuit for 
the measurement of the dielectric properties of oil is proposed using the capacitance- 
to-voltage conversion technique. The circuit measures the variation of capacitance of 
the sensor. The output voltage of the circuit is a DC voltage proportional to the capaci- 
tance of the sensor with a resolution of 5 mV/pF. The differential technique was used to 
eliminate stray capacitances caused by the circuit. Regression analysis has been car- 
ried out to describe a mathematical relation between the change of the capacitance 
and the output voltage. Simulations show a maximum relative deviation of 0.03 % for 
capacitance in the range of 100 pF to 1000 pF while experimental results show a max- 
imum relative deviation of 2% for the same capacitance range. The effects of stray 
capacitances and dielectric losses are also investigated. Stray capacitances decrease 
the measurement accuracy to 0.1%, while the parallel resistance, which represents 
the dielectric losses does not affect the accuracy, if its value is higher than 10 MO. For 
resistance values lower than 10 MO the circuit is not stable and it is threfore not pos- 
sible to estimate the capacitance value using the output voltage. 


Keywords: dielectric spectroscopy, capacitive sensor, capacitance measurement, 
capacitance-to-voltage conversion 


1 Introduction 


The measurement of small capacitance changes is very important in industrial appli- 
cations and in particularly for the measurement of dielectric properties of liquids such 
as oils. In fact, the measurement of a sensor capacitance dipped in oil can be used to 
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measure the dielectric properties of the oil. Many investigations used this technique 
to characterize edible oils [1], detect the adulteration [2] and estimate the water con- 
tent [3] or the total polar components in fried oils [4]. A typical measurement setup 
consists of a capacitive sensor connected to an impedance analyzer or an LCR meter, 
which measures the capacitance of the sensor and thus deduces the dielectric proper- 
ties. Such a measurement setup is accurate and simple to use but it has some limita- 
tion since it is not designed for online or mobile applications but only for laboratory 
measurement. 

Since new investigation are not focusing only in laboratory conditions but also 
in online industrial applications where the ideal conditions are not guaranteed and 
a smaller mobile setup with low energy consumption is required. It consists of an 
embedded measurement circuit connected to a well designed capacitive sensor. The 
measurement circuit converts the capacitance to measurable data (voltage, frequency, 
time, etc.), while the sensor is used in such a way that the change in the dielectric 
properties will cause a change in the capacitance. 

The capacitance measurement circuits can be classified into three main categories 
according to the type of conversion the circuit does: 

— capacitance to time, 
— capacitance to frequency, 
— capacitance to voltage. 


The capacitance-to-time conversion includes the charging/discharging and phase 
shift techniques. The working principle of the charging/discharging technique is sim- 
ple; the capacitance to be measured is charged during a known period of time and 
then the required time to discharge the capacitor is used to measure its value. It is also 
possible to measure the capacitance during the charging process by measuring the 
required time to charge the capacitor to a certain voltage. CMOS switches are usually 
used as switches to charge and discharge the capacitor; in fact the use of switches 
increases the stray capacitances and hence the measurement errors [5, 6]. The second 
technique is the capacitance-to-phase angle conversion; one of the capacitor elec- 
trodes is connected to an AC voltage source while the second electrode is connected 
to acharge amplifier and a compensating capacitor excited by a compensating signal. 
The use of this compensation permits a gradual phase shift in response to the change 
of the capacitance. This enables the circuit to measure very small capacitance changes 
but in a very limited range. The phase angle of the output signal of the charge ampli- 
fier is compared to that of the input signal and the difference in phase is converted to 
time using a Schmitt trigger [7]. 

Another techniques to measure capacitance is to convert its value to a periodic 
signal. Oscillation and resonance circuits are the main two types of the capacitance- 
to-frequency conversion technique. The oscillator produces a periodic signal with a 
frequency proportional to the capacitance C, [8]. The resonance technique is based 
on the use of a parallel or serial inductance and resistance to the capacitance to be 
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measured in order to create a resonator. The circuit is excited by an AC voltage with 
a variable frequency; the voltage and the current in the circuit should be measured 
at each frequency until the resonance happens where the current or the voltage is 
maximum, in case of a serial or parallel circuit, respectively [9]. 

The capacitance-to-voltage type includes mainly two techniques of conversion: 
the transducer and the I/Q modulation. The transducer circuit is simple and fast; it 
consists basically of a current-to-voltage converter, a rectifier and a low-pass filter in 
such a way that it provides a DC output voltage with an amplitude depending on the 
value of the capacitance of the sensor [10, 11]. This type of circuit does not consider 
the effect of dielectric losses existing in capacitive sensors. The I/Q modulation circuit 
considers the capacitive sensors as a complex impedance which includes not only a 
capacitance C but also a parallel resistance R representing the dielectric losses. In an 
I/Q modulation circuit, the sensor is fed by an AC voltage, the currents I. and I, flow- 
ing though the capacitance and the resistance, being converted to voltage and sent to 
two phase sensitive detectors which compare this voltage with two 90° phased volt- 
ages. The two output voltages of the phase sensitive detectors are proportional to the 
values of C and R. This type of circuit is complex and the accuracy of the measure- 
ment is influenced by multiple factors, such as the accuracy of the phase shifters and 
the phase sensitive detectors [12]. Another type of circuit uses only the Q modulation 
in order to compensate the effect of the parallel resistance. A feedback controls the 
phase of the input signal in such a way it compensates the phase shift caused by the 
sensor resistance [13]. 

The use of the capacitance-to-voltage type of circuit is also to be considered; 
in this type of circuit only an accuracy of mV is required, while the capacitance- 
to-time type requires nanosecond accuracy and the capacitance-to-frequency type 
requires MHz accuracy. Among the capacitance-to-voltage types the transducer cir- 
cuit is selected because of its simplicity and its high accuracy to measure capacitance 
changes compared to the I/Q type. The circuit proposed by Preethichandra in [10] is a 
capacitance-to-voltage type of circuit. A reference capacitance is used to compensate 
the effect of stray capacitance caused by the measurement circuit. It consists of two 
similar sub-circuits, one for the capacitance to measured C, and one for the refer- 
ence capacitance C,er. Each sub-circuit consists of a current-to-voltage converter, a 
rectifier and a low-pass filter. The difference between the two output signals is calcu- 
lated with a differential amplifier. The final output voltage V ut depends on C, and 
Cret. Since Cyor is kept constant, Vout depends in reality only on C,. The bias voltages 
caused by stray capacitances in both sub-circuits are eliminated by the differential 
amplifier. 

In this chapter, we study the possibility of using the capacitance-to-voltage type of 
circuit for the measurement of oil quality considering the effects of stray capacitances 
and dielectric losses. 
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2 Measurement circuit 


A simplified schematic of the measurement circuit is presented in Fig. 1. The working 
principle is to inject a sinusoidal signal to the sensor and a reference capacitance; 
this creates currents flowing through these two capacitances. These currents are then 
converted to voltage, rectified and filtered. The filters outputs two DC signals depend- 
ing on the sensor capacitance and the reference capacitance. The final output of the 
circuit, Vor, is the difference between both signals; Vut is maximum when C, is min- 
imum and it is theoretically O when C, = Cer. The circuit was adapted to oil measure- 
ment by adding a parallel resistance R, to C, and stray capacitances (C,,, C„ and C,;) 
connected to C, and Cyor. 


Active Rectifier 
LPF 


V2 

Fig. 1: Capacitance measurement 
circuit considering stray capaci- 
tances and dielectric losses. 


2.1 Active rectifier 


A precise full-wave rectifier is used in this circuit. This type of rectifier is considered 
more precise comparing with the half-wave rectifier, which keeps only half of the orig- 
inal input signal. In fact the full-wave rectifier reverses the negative part of the input 
signal and combines it with the positive part. The average (DC) output voltage is higher 
than for the half-wave rectifier and has much less ripples, producing a smoother out- 
put waveform. 

The input amplifier (Fig. 2) controls the biasing of the diodes D1 and D2 to change 
the signal path based on the polarity of the input signal. Diode D1 provides the positive 
half wave to the output amplifier while the diode D2 provides the negative half. The 
negative part is inverted and summed to the positive half at the output amplifier. 
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Fig. 2: Full-wave precision rectifier. 


2.2 Low-pass filter 


A low-pass filter is used to eliminate the ripples and high-frequency noises. A second- 
order filter based on Sallen-Key topology is used (Fig. 3). At low frequencies the ca- 
pacitors C, and C, are considered open circuit and the signal is buffered directly to 
the output. At high frequencies C, and C, cannot be considered as open circuit and 
the input signal is shunted to the ground. The cut-off frequency of the filter can be 
calculated as follows: 


1 


~ RRC 


A low cut-off frequency can guarantee a smooth DC output, but this will decrease the 
quality factor Q and the gain of the filter. 


fe (1) 


C2 


Fig. 3: Sallen-Key low-pass filter. 


2.3 Differential amplifier 


The function of the differential amplifier is to calculate the difference between the two 
output voltages of the low-pass filters, which represent the measured capacitance C, 
and the reference capacitance C,,r. 

The offset voltage caused by the non-ideality of the amplifiers should be low so 
that the output voltage is O when C, = Cyr. It is possible to measure capacitances 
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higher than the reference capacitance C,.r but only if the differential amplifier is sup- 
plied with +Vcc. In the case the amplifier is supplied with +Vcc and the ground, it is 
only possible to measure capacitances lower than Cyor. 


3 Simulation results 


The circuit blocks were first separately implemented in the simulator in order to study 
their behaviors and select the values of their components. The full circuit is then im- 
plemented for further investigation. The current flowing through C, is converted to 
voltage across the resistance R, of the rectifier (R, in Fig. 2). This voltage is then rec- 
tified and filtered by the full-wave rectifier and the low-pass filter; C, and R, forma 
high-pass filter with a variable gain and cut-off frequency changing with the change 
of C, and R,. For frequencies higher than the cut-off frequency the gain of the filter is 
equal to 1; in this case it is not possible to detect the change of the capacitance because 
the amplitude remains constant, but for frequencies lower than the cut-off frequency 
the gain will depend on the value of C, and the output amplitude will be changing 
with respect to the change of C,. The gain can be calculated as follows: 


(2) 


Gain = peeve 
VR + oo 

The sensitivity of the filter gain to the change of the capacitance is studied. The 
input frequency f,, and the resistance R, are varied in the ranges of 100 Hz to 100 kHz 
and 1kQ to 10 MO. The maximum frequency is set at 100 kHz due to the limitation of 
the amplifiers. In fact the amplifiers to be used in such measurement circuits should 
have a low offset voltage and low bias current, though amplifiers which have a low 
offset voltage, in the range of pV, and bias current, in the range of pA, have very limited 
bandwidth in the range of 100-1000 kHz. 

The maximum sensitivity (0.08 gain change for 100 pF C, change) was found at 
different combinations of f,, and R,. The input frequency fip is selected to be 10 kHz 
(much smaller than 100 kHz), corresponding to R; = 20kQ. 

In a first investigation, the circuit is simulated without stray capacitance and 
dielectric losses. The reference capacitance C, is kept constant at 500 pF and the 
voltages V,, V, and Vout are recorded (Fig. 1). The voltage V, changes linearly with 
the change of the capacitance C,, thus the voltage V ut changes also linearly with C, 
(Fig. 4). The simulations were done using the ADA4522 amplifier, which has a maxi- 
mum offset voltage equal to 5 uV and a maximum input bias current equal to 50 pA. 

Ideally the output voltage equals to 0 V when C, = Cyer, but due to the non-ideality 
of the amplifiers there is a voltage offset equal to 0.013 V at the output voltage. The 
relation between C, and V,,,; is linear with a resolution of 2 mV/1 pF. A curve fitting 
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Fig. 4: Output voltage of the measurement circuit. 


was used in order to find the relation between C, and Vout following the following 
equation: 


Cearc = 0-495 Vout + 492.3, (3) 


where C..ı. is given in pF and Vut in mV. 
To evaluate the goodness of the fit, the relative deviation was calculated using the 
following equations: 


Deviation = |Cy — Cealcl (4) 


Relative deviation = — ij wo. (5) 
X 


The relative deviation decreases with the increase of the capacitance C,; it has its 
maximum value 0.1% at 1pF and decreases to 0.02% at 1000 pF, which means high 
linearity between C, and V,,, in the ideal case. In fact, the assessment of oil quality 
using capacitance sensors suffers from two main stray effects: 

— the stray capacitances caused mainly by the sensor geometry and the interface 
between the sensor and the measurement circuit; 
— the parallel resistance which represents the dielectric losses. 


The effects of stray capacitances C,,, C,, and C,, were first studied (Fig. 1). The values 
of C,,, C, and C,; were set at 10 pF, 50 pF and 100 pF. 

A total of nine simulations have been made and the capacitance C,,),, the devi- 
ation and the relative deviation were calculated using equations (3), (4) and (5). The 
capacitance C,; has minor effects on the measurement since it is connected in paral- 
lel to the voltage source, which has a low input resistance. The stray capacitances C,, 
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and C,, also have minor effects when their values are equal; this is because the error 
in voltage caused by C,, is compensated by the error caused by Co (Vout = Va - Vj). 
The deviation in the former case was less than 1% in all cases (C,, = C, = 10, 50 or 
100 pF). The deviation caused by C,, + C,, is much higher than in the previous cases. 
It increases to a maximum of 8.5 % at 100 pF, decreases significantly around the refer- 
ence capacitance and increases again up to 2.5% at 1000 pF (Fig. 5). The importance 
of the reference sub-circuit is to compensate the stray effects. Using the reference ca- 
pacitance it is possible to guarantee a maximum deviation of 1% in the range of 400- 
600 pF in the presence of 100 pF stray capacitance. 


— Stray capacitance: Cai = 100pF and C2 = 100pF 
— Stray capacitance: Cai = 100pF and C2 =50pF 
TE —— Stray capacitance: Cai = 100pF and Cc. = 10pF 


—— Without stray capacitance 
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Fig. 5: Relative deviation caused by the stray capacitances. 


The second parameter to be considered is the parallel resistance to the sensor capaci- 
tance; this resistance R, represents the dielectric losses of the dielectric medium under 
investigation. The value of this parallel resistance varies depending on the sensor ge- 
ometry and the type of the medium. In the case of edible oils the value of R, varies 
depending on the type and quality of the oil. A typical value of R, is in the MQ range. 
For this purpose, a parallel resistance is added to the capacitance C, as illustrated in 
Fig. 1. The circuit is simulated with and without stray capacitances for resistance R, 
in the range of 1-100 MQ. 

With the presence of the parallel resistance, C,, R and R, form a high-pass filter 
expecting that the gain is constant at low frequencies (Fig. 6). We have 
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Fig. 6: Magnitude Bode plot of the filter formed by C,, R, and R,. 


The new transfer function depends on C,, R, and R, (equation (6)); at low frequencies 
(w = 0), H(w) = ree while at high frequencies (w = 00), H(w) = 1. The gain at low 
frequencies equals ie which is different from the case without R,. 

In the case Ry = 100kQ, the relative deviation is higher than 100 % and it was 
not possible to calculate C, using equation (3). In this case the current flows through 
R, only. This makes the circuit behave like an inverting amplifier with a constant gain 
independent of the value of C,. For higher values of R, the relative deviation decreases; 
at R, = 10 MQ it has its maximum (1%) at 100 pF while it has its minimum (0.6 %) for 
values near the reference capacitance (Fig. 7). In the case where both stray capacitance 
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Fig. 7: Relative deviation caused by the parallel resistance. 
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and parallel resistance are added to the circuit, the deviation effect is more by the 
parallel resistance; in fact, at a stray capacitance of 10 pF and when R, = 10MQ, the 
maximum deviation was 2% at 100 pF and it decreases for higher capacitance values. 


4 Experimental results 


The circuit was designed and implemented on a printed board for experimental inves- 
tigation. The amplifier ADA4522 was used in the circuit and it was selected because 
of its low offset voltage and bias current. Two amplifiers were used so that the two 
rectifiers and the two low-pass filters (measurement and reference branches) could be 
implemented in a way the error caused by stray capacitances of the measurement cir- 
cuit would be similar and hence compensated in the differential amplifier. The diode 
BASI6VY was used in the rectifier circuit; it was selected because of its very small re- 
serve recovery time (minimum time to switch from passing to blocking modes and vice 
versa). 

The measurement circuit is excited by a sine wave signal with a 2-Vpp amplitude 
at a frequency of 10 kHz. The waveform generator AD9933 was used because of its low 
cost, simple interface and stable output frequency and amplitude, especially at low 
frequencies. A micro-controller with an accurate ADC (12 bits resolution and 250 kS/S 
sampling rate) was used to measure the output voltage. 

A set of 30 capacitances in the range of 1-700 pF were accurately measured with 
a precision impedance analyzer for the calibration of the measurement circuit. The 
reference capacitance is chosen to be 500 pF; this limits the measurement range to 
500 pF, so that no negative voltage will be in the output of the circuit, which can dam- 
age the ADC of the micro-controller. 

An offset voltage of 60 mV was measured when C, = C,ef = 500 pF; this can be 
caused by the measurement circuit and the error of the analog-to-digital conversion. It 
was not possible to measure capacitances lower than 10 pF because the output voltage 
was not sensitive to the change of capacitance in this range, while it linearly decreases 
with the increase of capacitance in the range of 10-500 pF. The practical relation be- 
tween the output voltage and input capacitance was developed using the curve fitting 
tool, yielding 


Capacitance = 0.5304V ut + 468.3. (7) 


The maximum relative deviation is 30 % for capacitances smaller than 70 pF while 
it decreases to 2% at 100 pF and decreases further when the input capacitance in- 
creases. A minimum change of 3 pF was accurately detected, corresponding to a 
change of 6 mV in the output voltage. 
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5 Conclusion 


The proposed measurement circuit is suitable for the measurement of small capaci- 
tance changes in the pF range in the capacitance range from 10-500 pF. It was possible 
to reduce the effect of stray capacitances using the reference principle. The measure- 
ment accuracy increases for capacitance values near the reference capacitance and 
hence it is advised to used a reference capacitance equal to the nominal capacitance 
of the sensor. The dielectric losses have major effects on the measurement; in fact this 
circuit is not suitable for the measurement of high-loss materials, though it is possible 
to use this circuit for the assessment of oil quality since it has low dielectric losses. 
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Part Ill: Sensors 


Christian Weber, Markus Tahedl, and Olfa Kanoun 
Modeling for improved performance 
of non-contacting capacitive sensors 
for detecting aqueous solutions 


Abstract: In industrial applications, non-contacting capacitive sensors are used to de- 
tect conductive fluids in containers. During use, the material to be detected can stick 
to the inside of the container, leading to measurement deviations. Using analytical 
modeling, it has been found that an overall resistance can distinguish between the 
conductive film and the actual fill level. Impedance analyzers have limitations with 
respect to their maximum measurement frequency of approximately 100 MHz. This 
bandwidth is not sufficient for characterizing highly conductive media. In order to 
overcome this limitation, improved signal processing is proposed to determine the 
overall resistance. In this paper, two methods of parameter extraction using the total 
measured bandwidth and a limited frequency range are compared. Results show that 
the overall resistance can be extracted using only a limited frequency range. 


Keywords: impedance sensors, capacitive sensors, parameter extraction, modeling, 
impedance spectroscopy 


1 Introduction 


Limit levels of a container filled with a conductive fluid can be detected using ca- 
pacitive sensors. For easy installation, the sensors are mounted outside of the non- 
conductive container. Conductive films inside of the container influence the sensor 
signal and can produce a false positive detection of the medium [1]. 

Fig. 1 shows a capacitive limit level switch mounted outside a container. Since in 
many applications the exact geometry, i. e., container material and wall thickness, and 
material properties of the medium are unknown, the sensors have to be calibrated 
before use. This is often done by adjusting the switch point using a potentiometer. 
Simple capacitive proximity switches measure only the absolute impedance between 
their sensing electrode and ground potential [2, 3]. Depending on the measurement 
frequency, wall thickness and conductivity of the fluid, the measured value of such 
sensors is the same for both the actual fill level and a conductive film. 
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a. Fig. 1: Capacitive proximity switch mounted outside a 


container for detecting the limit level of the medium. 
medium Conductive films may stick to the inside of the con- 
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Fig. 2: Example plots to illustrate the difference between a thin layer and the actual full container for 
two different conductivities. 


To illustrate the problem, Fig. 2 shows example spectra for a full container and a thin 
conductive film of a conductive aqueous solution for two conductivities 2mS/cm and 
20 mS/cm. It is shown in Fig. 2 (a) that a thin conductive film cannot be distinguished 
from the actual fill level below a frequency of 2 MHz. A difference between a thin anda 
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thick layer can only be observed at frequencies above 2MHz. This effect is even worse 
at higher conductivities, as can be seen in Fig. 2 (b), where the difference between thin 
and thick layers is only observable at frequencies higher than 20 MHz. The necessary 
measurement frequency increases approximately proportionally with the conductiv- 
ity ofthe material to be detected. 

Fig. 3 shows the difference of the absolute impedance, A|Z| = |Zjninl — IZinicr)» Þe- 
tween a thin layer and a thick layer for several conductivities. It can be seen in the 
figure that the maximum difference decreases with increasing conductivity. Further- 
more, the measurement frequency at which the maximum difference occurs also in- 
creases with conductivity. 


Fig. 3: Absolute difference A|Z| = |Zjnin! — |Zthickl 
between a thin layer and a thick layer of the 
medium. Increasing the frequency reduces the 
observed difference and increases the necessary 
Frequency in Hz measurement frequency. 


Since measurement devices are limited in terms of measurement frequency, another 
method must be found. Figs. 2 (c) and 2 (d) show complex Nyquist plots of the same 
dataset. It is possible to distinguish thin films and actual limit levels using impedance 
spectroscopy and associated signal processing techniques. The low-frequency part, 
where the Nyquist plot forms a straight line, can be shifted to the right by a decrease 
in layer thickness or a decrease in conductivity. Therefore, for a given conductivity, dif- 
ferent layer thicknesses can be distinguished by inspection of the complex impedance 
plot even though there is no difference visible in the absolute value of the impedance. 

Using analytical modeling, it was found that the overall resistances of a capac- 
itively coupled aqueous solution can be used to distinguish between the conductive 
film and the actual limit level [1, 4]. To validate the proposed criterion experimentally, 
a measurement setup is needed to create thin layers of aqueous solutions. 

This setup has already been described in [5]. However, the fitting methods are 
reported to have some limitations. Relative fitting errors of the proposed fitting proce- 
dures are quite large, up to 20 % in the real part of the measured spectrum. 

In this contribution, improved fitting procedures are proposed to increase the 
quality of fit, which has been improved to below 0.5% and 7% for the imaginary 
and the real part, respectively. For very high conductivities, e. g., aqueous KCl solu- 
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tions with a conductivity of up to 208 mS/cm [6, 7], a linear fitting procedure [1] is 
introduced. In the next sections, we introduce and compare two fitting procedures. 

For low solution conductivity, the parameters of a behavioral equivalent circuit 
model are extracted using differential evolution. For high solution conductivity, low- 
frequency data are fitted to a first-order polynomial and the overall resistance is de- 
termined. 


2 Extraction of the overall resistance by non-linear 
fitting 


To adequately model the behavior of measured and simulated impedance spectra, the 
electrical equivalent circuit (EEC) shown in Fig. 4 can be used. 


CPE R; Ro 


o 


-A 4 Fig. 4: Behavioral model for the complex 
impedance of capacitively coupled conductive 
Cy Co media. 


Since the electrode is mounted outside a non-conducting container, the capacitive 
coupling into the material under test (MUT) can be modeled using a constant phase 
element (CPE). Its impedance can be parameterized as follows [8]: 


ZcPE = ž - jw)”, (1) 


where j is the imaginary unit and K denotes the amplitude of the CPE; n can be any 
real value between -1 and 1 and is a measure for the dielectric loss factor D. We have 


_ RZ) 


0 2) 


If n = 0, the CPE becomes an Ohmic resistance R = $. For n = -1, the equation 
yields an ideal capacitance C = K; n = 1 yields an inductance L = 2. 
To link n and D, we rewrite equation (1) using [8] and we obtain 


Zcpr = . - (cos(n5 } +j-sin(nZ)), (3) 


Substituting into equation (2) yields 


A cos(n5) RR 


(4) 


~ sin(n3) u tan(nZ) 


Modeling for improved performance of non-contacting capacitive sensors — 111 


The equations for the CPE phase angle p and quality factor Q are then 
p=n-, (5) 
1 1 
07, = tan( nZ ). (6) 


If an alternating electric potential is applied to an electrode, the equation of con- 
tinuity of the electric current density states that the current can be separated into a 
dielectric and a conductive part [9]. Thus, the behavior of a conductive aqueous solu- 
tion can be described by a series of resistor-capacitor (RC) elements. The impedance 
of a single RC element is 


R 
= $ 7 
RC 1+jwRC 2 
Substituting w, = x yields 
R 
eS ae (8) 
Re ay 2, 
Thus, the model impedance Z 10a is 
Ty = cnc Rı R 
Zmod = =: UW) + + ; 9 
mod K G ) ja 1472 (9) 


gr 


To extract the model parameters of the EEC shown in Fig. 4, a differential evolu- 
tion algorithm as proposed in [5] was used. Evolutionary algorithms benefit from a 
large search space and robust convergence behavior [10, 11]. Compared to gradient 
methods, which tend to get trapped in local minima of the cost function, especially at 
increased noise levels, stochastic methods are more robust because no gradients are 
calculated [12], though it is worth mentioning that the noise level of an optimization 
problem also depends on the cost function employed [13]. 

The first two points are of special interest here, since the impedance spectra of 
the EEC in question require a high-frequency bandwidth for all parameters to be ex- 
tracted [5]. 

For the fitting procedure to succeed, a suitable cost function must first be defined. 
Since the measured data represent a capacitively coupled conductive fluid, the phase 
angle of the measured impedance is always relatively close to —90”. It is therefore not 
recommended to use the absolute difference between measured and modeled data 
as a cost function since it is mainly sensitive to changes in the imaginary part of the 
spectrum. Instead, we propose to adopt a parametric cost function c as proposed by 
[14] to describe the modeling error for N data points and M parameters, i. e., 


1 u ( (R(Z; meas) = R(Zimoa))* (J(Zi meas) = Mai 
c= Y + 


(10) 
R(Zi meas)? ner 
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This cost function describes the sum of relative quadratic residuals of the mea- 
sured data. This way, the real and imaginary residuals at all frequencies contribute 
equally to the cost function, which allows for a better extraction of the unknown 
parameters. When using differential evolution, extracted parameters can depend on 
the parameter boundaries specified. Since model parameters must be found over a 
large frequency bandwidth and the algorithm should work on as many different spec- 
tra as possible, parameter boundaries also become very broad, ranging over several 
decades. Thus, we propose to employ logarithmic compression to shrink the search 
space, i.e., 


R 
Riog = logio( 10 ) (11) 


For example, if the search space for a single resistance ranges from 10 to 1MQ, 
one substitutes both of these values into equation (11), which yields unitless values of 
O and 6, respectively. 

At low frequencies, the real and imaginary parts of an impedance spectrum of a 
capacitively coupled aqueous solution have a high sensitivity with respect to the CPE 
parameters [5]. Small perturbations in these parameters will cause a large fitting error. 
To guarantee a reliable extraction of the CPE parameters, we suggest to first estimate 
the CPE parameters from the measured data and then choose the initial population 
boundaries to be within +50 % of the estimated values. The CPE values can be esti- 
mated by linearly fitting the low-frequency part of the spectrum to a first-order poly- 
nomial, i.e., 


Fig. 5 shows an illustration of the linear fitting procedure. For a successful esti- 
mation, the part of the Nyquist plot where the impedance curve forms a straight line 
must be used. The imaginary part of the spectrum is used as an independent variable; 
Ro is the intercept point with the real axis at J(Z) = 0. The loss factor D; represents 
the slope of the low-frequency line. 


—3(Z) inka 


o 10 20 30 Ao 50 60 70 Fig. 5: Illustration of the estimation of overall 
Ro RZ) in kQ resistance and CPE parameters. 
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Using equation (5), the exponent n,» of CPE can be approximated; K can be approxi- 
mated by rearranging equation (1), i. e., 


1 


zen ae = 


where Z(w,) is the measured impedance at the lowest frequency. With the estimated 
CPE parameters and logarithmic compression for all RC elements, the initial popula- 
tion is computed using latin hypercube sampling [15]. After convergence of the dif- 
ferential evolution optimizer, the implementation [15] employed in this work can op- 
tionally refine the stochastic result with a direct search optimizer [16, 17]. If the direct 
search optimizer yields, in terms of cost, a better result than the evolution algorithm, 
its solution is replaced. If not, the result of the differential evolution is returned. 


2.1 Evaluation of the proposed modifications 


To evaluate the proposed method, a noiseless test dataset using the EEC shown in 
Fig. 4 is generated in a frequency range from 1 kHz to 100 MHz with the nominal val- 
ues shown in Table 1. The modified extraction method with logarithmic compression 
and approximation of the CPE parameters was compared against a classical differen- 
tial evolution approach. For both methods the algorithm was set up with the same 
population size of 80 individuals, 10‘ maximum iterations and the same convergence 
criterion described in [15]. The algorithm terminates if 


stdev (population) -3 


1 
mean (population) (14) 


is true, where stdev (population) denotes the standard deviation of the population. 
All other parameters of the optimizer were left at the standard values provided by its 
implementation in the Python module Scipy 0.18 [15]. 

In the implementation used in this study, C, and C, are also treated as CPEs to 
make the algorithm easily expandable to more complex models, though their expo- 
nents n are chosen in such a way that their loss factor does not exceed 1.6 x 10°”. 


Table 1: Result of the parameter reconstruction. 


Parameter Nominal value Modified reconstruction Classical reconstruction 
n -0.9923 -0.9923 -0.9923 
Kin 424)” . 10°? 4 4 4 
C, in pF 150 149.99 104 
C, in pF 14.5 14.5 13.2 
R, 500.00 500.00 999633.70 
R3 6600 6599.99 7086.71 


Cost function c 0 2.61 x 102° 1.01 x 107° 
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The boundary values for the search space of the resistances is [0, 6] (logarithmic) 
for the modified method and [10,1MQ] for the classical method. For capacitances, 
intervals of [-13,-8] (logarithmic) and [10 °F, 1078F] were chosen. For the coupling 
CPE, the modified method employs boundary values of [0.5xK,p; 1.5xK,p] and [-1, 0.9x 
Nap], While the classical method uses [107 aaa, 1078 a caar] for the CPE ampli- 
tude and [-1,-0.8] for the CPE exponent n. Implementation was done using [15, 18]. 

Table 1 shows that a reconstruction of the nominal parameters is possible with 
the proposed method. With the classical method, the algorithm seems to converge to 
a local minimum, as can be seen in the rightmost column of the table. 

Furthermore, the result of the classical method is not always reproducible. Mul- 
tiple optimization runs yield different parameter vectors with cost function values, 
which in some cases are as high as c = 0.3. 

For the classical method to converge to the correct solution, the parameter bound- 
aries have to be narrowed by two decades to [10 "F, 10 ?F] for the capacitances and 
[10,10 KQ] for the resistances. This shows that using logarithmic compression in- 
creases conversion reliability and user convenience. 

However, at the chosen population size of 80 individuals, using linear fitting to 
estimate the CPE parameters does not seem to be necessary as long as logarithmic 
compression is used. The evolution algorithm also converges on the nominal pa- 
rameter values even if the parameter boundaries for the coupling CPE are set to be 
107? 424)", 10°84 (24)"] and [-1, -0.8]. 

Although the linear parameter estimation does not yield significant improvement 
to the non-linear extraction, it is still useful to characterize solutions where the mea- 
sured spectrum does not contain enough information about all RC elements. 


3 Extraction of the overall resistance by linear fitting 


Because of the limited measurement bandwidth, it is not always possible to extract 
all time constants of a given measured spectrum. However, one can still obtain in- 
formation about the overall resistance and the coupling CPE. Instead of fitting the 
impedance spectrum over its full bandwidth, one can use equation (12), which is used 
to extract CPE parameters before a non-linear regression. It has already been shown 
experimentally in a previous work that this parameter is related to layer thickness and 
conductivity of the fluid [1]. Using a behavioral model with a CPE and one RC element 
with time constant T = R,C,, it can be shown that Ro is an approximation of R,. Most 
symbolic calculations have been carried out using [19]. Separation of real and imagi- 
nary part of the impedance yields 


1 2p2, n+2 (=) n (=)) 
R(Z) = — | C;Rjw’ < cos| — |+ KR, +w cos| — } |, 15 
= Kerr oo 2 : 2 15) 
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1 2,2. 2, (AM 2 noi. [m 
Z) = = CUR — )-CıKR — |. 1 
AZ) Ren 1w sin( 5 ) C KR w + w sin( 5 )) (16) 


In the frequency region where the Nyquist plot forms a straight line, two frequen- 
cies w, and w, are selected. The line slope D; can be approximated as follows: 


R(Z(w,)) - R(Z(w,)) 


Dy(W,, w2) = ; (17) 
rv) = FZ) = E(w) 
By substituting equation (17) into equation (12) and solving for Ro we get 
Ro(W 1, w2) = R(Z(w,)) - Delw, w) - I(Z(w)). (18) 


We now let the difference between the two frequencies approach zero by calculat- 
ing the limit, i.e., 


Ro(W2) = lim Row, w2), (19) 
1 2 
and we get 


CR 
Ra +1) 


Ro(w;) = (20) 
CIRIw; + 
For simplicity, we also substitute RC} = a, so we obtain 
g 
R; w 
Ro(w,) = ( - + 1), (21) 
I 241 w,tan() 


s 


Equation (21) has a pole for n = 0. However, this pole is not relevant in the appli- 
cation discussed here, because the CPE becomes a pure resistance in that case, thus 
removing the capacitive coupling from the system; Ro is a good approximation of R4, 


w? 
if wy > w, holds. In this case, the terms > and my both approach zero and the 


Wg an 
equation becomes Ro = R4. Typically, w, k w is fulfilled for highly conductive media 


since the low resistances cause high cut-off frequencies. 
It can also be shown that D; is a good approximation for the impedance loss factor 
D. The loss factor of an EEC containing one RC element and one CPE is 


2am an m 
WwW cos(— w = cos (= 
KR, Gy Œ) 


2 2 
rae a a 
D = en P am it af. am * (22) 
KR,w ww sin (>) w sin (>) 
2 2 2 
ah) Wile) mar 


If W, >W holds, we can calculate the limit, i. e., 


—n 
D; = lim D= oom (RR, +w" cos (=) (23) 


ay” sin ‘ee 
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which we can write as 


= nit 
w "KR, cos(= 


E sin) 

The denominator of the first additive term of equation (24) is close to -1 since 
n is also close to -1 for a capacitively coupled medium. Also, w" = w holds for 
the same reason. However, since K is typically small, approximately 10 2 1 (Mayr 
to 1010 1.radın, it can be stated that for small values of R, the first additive term 
approaches zero. Thus, D; approaches D. 

Equations (21) and (24) show that a linear fitting routine can be used to extract the 
two CPE parameters and the overall resistance Ro, which is closely related to the ac- 
tual resistance of the fluid. Therefore, even though the bandwidth of the measurement 
setup is limited, it can still be used to characterize highly conductive media. 


4 Experimental investigations 


The proposed fitting methods are now validated using measured data. The media and 
conductivities used are listed in Table 2. 


Table 2: Conductivity of the measured media. 


Medium Conductivity o inmS/cm 
Tap water 0.5 
KCI #1 1.9 
KCI #2 10.6 
KCI #3 43.1 
KCI #4 92.5 


KCI #5 193.1 


Fig. 6 shows a concept drawing ofthe proposed measurement setup. It consists ofan 
outer and an inner cylinder fabricated from plastic. The inner cylinder is hollowed 
out to minimize its influence of the measured impedance. Its height can be adjusted 
to create a defined gap at the bottom of the setup, where the measurement electrode 
is mounted. The height of this gap will be referred to as layer thickness. The aque- 
ous solution to be measured is inserted between the outer and the inner cylinder. 
A detailed description is available in [2]. Impedance spectra were recorded at 801 fre- 
quency points, the maximum amount for a single sweep, in a range from 4 kHz to 
110 MHz. 

Measurement noise of the impedance measurement instrument depends on var- 
ious factors, i.e., the connected device under test, the selected frequency range and 
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Fig. 6: Measurement setup drawings [5]. The setup is connected to an impedance measurement 
instrument, which provides high, low and ground potential terminals. 


the test signal level [3]. Measurement noise for capacitively coupled resistances char- 
acterized in this work was approximated by connecting an inter-digital capacitor on 
an FR4 printed circuit board with CPE values K = 5.7 x 10? 1 man andn = -0.99 
connected in series to a 150-Q resistor. A total of 600 spectra with 50 points each in a 
frequency range from 1 MHz to 5 MHz were recorded. The maximum relative standard 
deviations for the absolute impedance and its phase angle are 0.07 % and 0.007 %, 
respectively. 

Measured data are smoothed using a linear fitting algorithm [21] with the modi- 
fications described in [1, 22, 23]. For all models, the relative residuals for real A,, and 
imaginary A,, part of the impedance are calculated as follows: 


4 R(Zmeas) = R(Zmoa) 


- 100 %, 25 
Rees) A = 
IZmeas) = AZ mod) 
Aim = - 100%, (26) 
ii Zus) i 


where Zmeas denotes the measured impedance while Z,,,q is the modeled impedance. 
A measured impedance is excluded from the spectrum, if its relative residual in either 
the real or the imaginary part exceeds 0.5%. However, even though the fitting pro- 
cedure itself is optimized, the selection of error thresholds and frequency bandwidth 
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must still be done by the user. Therefore, different media can have different frequency 
ranges. 


4.1 Extraction by non-linear fitting 


The filtered spectra for tap water and KC] #1 are then used to validate the improved 
fitting procedure. All other media are employed to validate the linear fitting method, 
since sensitivity analyses have shown [5] that the measurement bandwidth for high 
conductivities is insufficient to extract all parameters of the EEC. Fig. 7 shows the rel- 
ative residuals for tap water and KCl #1. A series inductance L, has been added to the 
EEC to model the influence of cables connected to the measurement setup. This sig- 
nificantly reduces the relative residual in the imaginary part for frequencies greater 
than 10 MHz. The relative error of the imaginary part, shown in Figs. 7 (d) and 7 (c), 
is now below 0.5 %. The extracted value does not change over varying layer thickness 


a < 
Frequency in Hz Frequency in Hz 
(a) Tap water real residual (b) KCI #1 real residual 
xL az 
c £ 
Ei £ 
l 0.25mm —0. —0. 75mm — 
— 0.5mm 1mm 
ch in Hz Frequency in Hz 
(c) Tap water imaginary residual (d) KCI #1 imaginary residual 


Fig. 7: Relative residuals after non-linear fitting for tap water and KCI #1. 
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and changes by about 10 nH, when the MUT conductivity is increased. This suggests 
that L, can be solely attributed to the measurement setup itself. 

The maximum residual observed for A,. is 7 %, which is a significant improvement 
in comparison to [5], where the maximum fitting error in the real part is up to 20%. 
Significant errors arise for frequencies below 100 kHz and above 10 MHz. At lower fre- 
quencies, the measured spectrum is dominated by the CPE and small changes to its 
parameters cause a significant change in the modeled impedance, as can be seen in 
Fig. 8. Also measured signals in this region are very small, causing a decrease in the 
measurement accuracy, especially in R(Zyeas) [20]. 


SS ET TY 


EE ee ed ee 
ian ee | 


— Model 


~3(Z) in kQ 
—3(Z) in kQ 


RZ) in kQ R(Z) in kQ 
(a) Tap water (b) KCI #1 


Fig. 8: Example plots of measured and fitted data for a layer thickness of 0.25 mm. For visibility rea- 
sons, only every third measurement point is shown. 


Since all frequencies have an equal weight in the employed cost function, measure- 
ment errors in this region contribute equally to the fitted parameters. If frequencies 
below 100 kHz are excluded, A,, is below 1% for frequencies between 100 kHz and 
10 MHz, as can be seen in Fig. 7 (a). Above 10 MHz high relative residuals are observed 
at single measurement points. These might be caused by electromagnetic coupling 
into the measurement setup. Since the data smoothing is based on a linear fit, it min- 
imizes the weighted quadratic sum of residuals, which might allow single points to 
pass through the smoothing scheme. To further increase accuracy, it is feasible to 
first extract parameters over a broad smoothed spectrum and then exclude frequency 
ranges with a high relative residual. 

Tables 3 and 4show the extracted values. As expected, the CPE parameters K andn 
do not depend on layer thickness and conductivity. Their maximum deviation is 1.72% 
and 0.08% and therefore they only represent the capacitive coupling into the MUT. 
Resistances decrease monotonously with increasing layer thickness and conductivity. 
Resistances decrease with layer thickness because of the increasing cell constant. For 
constant geometries the sum of R; and R, is in good approximation inversely propor- 
tional to the MUT conductivity; C} and C, increase monotonously with increasing layer 
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Table 3: Parameters extracted by non-linear fitting for KCl #1 for several layers. 


Parameter 0.25mm 0.5mm 0.75mm 1mm 5mm 


n -0.9928 -0.9927 -0.9926 -0.9926 -0.9927 
K 8.77 - 1071? 8.77 - 1071? 8.78 - 1071? 8.78 - 1071? 8.78.10"? 
C, in pF 62.0 65.1 71.1 82.9 142.0 
C in pF 13.0 13.8 15.0 15.3 17.1 
R, ind 165.72 126.44 97.41 73.24 51.76 
R ind 2982.90 2446.66 2076.06 1980.07 1651.90 
L; in nH 251 249 249 246 251 


Table 4: Parameters extracted by non-linear fitting for tap water for several layers. 


Parameter 0.25 mm 0.5 mm 0.75 mm 1mm 5 mm 


n -0.9923 -0.9922 -0.9923 -0.9922 -0.9924 
K 8.87 - 1071? 8.83 - 1071? 8.81 - 1071? 8.82 . 1071? 8.86 - 1071? 
C, in pF 77.8 77.2 80.5 81.9 167.0 
C in pF 12.7 13.5 14.8 15.2 17.1 
R, ind 319.28 290.71 234.67 228.05 84.55 
R ind 9868.40 8318.83 6892.72 6806.33 6060.37 
L, in nH 235 240 238 236 234 


thickness. This behavior can be expected because the mean permittivity increases 
when more water is present in the setup. 


4.2 Extraction by linear fitting 


In the previous section, broadband characterization has been used to extract all sys- 
tem parameters including the overall resistances. However, this approach fails for 
high conductivities because of insufficient sensitivity of the RC elements. The linear 
fitting procedure is independent of this sensitivity. Accuracy of the extracted CPE pa- 
rameters can also be increased using the linear fitting method, because it only utilizes 
a measurement bandwidth, where the impedance spectrum is dominated by the CPE 
and the overall resistance. 

Thus, for all other media the overall resistances are extracted using the linear fit- 
ting procedure. Results are shown in Table 5. Just like the resistances obtained with 
the non-linear fitting procedure, extracted values decrease monotonously with layer 
thickness and conductivity. Also, the sensitivity of Ro with respect to layer thickness 
decreases with increasing conductivity. This is consistent with modeling expectations 
formulated in [1]. For KCl #5, the sensitivity becomes small and therefore susceptible 
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Table 5: Extracted Ro values. 


Layer thickness Tap water KCl #1 KCl #2 KCl #3 KCL #4 KCI #5 
0.25mm 10465.03 3110.23 652.22 165.58 62.06 22.32 
0.5mm 8906.57 2562.32 484.37 115.95 52.42 28.17 
0.75mm 7436.25 2174.59 391.52 94.89 42.82 23.66 

1mm 7341.06 2065.30 377.92 89.09 40.56 22.32 

5mm 6442.48 1719.78 319.35 72.12 33.34 19.40 


to fitting errors caused by measurement noise. This can be circumvented by averaging 
over several spectra. 

For tap water and KC] #1 the extracted Rp values are within reasonable tolerances 
equal to the sum of R, and R, obtained by the non-linear fit. The maximum error for 
tap water is 4.8 % while the maximum error for KC] #1 is 1.2%. This is in line with the 
predictions made in Section 3. 

Also the relative change of the Rọ values for a constant layer thickness is approx- 
imately proportional to the change in conductivity, i.e., if the conductivity at 0.5mm 
decreases by 30% the Ro value increases by about 30 %. The mean deviation from 
this proportionality is 8.2% while the maximum deviation is about 33 % for KCl #5 at 
a layer thickness of 0.25mm. A possible cause for this is the mechanical tolerance of 
the measurement setup, especially at small layer thicknesses. 


4.3 Comparison of the two fitting methods 


Fig. 9 illustrates the difference between the two fitting procedures. The non-linear fit- 
ting method employs the entire available frequency range, the full width of the Nyquist 
plot, to extract the parameters of the EEC shown in Fig. 4. However, for high conduc- 
tivities the non-linear extraction fails, because of the short time constants formed by 
low resistances. 

Using linear fitting on a limited frequency range, where the Nyquist plot is a 
straight line, one can still obtain the CPE parameters and the overall resistances 
Rọ and use them to distinguish between a thin film and the actual fill level, even 
though the non-linear fitting method fails due to insufficient sensitivity with respect 
to the RC elements. 

In the previous sections, it has been shown mathematically and experimentally 
that Ro is related to the resistances of the EEC. 

Qualitative behavior of the extracted parameters is within theoretical expecta- 
tions. However, the linear fit is only valid if the frequencies used for fitting are much 
smaller than the cut-off frequencies formed by the conductive and capacitive part of 
the current. This is usually the case for highly conductive media, for which this method 
is intended. 
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5 Conclusion 


Improved techniques for data processing have been proposed for measurement setups 
used to characterize thin layers of aqueous solutions. In contrast to [2], where fitting 
was done using purely gradient-based or direct search methods, broadband fitting of 
measured data to an EEC is now done using differential evolution and logarithmic 
compression of the parameter space. 

Using differential evolution with logarithmic compression has been shown to 
increase reliability of conversion. The evolution algorithm can be called with larger 
search spaces, thus improving convenience for the user. 

Also, a series inductance is introduced into the model to improve quality of fit 
at frequencies greater than 10 MHz. The improved method yields maximum relative 
errors of 0.6 % and 7 % in the imaginary and real parts of the spectrum, respectively. 
This is a significant improvement with respect to [5], where relative errors in the real 
part of the impedance were reported to be up to 20 %. 

If characterization using linear fitting is sufficient, solutions with a conductivity 
of up to 190 mS/cm can be characterized. 
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Eddy current corrosion measurement of steel 


Abstract: Contactless eddy current sensors allow efficient detection and measurement 
of the corrosion of metals (like carbon steel), also under the paint layer, as demon- 
strated in the described experiments and analysis. The corrosion was estimated by a 
Dodd-Deeds eddy current air-core single-coil model, according to the variation of the 
impedance in the sensor coil. In the case of magnetic materials (like carbon steel) the 
sensor signal depends not only on the conductivity of the specimen, but also on the 
magnetic permeability. Furthermore, the sensor lift-off influences the results, so ap- 
propriate lift-off compensation schemes are beneficial. In the current study a planar 
sensor coil and another ferrite-core probe were used. The estimated resolution of the 
corrosion layer measurement is about 1 pm (or better). A PC-based application front- 
end for estimation of the corroded area was developed. 


Keywords: eddy current, NDT, Dodd—Deeds model, corrosion measurement, carbon 
steel 


1 Introduction 


Eddy current sensors are widely used for non-destructive testing and characterization 
of metals and other electrically conductive materials [1]. The analysis of the complex 
impedance of the single-air-core measurement coil as a probe above the metal under 
test has been proposed in the works of Dodd and Deeds from the 1960s [2, 3]. Alter- 
natively, two-coil measurement setups, where one coil is transmitting the AC signal 
and the other is receiving, have been proposed, in one case as the receiving sensor 
coil is on the same side of the metal plate (“reflection-type measurement”) [4] or on 
the other side of the metal plate under test (“through-the-metal setup”) [5]. Anyway, in 
the current work, the single-coil configuration [2] was used, as being simple and well 
matching with the theory in the wide frequency range (up to 10 MHz) [6-8]. In [9], a PC- 
based simulation tool has been demonstrated for solving the Dodd-Deeds model for 
one- or two-layer metal plates with measurement coil above the metal plate. An im- 
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proved, computationally faster “series expansions” (instead of infinite integration in 
the space) version of the single-coil model is also available [10]. The Dodd—Deeds 
model and other similar algorithms solve the “forward problem” - that means, how 
the coil impedance depends on the parameters of the measured object (e. g., metal 
plate) and the measurement setup (parameters of the measurement coil). But in real 
life the solving of the inverse problem is needed - calculation of the properties (elec- 
trical conductivity and/or magnetic permeability and sensor coil lift-off) of the object 
under test from the measured complex impedance of the measurement coil. The lift-off 
effect and compensation of it have been discussed in [11]. Full inverse solutions have 
been proposed, e. g., in [12, 13]. One possible problem in eddy current measurements 
is the dependency of the results not only on the electrical conductivity of the material 
under test, but also on the magnetic permeability [14]. 

Some ideas on how to solve both these variables have been shown in [15, 16]. In 
[17] it has been shown (for the carbon steel AISI 1018) that the magnetic permeability is 
a frequency-dependent complex-valued variable. The eddy current corrosion process 
investigation has been done in [18, 19]. In [18] ithas been shown that both the electrical 
conductivity and magnetic permeability of the corrosion area are much (by an order) 
smaller than the values for the carbon steel (for SS400 in their example). 


2 Research performed 


2.1 Measurement setup with air-core coil 


The single-coil measurement setup is shown in Fig.1. The setup includes the pla- 
nar (PCB) measurement coil (shown in Fig. 2) with internal and external diameters 
D1 (2.5mm) and D2 (16mm), height h (20 um), number of turns N = 50 and lift-off z 
and the metal (steel) plate with thickness of d (about 2mm in the described exper- 
iments). The expected electrical conductivity of the material is about 6MS/m. The 
impedance measurement instrument Wayne-Kerr 6500B has been used to get the coil 
impedance and a specially designed 3D scanner moving the probe (coil) has been de- 
veloped (Fig. 3). Also the Windows-based software has been developed for reverse 
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Object under test Fig. 1: Single-coil measurement setup. 
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Fig. 2: Measurement coil — a photo (on the 
left) and PCB layout. 


Fig. 3: Developed eddy current 
scanner. 


eddy current modeling, using a combination of the forward model implemented as 
C/C++ code in the form of DLL and R-language for reverse solving, analysis and plot- 
ting of the features of the impedance measurement results. The single coil above the 
single-layer metal with variable lift-off and electrical conductivity and/or magnetic 
permeability approach has been used in this modeling and analysis software. 


2.2 Experiments with initial specimen 


Allspecimens used in the current research have been produced from the steel EN 10130 
class DCO1 (similar to the American AISI 1008 grade). On the first specimen (Fig. 4) 
some round-shape areas were corroded and the whole plate was painted afterwards. 
The region no. 2 was scanned by the measurement coil. The inverse-modeled results 
of the apparent relative magnetic permeability (ur) at 400 kHz, as given in Fig. 5, show 
that the corroded part is clearly distinguishable. Also, one profile line has been taken 
for the lift-off plots, of the same region under interest, shown as a brown vertical line in 
Fig. 5. Lift-off has been taken for four scans on different lift-off values, shifted vertically 
by 50 um/s from each other, as shown on an upper plot of Fig. 6 (a). The same four lift- 
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Fig. 4: The initial specimen. 


Fig. 5: Apparent relative mag- 
netic permeability of region 
no. 2 at 400 kHz. 


off profiles, but shifted vertically onto each other, are shown on the lower plot on the 
same Fig. 5. The quite exact matching of these four curves indicates the validity of the 
eddy current modeling and confirms the hypothesis that the corrosion depth can be 
characterized as change of the apparent lift-off of the material (carbon steel), in the 
total range of 55 um/s for the current example. 


2.3 Experiments with specimens no. 1 and 2 


Two specimens have been prepared for the experimental research (no. 1 and no. 2), as 
shown before and after painting in Fig. 7 and Fig. 8, with round- and strip-shaped cor- 
rosion sub-areas, correspondingly. It was noticed that, if using the described air-core 
planar coil, the resolution of determining of the shape and area of the corroded (rust) 
area is limited. Still, as seen from Fig. 9, for specimen no. 2, a more or less reasonable 
corrosion image is acquired. 
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Fig. 6: Four lift-off profiles, 
inversely calculated from the 
measurements at 400 kHz; 
absolute lift-off values on the 
upper plot and relative values 
on the lower plot. 
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Fig. 7: Specimen no. 1, before 
and after painting. 


3 Application to estimate the rust area 


To achieve better resolution in the estimation of the rust area shape and size, a mag- 
netic-core measurement coil (diameter of the core is 15mm and N = 50 turns) has 
been introduced (Fig. 10) to work together with the mentioned 3D scanner. As, differ- 
ently from the air-core approach, in real life it is hard to produce a precise eddy cur- 
rent model, due to several uncertainties (magnetic and geometrical properties of the 
coil, magnetic core, etc.) simply a complex impedance change has been taken into ac- 
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Fig. 10: Sensor coil with ferrite core. 


Fig. 8: Specimen no. 2, before 
and after painting. 


Fig. 9: Apparent magnetic 
permeability at 400 kHz of 
specimen no. 2. 
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count for rust measurement. Experimentally it has been found, for the current setup 
and for a working frequency of 400 kHz, the linear combination of the real and imagi- 
nary components of the impedance of the coil can be used. This combined impedance 
values are scanned in 2D at fixed lift-off and acquired as the 2D data array; further 
the specially developed PC Windows-based application is segmenting this 2D plot as 
“corroded” and “non-corroded” parts, by selectable threshold, showing the resultant 
image on the screen (Fig. 11 and Fig. 12). Additionally the user can mark the rectangular 
region of interest on the image, where the area of the corroded segment is calculated 
and displayed, in absolute (mm?) and relative (%) units. 


Threshold: 112 


MinLevel: 0 


MaxLevel: 220 


Fig. 11: Windows application 
for rust area measurement: 
specimen no. 1. 


4 Conclusion and future work 


The corrosion of metals, specifically the corrosion of the carbon steel, can be efficiently 
detected and measured by the eddy current approach with single-coil setup, also un- 
der a paint layer or similar coatings. 

The proposed work is similar to [19], but instead of a pulsed sensor signal the 
harmonic excitation has been used, with a single, selectable frequency at a time. Also, 
the usage of the air-core coil in combination with the Dodd—Deeds model enables to 
get rid of dependency on the sensor lift-off effect and also to estimate directly, without 
special calibration, the depth of the corrosion. 


132 — 0.Märtens etal. 


Threshold: 218 


MinLevel: 0 


MaxLevel: 254 


Threshold: 171 


MinLevel: 51 


MaxLevel 190 


Fig. 12: Windows application 
for rust area measurement: 
specimen no. 2, at lower (left 
image) and higher (right im- 
age) thresholds. 


The measure of the corrosion can be the variation of the apparent lift-off or change 
in the apparent magnetic permeability of the material under test. A resolution of 
1pm of the apparent lift-off variation can be achieved. A simple Windows applica- 
tion was implemented to measure the area of the corroded part(s). Further low-cost 
handhold instrumentation can be developed for such measurements. One future re- 
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search challenge is to develop the improved absolute models for eddy current-based 
corrosion measurement. Also, fusion of larger planar air-core coils with smaller fer- 
rite core probes into the same measurement procedure can combine the absolute 
(model-based) accuracy of the first and good spatial resolution of the second. 


Bibliography 


[1] 


[2] 


[3] 


[4] 


[5] 


[6] 


[7] 


[8] 


[9] 


[10] 


[11] 


[12] 


[13] 


[14] 


J. Garcia-Martin, J. Gömez-Gil, and E. Vazquez-Sanchez, “Non-destructive techniques based on 
eddy current testing”, Sensors, vol. 11, no. 3, pp. 2525-2565, 2011. 

C. Dodd and W. Deeds, “Analytical solutions to eddy-current probe-coil problems”, Journal of 
Applied Physics, vol. 39, no. 6, pp. 2829-2838, 1968. 

C. V. Dodd and W. E. Deed, “Calculation of Magnetic Fields from Time-Varying Currents in the 
Presence of Conductors”, NASA STI/Recon Technical Report N, vol. 76, Jul. 1975. 

S. K. Burke and M. E. Ibrahim, “Mutual impedance of air-cored coils above a conducting plate”, 
Journal of Physics D, Applied Physics, vol. 37, no. 13, pp. 1857-1868, 2004. 

C. Dodd and W. Deeds, “Absolute Eddy-Current Measurement of Electrical Conductivity”, Oak 
Ridge National Lab., TN (USA), Tech. Rep., 1981. 

O. Martens, R. Land, R. Gordon, M. Min, M. Rist, and A. Pokatilov, “Precise eddy current 
measurements: Improving accuracy of determining of the electrical conductivity of metal 
plates”, In: Lecture Notes on Impedance Spectroscopy: Measurement, Modeling and 
Applications, vol. 4, 2014, pp. 109-115. 

M. Parker, A. Pokatilov, T. Kubarsepp, O. Martens, and A. Kolyshkin, “Investigation of planar 
coil for eddy current conductivity measurements in wide frequency range”, In: 19th IMEKO TC4 
Symposium - Measurements of Electrical Quantities 2013 and 17th International Workshop on 
ADC and DAC Modelling and Testing, 2013, pp. 55-59. 

A. Kolyshkin, T. Kubarsepp, O. Martens, and A. Pokatilov, “Determination of electrical 
conductivity of metal plates using planar spiral coils”, In: Proceedings of the 13th WSEAS 
International Conference on Systems - Held as part of the 13th WSEAS CSCC Multiconference, 
2009, pp. 545-550. 

T. P. Theodoulidis and M. Kotouzas, “Eddy current testing simulation on a personal computer”, 
In: 15th World Conference on Nondestructive Testing, Roma (Italy), 15-21 October 2000, 
Retrieved 26 May 2016, 2002. [Online]. Available: http://www.ndt.net/article/wcndt00/ 
papers/idn228/idn228.htm. 

T. Theodoulidis and E. Kriezis, “Series expansions in eddy current nondestructive evaluation 
models”, Journal of Materials Processing Technology, vol. 161, no. 1-2, pp. 343-347, 2005. 

M. Fan, B. Cao, P. Yang, W. Li, and G. Tian, “Elimination of liftoff effect using a model-based 
method for eddy current characterization of a plate”, NDT and E International, vol. 74, 

pp. 66-71, 2015. 

A. Pokatilov, M. Parker, T. Kübarsepp, O. Martens, and A. Kolyshkin, “Grid-based computational 
algorithm for accurate ac conductivity measurements”, In: 16th International Congress of 
Metrology, EDP Sciences, 2013, p. 11007. 

S. Norton and J. Bowler, “Theory of eddy current inversion”, Journal of Applied Physics, vol. 73, 
no. 2, pp. 501-512, 1993. 

J. H. Rose, E. Uzal, and J. C. Moulder, “Magnetic permeability and eddy-current measurements”, 
In: Review of Progress in Quantitative Nondestructive Evaluation, Springer, 1995, pp. 315-322. 


134 — 0.Märtens etal. 


[15] Y. Nonaka, “A double coil method for simultaneously measuring the resistivity, permeability, 
and thickness of amoving metal sheet”, IEEE Transactions on Instrumentation and 
Measurement, vol. 45, no. 2, pp. 478-482, Apr 1996. 

[16] X. Ma, A. Peyton, and Y. Zhao, “Eddy current measurements of electrical conductivity and 
magnetic permeability of porous metals”, NDT and E International, vol. 39, no. 7, pp. 562-568, 
2006. 

[17] N. Bowler, “Frequency-dependence of relative permeability in steel”, In: AIP Conference 
Proceedings, vol. 820 Il, 2006, pp. 1269-1276. 

[18] Y. Gotoh, H. Hirano, M. Nakano, K. Fujiwara, and N. Takahashi, “Electromagnetic nondestructive 
testing of rust region in steel”, IEEE Transactions on Magnetics, vol. 41, no. 10, pp. 3616-3618, 
2005. 

[19] Y. He, G. Tian, H. Zhang, M. Alamin, A. Simm, and P. Jackson, “Steel corrosion characterization 
using pulsed eddy current systems”, IEEE Sensors Journal, vol. 12, no. 6, pp. 2113-2120, 2012. 


Ingo Tobehn-Steinhäuser, Susanne Kersten, 
Jaqueline Stauffenberg, Heike Wünscher, Christian Zeitnitz, and 
Thomas Ortlepp 


First measurements with a radiation hard 
relative humidity sensor for the ATLAS 
Experiment 


Abstract: A radiation hard humidity sensor, based on an inter-digital structure is de- 
veloped for the upgrade of the Large Hadron Collider (LHC) to high-luminosity LHC 
(HL-LHC). In this context, the interior of the detector of the ATLAS Experiment will 
be replaced. To guarantee a failure-free operation over the whole time, it has to be 
proven that the humidity in the detector region is negligible. To achieve the required 
very small size of the humidity sensor, double-sided wafer processing is used; fur- 
thermore the inter-digital structures consist of 3D fingers to decrease the necessary 
area. The otherwise customary but not radiation hard polymer as a sensitive layer is 
replaced by an SiO, matrix. In order to overcome the great distance between sensor 
and data processing, very compact signal preprocessing is also necessary. 


Keywords: capacitive humidity-sensor, radiation hard, low residual moisture 


1 Introduction 


Currently the ATLAS collaboration develops a concept for the new interior detector 
which is completely made of silicon semi-conductor detectors [1]. ATLAS was origi- 
nally an acronym for “a toroidal LHC apparatus”, but is now only used as a proper 
name. The innermost area is formed of multiple layers of pixel detectors, further out 
follow strip detectors. The coolant of -40°C requires an extremely low relative humid- 
ity (RH) of the atmosphere in order to avoid condensation on the electronics in all 
circumstances. Therefore, a permanent monitoring of the RH and the dew point by 
determining moisture sensors is essential for the safe operation of the detector. (The 
dew point of the environment should be < -40°C). In addition to the measurement of 
very low amounts of residual moisture at the low temperatures, the radiation hard- 
ness of humidity sensors is a further requirement. Assuming that the ATLAS detector 
will operate on the high-luminosity large hadron collider (HL-LHC) for ten years, one 
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would expect an integrated luminosity of 3000 fb’. Commercially available moisture 
sensors use polymers as a wet-sensitive layer and are unsuitable for use because of 
the lack of radiation hardness. Sensors used so far always had to be replaced after a 
relatively short time because they changed their characteristics [2]. 


2 Experimental setup 


2.1 The ATLAS Experiment 


For the upgrade of the LHC to the HL-LHC, the ATLAS experiment will replace the en- 
tire internal detector to handle both the high radiation load and the increased data 
rates. This is currently planned for the year 2022, but it is also possible that parts of 
the internal detector have to be exchanged in 2018. The ATLAS Collaboration is cur- 
rently working on a concept for this new internal detector, which will consist entirely 
of silicon semi-conductor detectors. The innermost region is formed from a plurality 
of layers of pixel detectors, and strip detectors follow further outwards Fig. 2. Though 
to a different extent, both detectors have an extremely high power density of the read- 
out electronics, which requires a correspondingly effective cooling (cooling power of 
240 kW). 

From Fig. 1, the requirement non-ionizing energy loss (NIEL) value is 4 x 10” par- 
ticles/cm?. The NIEL value is a theoretical quantity, which makes it possible to easily 
compare different effects, which play a particular role in the damage of silicon detec- 
tors. 
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Fig. 1: The calculated non-ionizing energy loss. 
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Fig. 2: The ATLAS Experiment at CERN. 


An evaporative temperature of the coolant of -40°C requires an extremely low RH of 
the atmosphere in order to avoid condensation on the electronics in all circumstances. 
Therefore, a permanent monitoring of the RH or the determination of the dew point 
by moisture sensors is an indispensable prerequisite for the safe operation of the de- 
tector (the ambient dew point should be < -40°C). In addition to measuring the very 
low amount of residual moisture at low temperatures, the radiation hardness of the 
moisture sensors is another requirement. Based on the fact that the ATLAS detector 
will be operating for 10 years on the HL-LHC, an integrated luminosity of 3000 fb * is 
expected. The atmosphere will usually be steamed off nitrogen or highly filtered clean 
air without dust or oil residue, so that no contamination of the sensor by deposits is 
to be feared. In order to influence the radiation length as little as possible, the sensor 
should introduce little additional material into the detector volume or preferably con- 
sist of materials with a low mass number. For the same reason, the number of leads 
for supply and readout should be as low as possible. Due to the extreme packing den- 
sity of the detector, the moisture sensor should also be as small as possible. Since 
ultimately the dew point is the critical variable, a temperature measurement should 
be integrated. The signal strengths must be sufficient to allow for transmission over 
at least 15-m cables. Depending on the type of the measured variable, local data pro- 
cessing is to be considered. The entire detector is located in a magnetic field of 2 T. 

The whole arrangement is designed for a duration of ten years. After the installa- 
tion, there is no access capability, which consequently places the highest demands on 
the reliability of all installed components. Maximum stability of all calibration factors 
is necessary, unless the sensor system has a self-calibration. 
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2.2 Design of the radiation hard humidity sensor 


2.2.1 Layout 


In order to achieve the required small size, silicon-on-insulator (SOI) material is used 
in order to be able to perform double-sided processing. This allows to use both the 
front and the rear of the chip for components. On the device side of the wafer with a 
thickness of 10 um, the 3D inter-digital structure is produced by deep reactive ion etch- 
ing (DRIE). The use of a 3D structure results in two advantages for the detector. First, 
the effective detector area is thereby increased with the same base area. Second, it 
is advantageous to use highly doped silicon instead of metal for the conductive struc- 
tures. This is in accordance with the requirement to use materials with a small number 
of nuclear charges. This reduces the activation of the material by the irradiation and 
the radiation length is influenced as little as possible by the sensor. 

Three different variants have been implemented with regard to the geometric lay- 
out. The combinations listed in Fig. 3 are used for the finger widths and distances be- 
tween the fingers. 


gap width [um] gap depth [pm] finger width [pm] 


1.1 10 3.5 
1.5 10 2.5 
3.0 10 4.0 


The smaller the distance between the fingers, the higher the total capacitance of 
the sensor. On the other hand, small finger distances make it difficult to insert the 
moisture-sensitive material between the fingers. Since this is an RH sensor, the tem- 
perature must be measured simultaneously to determine the dew point temperature. 
This is done according to the formula 


1 
Tr = f” (109.8 + T) - 109.8, (1) 


where Ty is the dew temperature point and fp is the RH. 

Like all other components, the diode, used to measure the temperature, is placed 
on the handle side of the wafer. In the first design, a micro-controller is used, which 
has very small geometric dimensions [3]. This concept is followed because the sensor 
must first be as small as possible, but on the other hand the output signals should be 
given via a 16-m long line (Fig. 4). In the case of the inter-digital structure, a square 
pulse is applied to the structure via a fixed resistor. On the basis of the capacitance, 
which has been adjusted due to the humidity, there is a decay time due to the RC con- 
stant. 

This is analyzed by measuring the voltage after a fixed time t. In order to provide 
the information of this analog signal via the signal line, it is digitized with an ADC inte- 
grated on the micro-controller. The micro-controller used is commercial standardware 
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Fig. 3: Layout ofthe sensor. 


tion of the wafer. 


— 139 


Fig. 4: Double-sided prepara- 


and therefore does not meet the increased requirements at CERN, in particular the 
very high radiation load. Therefore, this part is only a conceptual design and should 
be transferred from a sub-center of CERN to a layout with radiation-proof electronics. 
The large solder pads in the layout are only for programming the micro-controller and 
do not belong to the definitive concept of the sensor. For the online calibration of the 
system, the installation option of a capacitor is also provided as a reference on the 


circuit board. 
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All previous RH sensors used in CERN failed due to their insufficient radiation 
resistance. The moisture-sensitive layer consists of polymer chains on all commer- 
cially available sensors. These are split at high beam exposures. In this way, the sensor 
changes its response behavior and thus its characteristics. In the concept pursued at 
CiS, SOL-GELs are used. After processing, these form an SiO, matrix in which some 
water binding substances are additionally incorporated. Furthermore, in the new de- 
velopment of the SOL-GEL used, no organic solvents were used but only water [4]. 

In dip coating, a sample is placed in a liquid medium to coat it. The sample is 
dipped into the substance, then remains there for a certain time t and is then pulled 
out at a constant velocity v,. After this, the sample has to dry for a certain time and 
may have to be dried again in an oven at higher temperatures. For this purpose, the 
layer thickness can be calculated using the Landau-Levich, i.e., 


d = 0.8, 2. (2) 
\pz 


The layer thickness d is determined by this formula from the viscosity n, the draw- 
ing speed vo, the density p of the liquid and the gravitational constant g [5]. 

For applying the SOL-GEL to the sensor layer, the spin coating method was un- 
suitable since the 3D sensor structures were not filled (Fig. 5). The best results were 
obtained by the dip coating method, whereby the structure to be coated is slowly 
drawn out of the solution and subsequently annealed. It was also found to be advan- 
tageous during the immersion phase to sonicate the solution with ultrasound in order 
to achieve a true homogeneous coating. 


Fig. 5: 3D structure. 


2.3 Measurements 


2.3.1 Uncoated sensors 


The sensors were initially measured dry and with de-ionized water using an LCR meter 
(Agilent 4294A). The results are shown in Fig. 6. 

Significant differences of the sensor signal at room temperature in the dry state 
and with the applied drop of de-ionized water can be seen. The water quantity 
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Fig. 6: (a) Uncoated Sensor. (b) Coated Sensor. 


was 5uL. The difference corresponds to about 350 pF at 100 Hz and is thus easily de- 
tectable. In this way, the general functionality of the first built-in sensors was checked. 
After this function test, the sensors were provided with the moisture-sensitive layer 
for further measurements. 


2.3.2 Coating with dip coating 


The temperature was kept constant at 25°C during the measurements and the RH var- 
ied between 10% and 60%. First, the capacitance was measured as a function of 
the frequency. The measured frequency range was between 100 Hz and 1 MHz. Fig. 7 
shows the difference between the two moisture values. 

Based on these measurements, a value of 5 pF/ % RH at 100 Hz can be calculated. 
The capacity at 100 Hz and 10% RH is 317 pF and with a humidity of 60 % 575 pF is 


C1 “T 
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0 50.000 100.000 
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Fig. 7: Equivalent circuit. 
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achieved. In addition, the impedance and the phase were measured in a frequency 
range of 25 Hz to 1 MHz. Again in this case, the difference between the different mois- 
ture values can be clearly seen. At 10 % RH, the impedance is 6.89 MOhm and with a 
humidity of 60 % 3.71MOhm. The course of the measurement curves corresponds to 
that of RC elements. 


2.3.3 Irradiation of samples 


To investigate whether the new layer used will withstand future radiation exposure 
without changing its required properties, several samples were irradiated at the 
27-MeV proton linear accelerator in Birmingham [6]. Prior to irradiation and after 
irradiation, FTIR spectra were recorded to investigate whether the chemical bonding 
structures detectably altered. In order to obtain a better indication of the stability of 
the coating, different samples were irradiated at different doses. These had the values 
0.77 MGy, 1.11 MGy, 4.84 MGy and 6.12 MGy. Even at the highest dose, no significant 
differences were found between irradiated and unirradiated specimens. In a next 
step, a mounted unit under irradiation is to be subjected to a functional test. 


3 Summary 


The experimental results so far suggest that the concept underlying the sensor works 
reasonably. The vias proved to be critical; the technology still requires a modification 
to make it less susceptible. Moreover, it was found that the method of dip coating for 
applying the sensitive layer appears to be the only practicable method. For this pur- 
pose, the process parameters such as pulling speed and tempering times had to be 
found by means of corresponding test series. Further, an additional application of ul- 
trasound during the coating process was useful to achieve a more homogeneous coat- 
ing image. Experiments with spin coating showed unsatisfactory results because the 
3D structures are not filled sufficiently. 
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Frank Wendler, Paul Aurosmit, and Olfa Kanoun 
CNT-enhanced carbon electrode for nitrite 
detection and water analysis 


Abstract: In this contribution we present a novel electrode based on a bare pencil lead 
electrodes (PLEs) coated by carbon nanotubes. We investigate the sensitivity of the 
novel electrode during electrochemical analysis of aqueous solutions. For this pur- 
lose, a phosphate buffer having different nitrite concentrations has been subjected to 
meaurements with enhanced and bare PLEs over a wide potential range. The proposed 
coating increases the sensitivity to nitrite by increasing the electrode surface and its 
reactivity and shows therefore a significantly observable and sensitive signal. 


Keywords: square wave voltammetry, nitrite, pencil electrode, carbon nanotubes 


1 Introduction 


Nitrite ions (NO3) can be found in food, water and soils [1]. In the food industry, nitrite 
is used in the curing of meat by preventing bacterial growth [2]. It is found as a natural 
component in vegetables and water due to chloramination [3]. High concentrations of 
nitrite ions are toxic to humans and animals. Especially during cooking, excess ni- 
trite in meat can react with degradation products of amino acids to form carcinogenic 
nitrosamines [1, 4-6]. It also oxidizes hemoglobin to methemoglobin, which causes 
health problems like methemoglobinemia [7]. 

To detect harmful concentrations of nitrite in the food chain of human and ani- 
mals, a reliable and cost-effective way of determining the concentration of nitrite in an 
aqueous solution is needed. The method of square wave voltammetry (SWV) provides 
excellent results in the quantitative analysis of chemical compounds in low concen- 
trations [8]. The selective nature of this method is very useful to detect chemicals that 
are involved in electrochemical redox reactions. The sensitivity and the selectivity of 
the method can be further enhanced by specialized electrodes. 

This enhancement is usually performed by coating a bulk electrode material 
within an electrocatalytic layer. Often used bulk materials for the electrodes are 
glassy carbon and pyrolytic carbon. These are expensive and are usually cleaned 
and re-used in several experiments. To evaluate a large amount of coatings in the 
optimization process a more cost-effective base electrode is needed. 

In this work pencil lead electrodes (PLEs) are proposed as cheap and accessi- 
ble base material for electrode prototyping. In an additional stage the electrodes are 
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coated with multiwalled nanotubes (MWNTs) to enhance the sensitivity to the concen- 
tration of nitrite ions as they have a strong electrocatalytic property [9]. The aim of the 
paper is to invest the feasibility of this type of electrodes and to evaluate its sensitivity 
for measurements. 


2 Experimental 


2.1 Solution preparation 


For the determination and graphical representation of the reaction kinetics, SWV was 
used. 

SWV is known for well-defined current peaks at characteristic redox potentials 
[8]. It is also very sensitive to the electrode kinetics [8]. The redox reaction that takes 
place at the electrode is 


NO, — NO, +e. (1) 


A 0.65-M (pH = 7) phosphate buffer solution (NaH,PO,/Na,HPO,) from Grüssing 
GmbH and sodium nitrite of 99% purity was used to prepare the investigated solu- 
tions. The buffer solution was diluted to 0.325 M; 100 mL of the diluted buffer solution 
was contaminated with 0.7 mg sodium nitrite and 100 mL of the buffer solution was left 
uncontaminated. Both base solutions were used to create samples with mixing ratios 
of 0%, 25%, 50%, 75 % and 100 %. In this 100 mL of base solution the nitrite concen- 
tration was 100 uM. So in this investigation the different samples of 100 %, 75 %, 50 %, 
25 % and 0 % correspond to concentrations of 100 uM, 75 uM, 50 uM, 25 uM and 0 uM. 


2.2 Preparation of electrodes 


The pencil lead consists of graphite, clay and wax has a resistance of approximately 6 
Ohms. The whole pencil was heated to reduce the resistance and evaporate the excess 
binding material so that the properties of graphite are predominant. The pencil lead 
is covered by an insulating cover so that only the tip of the pencil is exposed to the 
solution. 

For the enhancement, multiwalled carbon nanotubes (MWCNTs) of 2% w/w in 
2-propanol (stabilized) were used from Future Carbon GmbH. The pencil lead was 
cleaned and then the tip was dipped in the prefabricated MWCNT solution. 

The pencil lead was dried in a heater for 5 hours to evaporate the propanol and 
fixate the MWCNTs on the surface of the pencil lead. 


CNT-enhanced carbon electrode for nitrite detection and water analysis —— 145 


2.3 Instrumentation and experimental setup 


A three-electrode cell was used for the electrochemical experiments. The measure- 
ments were conducted using the HEKA 310 potentiostat. The pencil lead was used as 
the working electrode. Ag/AgCl 0.1 M KCl from CH Instruments, Inc. was used as the 
reference electrode. The counter electrode consisted of the platinum wire. The mea- 
surement setup is described in Fig. 1. 


Reference N s 


electrode Potentiostat 


Auxiliary 
electrode 


Fig. 1: Experimental setup [10]. 


2.4 Instrument settings 


The start potential was set at —0.5 V and the end potential at 1.2 V. During the measure- 
ments in the negative potential range the samples were able to reach a steady state. 
The potential was increased in steps by 10 mV, and the square wave frequency was set 
to 1 Hz. The scan rate was set at 10 mV/s. The current was sampled from 98 %-100 % 
of the period so that the double layer capacitive effects can be avoided and only the 
Faradaic currents are represented in the data. 


3 Results and discussion 


The electrode was immersed in all five samples solutions and the SWV was conducted. 
The bare PLE in Fig. 2 shows a small current peak in the region of 850 mV. This peak 
is only visible at high concentrations of 100 % and 75% of the base nitrite solution. 
We conclude that the bare PLE has a limited use in the detection of nitrite, but can be 
used as base electrode to observe the effects of enhancing coatings. 
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Fig. 2: SWV of pencil lead at different concentrations of nitrite. 
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Fig. 3: SVW of CNT-coated PLE at different concentrations of nitrite. 


The PLE was coated with the MWNT film and immersed in the same five sample solu- 
tions. The coated electrodes show an increased current at all concentrations as seen 
in Fig. 3 due to the increase in surface area of the nanomaterial coating. The current 
peak at 0.85 V for the nitrite is enhanced and also observable at lower concentrations 
of nitrite, as shown in Fig. 4. This increases the sensitivity and the detection limit for 
nitrites of this electrode type. Furthermore, a second peak can be observed at 0.275 V. 
This peak has a constant level at all concentrations of nitrite. All the five samples al- 
ways contain the same amount of phosphate buffer and de-ionized water irrespective 
of the nitrite concentration. So we conclude the peak at 0.275 V is related to the con- 
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Fig. 4: SWV of CNT-PLE showing the peak currents at different concentrations of nitrite. 
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centration of the phosphate buffer. Since the SWV waveform for the bare PLE does not 
show any prominent peak at any potential, we assume that CNTs have a catalytic effect 


on the phosphate in the buffer solution. 


The direct comparison between the voltammograms of the bare PLE and CNT- 
enhanced PLE (CNT-PLE) for base solutions with concentrations of 0 % and 100 % of 
nitrite can be seen in Fig. 5 and Fig. 6, respectively. The peak related to the base solu- 
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Fig. 5: SWV of PLE and CNT-PLE in 0 % nitrite solution. 
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Fig. 6: SWV of PLE and CNT-PLE in 100 % nitrite solution. 


tion at 0.25 V is observable in both voltamogramms for the enhanced electrode. The 
nitrite peak at 0.85 V can be observed with both electrodes in solutions containing 
nitrite. 


4 Conclusion 


In this chapter, a novel CNT-coated PLE is proposed for nitrite measurements by SWV. 
The results show that the PLE can be used as an excellent base electrode for modifi- 
cation and evaluation. The modification by coating with CNTs show good results for 
concentrations of nitrite from 100 pM to 25 uM. 

The CNT-coated electrode can be made sensitive to other chemicals, like phos- 
phate, to adapt to other chemicals and fields of application. The electrode material is 
cheap and can be enhanced with minimal effort. Therefore PLEs support a screening 
of large quantities of coating materials. For future work the miniaturization and opti- 
mization will be the main aim along with checking the sensitivity to other reagents, 
so that cheap and reliable electrochemical sensors can be created. 
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Part IV: Bioimpedance measurements and methods 


Rauno Gordon, Mart Min, and Raul Land 
Energy-efficient pulses for deep brain 
stimulation 


Abstract: Energy-efficient excitation of neural tissue is investigated in this paper 
through the generation of optimal shape pulses. It is essential for neural stimulation 
to have a large percentage of the energy in the most useful frequency range. The useful 
range is above the frequency of the tissue—electrode interface capacitance and below 
the frequency of other parasitic capacitances (small capacitances of cables, connec- 
tors). One powerful observation to determine that useful range is that it centers on the 
minimum phase angle of the impedance. In this study we have generated stimulation 
signals in the form of short chirp pulses that have 95% of energy in one decade of 
frequency range. 


1 Introduction 


Electrical stimulation of neural tissue (deep brain stimulation) is an established 
method for the treatment of neural diseases. The American Food and Drug Admin- 
istration (FDA) approved deep brain stimulation as a treatment for essential tremor 
in 1997, for Parkinson’s disease in 2002 and for dystonia in 2003 [1]. For treatment, 
stimulation leads are placed in the brain according to the type of symptoms to be 
addressed. 

Practical application of the method needs extremely high energy efficiency. The 
power of the pulses has to be exactly controlled because the tiny electrodes with high 
electrical impedance can degrade themselves but they can also heat up the nerve tis- 
sue, damage it and modify the neural functioning. Nowadays mostly fully implanted 
devices are used for neural tissue stimulation (battery-powered heart and brain pace- 
makers). Their functions include continuous monitoring, signal processing, decision 
making and stimulation when required. They should be able to operate for years. 
Therefore simple and extremely low-power electronics is required for the generation 
of stimulation pulses with optimized waveforms. 

The spectra of the stimulation pulses need to be matched with the spectral prop- 
erties of the complex electrode-tissue impedance Z(jw) for the desired neural tissue 
stimulation effect. The energy of the stimulation signal should desirably be concen- 
trated into a spectrum band between indicative frequencies f} and f,. We should avoid 
generating energy with frequencies under the lower bound f, where the stimulation 
signal can poorly penetrate into the tissue due to double-layer impedance on the elec- 
trode interfaces. The frequency above f, should again be avoided because there are 
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parasitic losses due to small capacitances and inductances in the whole stimulation 
pathway. 


2 Objectives 


Sequences of pulses with the repetition rates from 100 Hz to 200 Hz and durations be- 
tween 100 us to 500 us are typically used for stimulation [2]. Our aim is to find specific 
waveforms for the stimulating pulses that have most of the energy in the most sensi- 
tive frequency range for the neural tissue. We try to focus as much energy as possible 
to anarrow energy spectrum band. The center of the frequency band would be the fre- 
quency where the impedance has its maximum imaginary part — also called the peak 
resistance frequency f, [3, 4]. Our metric for measuring the focussed energy is the per- 
centage of the energy that falls into one decade of the spectrum (Fig. 1). For the most 
promising excitation signal we also look at the spectrum of repeating pulses with var- 
ious repetition rates to see if a practical usage scenario can affect the spectrum and 
therefore the energy efficiency of the stimulation. 


Oe ee 


| | pulse spectrum, V(w) 


Z(jw) - magnitude 


Fig. 1: Simplified magnitude and 
phase of tissue impedance (blue solid 
lines) with desired pulse spectrum 
(red dashed line). Frequencies f} and 
f are bounds between which we want 
to focus the energy of the pulse; f, is 
the frequency where impedance has 
its minimum phase angle (peak). 
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3 Methods 


Excitation pulses for nerve tissue must be exactly balanced - for a single excitation 
pulse there should be a positive and a corresponding negative part of the signal where 
electric charge given into the electrode is later balanced. If unbalanced charge is left 
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into the tissue near the electrode, tissue damage can occur over time and electrode 
degradation can follow too, which would limit the longer-term useability of the de- 
vice [5]. 

We have previously studied square pulses (also balanced) with and without a gap 
between the excitation and the balancing pulse [6]. Square pulses have the advantage 
of being easy to generate with simple low-power electronics. We looked into square 
pulses with various ratios between the amplitude and longevity of the excitation pulse 
and balancing pulse. We found that square pulses can be generated with up to 85 % 
of energy in the desired band within one decade. A good and typical example is given 
in Fig. 2. A gap is often added in practice between the excitation and balancing part of 
the pulse [7]. Based on our results it did not improve the energy efficiency of the pulse 
but degraded it, depending on the length of the gap [6]. 


A1 = 0.707 


Fig. 2: Above: balanced square excita- 

tion pulse with an amplitude ratio of 

0.707. Below: spectrum of the signal 

calculated with discrete FFT showing 

10 10 10° 10° 10 10° the bounds with gray lines where useful 
frequency [Hz] signal energy is calculated. 


The stimulating pulses with the predetermined repetition rate and duration are syn- 
thesized mathematically for obtaining the focused and controllable energy spectrum. 
Several pulse signals with a frequency that rises following square function are exper- 
imented with. 

Pulses of the first type are produced with frequency rising following a square func- 
tion and when the sinus signal reaches the quarter-period, the frequency is lowered 
following the same square function. This however produces a signal that is not sym- 
metrical with respect to the central line where the frequency rise was reversed. The 
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Fig. 3: Above: part of balanced signal 
magnitude (solid line) and rising fre- 
quency (dashed line). Below: spectrum 
ofthe signal calculated with discrete 
FFT shown by gray lines where the sig- 
nal energy is calculated. 


frequency function and an example ofthe pulse can be seen in Fig. 3. Please note that 
the pulse includes a balancing, negative part as well that is identical but not shownin 
the figure. After generating the balanced parts of the pulses, the frequency spectrum 
was calculated with discrete FFT in MATLAB. Energy within one decade of frequency 
was calculated for comparison purposes. The selection ofthe actual desired frequency 
band depends on the application but here the one-decade band is chosen based on 
typical tissue conditions. 

Pulses of the second type use the rising frequency as well, but the half-period of 
the signal is used. The same signal in reverse sequence is then added for the balance 
on the negative side. The rising frequency function and the best sample pulse with 
both positive and negative parts is shown in Fig. 4. The spectrum and energy in one 
decade were calculated after that. 
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Fig. 4: Above: balanced signal magni- 

tude (solid line) and rising frequency 

(dashed line). Below: spectrum of the 

f | signal calculated with discrete FFT 

10° 10° 10° 10° i 40° shown by gray lines where the signal 
frequency [Hz] energy is calculated. 


4 Results 


The spectrum and energy efficiency of the first type of balanced signal are given in 
Fig. 3. We selected one decade in frequency to be the band where we desire to concen- 
trate the energy. It can be seen that 95 % of the energy resides in the desired region. 

The spectrum and energy efficiency of the second type of the pulse are given in 
Fig. 4. The frequency spectrum exhibits a smooth and narrow curve that also concen- 
trates around 95 % of energy into the desired band. 

The second type of pulse - half-cycle chirp pulse with added identical but re- 
versed balancing part — shows a remarkable and smooth drop in energy at higher 
frequencies (Fig. 4). This makes it especially suitable for stimulation when parasitic 
losses at higher frequencies are limiting. We also investigated the spectrum of this 
second type of pulse in situations where the stimulation needs to be repeated. The 
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Fig. 5: Above: a signal with two bal- 
anced half-cycle chirps with a gap 

in between. Below: spectrum of the 
signal calculated with discrete FFT 
shown by gray lines where the signal 
energy is calculated. The signal re- 
tains 96 % of the energy in the desired 
frequency [Hz] one-decade band. 


results (Figs. 5 and 6) show that, even if the half-cycle balanced chirp is used for re- 
peated stimulation, the energy efficiency is still retained above 90 % in the desired 
one-decade band. 


5 Conclusion 


It has been found that short balanced chirp pulses (chirplets or titlets) with a dura- 
tion of only a quarter or half of the full cycle concentrate about 95 % of the generated 
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UM. f [He] 


Fig. 6: Above: a signal with four bal- 
anced half-cycle chirps with a gap in 
between. Below: spectrum of the signal 
calculated with discrete FFT shown by 
gray lines where the signal energy is 
calculated. The signal retains 93 % of 
the energy in the desired one-decade 
frequency [Hz] band. 


energy into the most desirable frequency range from f; to f,. The optimal frequency fp 
[3, 4] is the peak value of the phase spectrum of Z(jw). The chirp-based stimulation 
pulses with optimized waveforms enable to concentrate the stimulation energy pre- 
dominantly into the most sensitive frequency band. Short half-cycle chirps with an 
identical but reversed balancing part show remarkable concentration to the limited 
frequency band and are the preferred waveforms. Even if used for repeated stimula- 
tion, the spectrum stays limited to the desired band with above 90 % efficiency. The 
true effectiveness of the balanced chirp pulses for deep brain stimulation remains to 
be tested in vivo. 
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Jaan Ojarand, Raul Land, Mart Min, and Marek Rist 
How many frequencies to use in electrical 
bioimpedance measurements 


Abstract: The required frequency range of electrical bioimpedance (EBI) measure- 
ments depends on the dielectric properties of the object and covers about three to 
four decades typically within the kHz to MHz range. The accuracy of fitting of the 
equivalent circuit-based impedance model with the measurement results depends on 
the number of frequency components used. In general, the use of more frequencies 
improves the fitting accuracy but complicates the measurements. However, the fitting 
accuracy also depends on the signal-to-noise ratio (SNR) of measurement signals. If 
covering the frequency range with a simultaneous multifrequency excitation, the root 
mean square values of its frequency components are decreasing when their number 
increases. It follows that the number of used frequency components has two contra- 
dictory effects on the accuracy of estimation of parameters of the equivalent circuit 
model. This chapter deals with the analysis of typical EBI model-fitting situations 
and shows how the number of frequency components influences the fitting error. Our 
research shows that in the case of multifrequency excitation the accuracy of fitting 
of the typical equivalent circuit-based EBI model may not increase when the num- 
ber of frequencies is increasing. The results described here could be useful in other 
implementations of impedance spectroscopy. 


Keywords: electrical bioimpedance, impedance measurement, impedance spec- 
troscopy, equivalent circuit model, estimation of parameters 


1 Introduction 


Electrical bioimpedance (EBI) spectra are widely used to characterize the structure of 
tissues and cell cultures [1]. For the characterization of biological matter, wide band 
impedance spectroscopy is required that typically covers three to four frequency 
decades. 
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An important efficiency criterion of the of impedance measurements is the signal- 
to-noise ratio (SNR) which is proportional to the power of measurement signal. In EBI 
measurements, unfortunately, there is no possibility to improve the SNR by increas- 
ing the overall amplitude of the excitation signal since it is limited due to different 
reasons, such as non-linearity and safety requirements, as well as power consump- 
tion and some other technical issues [2, 3]. 

In situations where the properties of the measurement objects are changing in 
time, e.g., a beating heart muscle or fast moving cells in a high-throughput micro- 
fluidic channel, the frequency range of interest must be covered within a short time 
frame. The use of multifrequency excitations containing several frequency compo- 
nents simultaneously as a multisine signal does [3] gives us a possibility to perform the 
measurement procedure in a much shorter time span than in the case of subsequent 
measurements at a single hopping or sweeping frequency. We consider such cases in 
this chapter. However, the disadvantage is that energy of the multifrequency excita- 
tion spreads among multiple frequency components, which decreases the power of 
each component in the case of a limited total amplitude of the signal. 

Another factor, characterizing the power limitation of individual signal compo- 
nents, is the crest factor (CF) of a complex excitation signal — the ratio of its peak and 
root mean square (RMS) values. In the current analysis, we use a multisine signal (a 
sum of several sine waves with different frequencies) as the multifrequency excitation. 
For a graspable comparison of results, some variables must be predetermined. Since 
the overall amplitude of signals is limited, it is reasonable to normalize the amplitude 
(peak value Apk) of the sum of k components of the multisine Max|s(t)| = Apk = 1. In 
the case of equal amplitudes of k multisine components, their individual RMS values 
are also equal and depend only on a number of components k and the CF [3], i.e., 


Apk = 1 
Vk x CF Vk x CF 


Srms(i) = (1) 

The units of signal S can be either V or A, depending on the choice to use voltage 
or current excitation. In Fig.1, the essentials of formation of the EBI described in [4] 
and a simplified equivalent circuit of EBI are illustrated. 

At lower frequencies, current passes mostly through extracellular fluid. At higher 
frequencies, when the reactance of the membrane capacitances becomes lower, it 
also passes through intracellular fluid. A more accurate electrical model also includes 
both, a constant phase element (CPE) acting at lower frequencies (a dispersion range) 
and a stray capacitance of the measurement circuitry. The required frequency range 
depends on dielectric properties and the structure of the sample under test (SUT) 
covering decades within the kHz to MHz band [1]. 

The accuracy of fitting of the frequency response curve with the equivalent cir- 
cuit composed on the bases of measurement results depends on the number of com- 
ponents in the equivalent circuit, the shape of the frequency response of impedance 
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Fig. 1: (a) Essentials of formation of the EBI. (b) The three-element equivalent circuit. 


and, finally, also on the number of frequency components used. The use of more fre- 
quencies improves the accuracy because the shapes of the curves of the spectra are 
followed better. However, the shapes of EBI spectra are relatively smooth since the 
elements of the equivalent circuit are not producing resonances. 

Furthermore, the power of noise and its spectral properties have a substantial neg- 
ative impact on the fitting accuracy. 

Not many results of the systematic study of the number of required frequencies are 
available. Analysis of the three-element equivalent circuit obtained with sequential 
single-frequency excitation is described in [5] and [6]. 


2 Methodology 


A method used for testing of the influence of the number of frequencies k on the ac- 
curacy of estimation of values of equivalent circuit elements is illustrated in Fig. 2. 
Computer simulation of all the parameters and the equivalent circuit at the predeter- 
mined known values of its model elements are used instead of a real SUT. In real ex- 
periments, there are more variables, e. g., parameters of electrodes, internal noises of 


Spectrum ot the multisine excitation F v(t) 
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adjustment elements 


Fig. 2: A simplified structure of calculations used for testing of the influence of the number of fre- 
quencies k on the accuracy of the estimation of values of equivalent circuit elements. 
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the SUT, stray impedances of the connection wires and measurement circuitry. More- 
over, an equivalent circuit model of real SUT is not accurate. The use of synthetic data 
fixes the number of variables and allows focusing on the dependence of the fitting 
accuracy on a number of excitation frequencies. 

It is assumed that the amplitude of the multisine excitation stays on thelevel of+1 
V, and the initial phases of frequency components are optimal for achieving the best 
CF. An optimization algorithm described in [3] is used. 

We have used five quasi-logarithmically distributed frequency sets in the current 
analysis. Three decades are covering the frequency range from 1 kHz; see Table 1. The 
number of frequency points per decade varies from 1 to 5. Amplitudes of all the sinu- 
soidal components of the excitation voltage are kept equal, but their value is V(i) = 
v2 x Vems(i), the value that decreases when increasing the number of components; 
see equation (1). 


Table 1: Frequency sets of multisines and corresponding crest factors. 


k List of frequencies (kHz) CF 
3 1,19,401 2.43 
6 1, 3,13, 47,167, 607 2.77 
9 1,3,7, 19, 41, 83, 179, 379,797 2.58 
12 1,3,5, 7,17, 29, 47, 79, 139, 241, 421, 727 2.41 
16 §1,3,5, 7,11, 17,23, 31, 47, 73, 109, 163, 241, 359, 523,787 2.22 


A simple three-element equivalent circuit model (Fig. 1) and a more complicated six- 
element equivalent circuit model (Fig. 3) containing a constant phase element (CPE) 
were both used. 

Magnitudes of each frequency component of the response current were obtained 
as a ratio of the magnitudes of the excitation voltage and impedance of an equivalent 
circuit, as shown in Fig. 3. The phase of each component of the response current was 
found similarly by subtracting the phase values. A waveform of the response current 


G Rint 
R ext 


(a) (b) 


Fig. 3: (a) Three-element and (b) six-element equivalent circuit models with values of C = 5nF, 
Rint = 1kQ and Rext = 560 Q. The additionally introduced elements have values Cp = 10 pF, n = 0.8 
and Q = 10°°5, 
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i(t), in Fig. 2, is asum of k sinusoidal components with magnitude and phase values 
determined by: (a) - known multisine excitation waveform; and (b) - primarily un- 
known values introduced by the impedance of an equivalent circuit. 

Uniformly distributed Gaussian noise adds to the response signal in the next step, 
and the spectrum of the summary signal is calculated using the Fourier transform (FT) 
with 4096 points of the waveforms. For a graspable comparison of the results, the RMS 
value of added noise remains constant when the number of a frequency component 
varies. The impedance spectrum as the ratio of spectra of excitation and response sig- 
nalsis calculated, and the fitting of noisy impedance curves with a frequency response 
of the equivalent electrical circuit is performed. Then the differences between the ob- 
tained and ideal values are calculated, and the data for histograms are collected. 

The fitting algorithm minimizes the mean square error (MSE) as a relative measure 
of the residuals between the modeled curve values and the actual observed values. Pa- 
rameters of all components of the equivalent circuit are fitted at once. Along with the 
residuals also the relative importance of each component of the equivalent circuit at 
all the frequencies is shown, which allows the user to judge the adequacy of the model 
and to adjust its parameters if necessary; see the example in Fig. 4. In the current ex- 
periments, however, a priori knowledge of the model’s topography and the range of 
parameters simplified the task. 

Finally, the errors of mean values and their standard deviations are compared for 
a different number of frequency components. 

The SNR of the response signal is kept at the level of 10 dB, and 500 readings of the 
response and excitation signal spectra are used for each fitting of model parameters. 
Then, 500 fitting results are chosen for calculating the standard deviation of the error 
level for each number of frequencies. 
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Fig. 4: (a) A snapshot of the magnitude spectrum fitting error and (b) relative significance of com- 
ponents in formation of the magnitude spectrum, in the case of a three-element equivalent circuit 
(described in Fig. 3 (a)) with multisine excitation with six components. The upper dotted line in (b) 
represents Rext, the dashed line represents Rint and the dash-dotted line represents C. 
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3 Results 


In Figs. 5 and 6, snapshots of noisy spectra used in experiments for the fitting of val- 
ues of three-element and six-element equivalent circuits are shown.The SNR of the 
response signal is kept at the level of 10 dB in all cases and histograms of the 500 fit- 
ting results have been collected. 

In Fig. 7, histograms are shown of the 500 fitting results collected in the case of 
the three-element equivalent circuit for different values of C and Rint using multisine 
excitation with three components. 


Mag. Z, Ohm Phase, degree 
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Frequency, Hz Frequency, Hz 


(a) (b) 
Fig. 5: (a) Snapshots of the noisy magnitude and (b) phase spectra, in the case of a three-element 
equivalent circuit (shown in Fig. 3 (a)) and multisine excitation with six frequency components. The 


dotted lines represent calculated ideal spectra and the circles denote instant values of spectra at 
used frequencies. 
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Fig. 6: (a) Snapshots of the noisy magnitude and (b) phase spectra, in the case of a six-element 
equivalent circuit (shown in Fig. 3 (b)) and multisine excitation with nine frequency components. The 
dotted lines represent calculated ideal spectra and the circles denote instant values of spectra at 
used frequencies. 
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Fig. 7: Histograms of 500 fitting results of the values of (a) three-element equivalent circuit elements 
C and (b) Rint given for the multisine excitation having three components. 
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(a) (b) 
Fig. 8: The standard deviation of the fitting error of circuit element values of three-element equiva- 


lent circuit elements in the case of (a) sequential single frequency excitation and (b) multifrequency 
excitation; circular dots represent C, rectangles represent Rint and circles represent Rext. 


The dependence of the standard deviation of the fitting error from the values of three- 
element and six-element equivalent circuit elements on the number of frequencies is 
illustrated correspondingly in Fig. 8 and Fig. 9. 


4 Conclusion 


Increasing the number of frequencies has a low impact on the accuracy of fitting re- 
sults in the case of using multifrequency excitation and simple three-element equiva- 
lent circuits. Even three frequencies could provide sufficiently close results. Reduced 
energy of individual signal components appearing when adding the frequencies re- 
strains the effect of using a larger number of frequencies. 
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Fig. 9: Standard deviation of the fitting error of the values of elements of six-element equivalent cir- 
cuit in the case of (a) multifrequency excitation and (b) a corresponding equivalent circuit. Values 
of C, Rint and Rext are the same as of the three-element circuit in Fig. 1, but the values of addition- 
ally introduced elements are Cp = 10 pF, n = 0.8, and Q= 10%, 


It is reasonable to keep the number of frequencies as low as possible at the ac- 
ceptable errors. Expectedly, the number of frequencies should be larger, if the noise 
distributes non-uniformly over the frequencies or some concentrated disturbances ap- 
pear at some frequency sub-ranges. Every specific case deserves special handling. 

The positive effect from the use of more frequencies is larger in the highly noisy 
environment if the single-frequency signals are applied sequentially for excitation, 
but then the measurement process becomes too long for the measurements of time- 
varying impedances. 
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Andres Kink, Marek Rist, Raul Land, Hip Köiv, and Mart Min 

A method of intra-ventricular bioimpedance 
spectroscopy to estimate the dynamic volume 
of right ventricle 


Abstract: Dynamic volume measurements of the heart are primarily done in the left 
ventricle, but as the right ventricle is a frequent place for pacing electrodes it could 
be easier to estimate cardiac output through the right ventricle. Our goal was to find 
the correlation between intra-cardiac bioimpedance deviations and the variations re- 
lated to changes of the right ventricle volume. The variations of right ventricle volume 
have been evaluated experimentally using a setup with an isolated pig heart. Mea- 
surements were made with an intra-ventricular catheter and using multifrequency 
impedance spectroscopy. 


Keywords: bioimpedance, cardiac output, right ventricle volume, estimation 


1 Introduction 


In Europe, 45 % of all deaths are a result of cardiovascular disease, and the costs that 
derive from heart failure exceed the costs of all other diseases together [1]. Therefore, 
the motivation to develop a method to improve diagnostics and prevent heart-related 
problems is considerable. The importance of dynamic volume changes and the func- 
tion of the right ventricle have increased in the last decade because they provide nec- 
essary information to assess cardiovascular performance [2]. Measuring blood volume 
changes in the right ventricle is essential to estimate total cardiac output and it is an 
important parameter to characterize the heart and its pump function in experimen- 
tal studies. Cardiac output is usually measured from the left ventricle [3] because its 
geometry is relatively simple - a sphere - while the right ventricle is like an unequal 
layer around left ventricle in 3D space. Also, the right ventricle is thin-walled com- 
pared to the left ventricle, making the measurements difficult and more sensitive to 
changes from outside the right ventricle (Fig. 1). 
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Fig. 1: Structure and shape of the adult right and left ventricles. 


Although the left ventricle has a better shape for measurements, in real life the right 
ventricle is a basis for artificial pacing and a common place for pacing electrodes, 
thus making cardiac output estimation in the right ventricle potentially more conve- 
nient. Modern heart assisting devices like pacemakers are “smart devices” that can 
regulate the heart according to its condition. To help in the management of the unsta- 
ble patient, measurements of cardiac output need to be taken repetitively [4]. Pacing 
electrodes are inside the right ventricle permanently, hence making constant measure- 
ments possible. Our main goal was to estimate the absolute or relative right ventricle 
volume changes experimentally with impedance spectroscopy and a multielectrode 
intra-ventricular catheter using a setup with an isolated pig heart. Multiple hypothe- 
ses have been developed and evaluated. 

1) Bioimpedance can be used for the measurement of the amount of blood in the 
vicinity of the electrodes, however due to unknown position and shape (bending) of 
the electrode the measured result may not correlate with the actual volume change in 
the whole ventricle. 

2) By using bioimpedance spectroscopy, the distance between the electrode and 
the cardiac tissue could be evaluated based on the changes in £ dispersion area, since 
impedance of blood has negligible reactive components in the used frequency range 
from 1 kHz up to 349 kHz and heart muscle has a much larger reactive component due 
to the dense packing of muscle cells [5]. 

3) With compensation based on the estimated position of the electrode, the 
segment-based volume changes could be measured with better precision. 

The use of a multielectrode intra-ventricular catheter for impedance spectroscopy 
is a potentially useful technique for repetitive measurements to estimate the right ven- 
tricular volume. 
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2 Methods and general principles of design 


We evaluated the right ventricle volume during multiple working cycles with differ- 
ent loading conditions. Usually intra-cardiac impedance is measured with single fre- 
quency but we believe that by using multifrequency measurements, additional infor- 
mation could be gained. The change of volume is measured by inserting an Inquiry 
1110-5-25-M decapolar electrophysiology catheter, shown in Fig. 3, from St. Jude Medi- 
cal into the right ventricle and monitoring the change in impedance between the elec- 
trodes using the four-wire connection method. The electrode positioning is uneven: 
2-5 mm, giving four 2-5-2 (distances between electrodes in mm) measurement combi- 
nations or seven when we also include 5-2-5 combinations. Measured current and volt- 
age values from 2-5-2 and 5-2-5 need separate correlation parameters when translating 
measured impedance data to blood volume (Fig. 2). 
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Fig. 2: Switching between electrode 
modes. 


Fig. 3: Catheter that is inserted into the right ventricle through the pulmonary artery to measure 
impedance change (INQUIRY™ electrophysiology catheter). 


The relationship between volume (V), electrode distance, impedance and blood resis- 
tivity is given as follows: 


1 2 
Yeo oe - (Gx - Gy), (1) 


where p is the resistivity of the blood, L is the distance between measurement elec- 
trodes, a = SV conductance is a stroke volume (SV) correction factor, G, is a total measured 
conductance and G, is the conductance of heart muscle tissue [3]. This equation is 
important for left ventricle volume measurements but research [6] has shown it is ac- 
curate also in the right ventricle. Our experiments show that it is not reliable to use 
the given equation inside the right ventricle due to the unknown electrode position. 
In order to obtain information about the shape and volume of the ventricle during 
operation of the heart, the measurement must be made relatively fast. When using a 
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single frequency, fast measurement is relatively easy to achieve. When using spectrum 
measurement, obtaining real-time measurements becomes more difficult. We have de- 
veloped a real-time impedance spectrum analysis device, Quadra, which uses binary 
multifrequency signals for excitation and measures 15 frequency points, from 1 kHz to 
349 kHz, with a period of 1 ms, giving 1000 measurements per cardiac cycle when the 
heart rate is 60 bpm [7]. 

Excitation in the form of voltage was used with Vpys of 200 mV resulting in root 
mean square (RMS) currents in the range of 0.8 mA to 1.2mA. Since we need to an- 
alyze four or seven sets of electrodes, a fast switching of the measurement device to 
electrode sets is needed. This could be done with 28 relays or one integrated cross- 
point switch. The best results were achieved by using jumping excitation along with 
measurement. Changing excitation is related to longer stabilization time, thus reduc- 
ing the measurement rate. 

We noticed a large measurement noise at frequencies < 10 kHz when multiplex- 
ing both excitation and sense electrodes. The noise decreases exponentially in time 
and can be considered negligible after 100 ms of settling time. This could be due to 
the formation of an electrical double layer since we are measuring in an electrolytic 
environment [8]. 

When the excitation was fixed to the first and last electrodes, and only the position 
of the sense electrodes was multiplexed, the measurement noise from lower frequen- 
cies was removed. The current distribution, in an already difficult environment, gets a 
lot more complex and interpretation of measured impedance would be more difficult. 

For the study of our results we use custom-made analysis software made in Lab- 
VIEW. It shows impedance magnitude values, the pressure value and the phase angle 
in real time. 


3 Measurement setup 


For experiments, we obtained a pig’s heart from the slaughterhouse, which was in 
the condition to be resuscitated. The electrophysiology catheter was placed inside the 
right ventricle through the pulmonary artery together with a silicone tube to manipu- 
late the loading condition and the blood volume. In addition, a tube for the pressure 
sensor was added to measure the blood pressure in the right ventricle. The tricuspid 
valve was then closed to prevent the tubes from moving too much in the ventricle. The 
coronary arteries were kept connected for the external blood supply to keep the heart 
alive. Oxygenated warm blood was used for heart perfusion, since living heart needs 
oxygenated blood of a correct temperature (around 37°C) (Fig. 4). 

For measurements, an intra-ventricular catheter was attached to the real-time 
impedance spectrum device which was connected to the PC and analysis software 
made in LabVIEW (Fig. 5). The right ventricle was then filled with blood and the heart 


A method of intra-ventricular bioimpedance spectroscopy —— 173 
OXYGENATOR / 
ARTIFICIAL LUNG ig 

HEATER / COOLER 

set to 36°C 
02 
BLOOD BATH BLOOD PUMP/ i . 
ARTIFICIAL HEART Fig. 4: Isolated retrogade perfusion heart 
setup. 


IMPEDANCE MULTI- 
MEASUREMENT DEVICE PLEXER 


ECG AMP. 


Fig. 5: Measurement setup with pressure sensor 
and impedance measurement device. ECG is useful 
to monitor heart electrical activity. 


was given an electrical impulse or was contracted by manual squeezing to revive it. 
Changing the blood volume inside the right ventricle gave us impedance results we 
could expect. We used four-wire measurement. Method A consists of switching the 
excitation position and switching the sense position, which gives us a more uniform 
and localized current distribution but also more measurement noise. Method B con- 
sists of a fixed excitation position and switching the sense position, giving a shorter 
acquisition time but a wider and less uniform current distribution. 


4 Results 


A versatile, compact and portable laboratory was constructed for isolated experiments 
together with measurement devices. The isolated heart filling was modulated accord- 


174 — A.Kinketal. 


ing to the level in the tube inserted into the pulmonary artery. The pulmonary valve 
was removed to adapt different loading conditions in the right ventricle. Under low 
loading conditions, prominent changes in bioimpedance were registered and evalu- 
ated in tissue distance changes. Our research shows that phase angle differences were 
the most noticeable. Under mean or high loading conditions, variations in the phase 
angle do not correlate well with value changes inside the right ventricle. 

The initial position and changes in position of the electrodes can be detected with 
bioimpedance measurements. It can be seen from Fig. 6 that the modulus of measured 
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Fig. 6: Magnitude (M1, M3, M5, M7 in Ohms) and pressure (in mmHg) curves measured from the right 
ventricle. 
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Fig. 7: Electrode positions 1and 2 with sensitivity spaces (spheres). 
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Fig. 8: Electrode positions 3 and 4 with sensitivity spaces (spheres). 


impedance shows good correlation with the cardiac pressure. The changes in magni- 
tude of different segments represent the volume of blood and the proximity of the my- 
ocardium. We are able to assess segment-based volume changes in the right ventricle. 
Fig. 7 and Fig. 8 show segments of the heart and how electrode positions can change 
during the heart cycle. Four different positions are proposed with spheres marked as 
electrode sensitivity spaces. 


5 Discussion 


In the present study, we have examined blood volume changes in the right ventricle us- 
ing impedance spectroscopy and a decapolar catheter. We measured electrode shape 
changes, segment-based volume changes and the distance between the measurement 
electrode and the cardiac muscle tissue in the right ventricle by using bioimpedance 
measurements with an intra-ventricular electrophysiological catheter. The proposed 
technique is affected by the complex geometry of the right ventricle and the conduc- 
tivity of surrounding tissues. Using a Baan catheter [3] there is an assumption that the 
electric field produced by the catheter is homogeneous, parallel to the long axis of the 
ventricle, and does not extend to other cardiac chambers. If these presumptions are 
not considered, the errors of ventricle volume estimation are significant. The outcome 
shows a noticeable correlation between the changes of bioimpedance and the varia- 
tions of right ventricle volume, but the interpretation of results can be quite compli- 
cated due to the different electrode positions inside the right ventricle. For acceptable 
results, it is necessary to estimate the electrode configuration inside the right ventricle 
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and separate impedance change contributions from blood and cardiac tissue. This is 
a subject of ongoing research. 
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Mahdi Guermazi, Olfa Kanoun, and Nabil Derbel 
Meat quality monitoring by classification 
based on impedance spectra 


Abstract: In this contribution, we investigate the possibilities of meat state assessment 
based on bioimpedance spectroscopy and fuzzy methods. We report on the procedure 
of feature selection from impedance spectra of meat aiming to reduce measurement 
data to features providing sufficiently information to build the basis for meat state 
estimation, classification of meat type and estimation of the level of freshness. The 
developed method is suitable for signal processing of impedance spectra in a devices 
for real-time meat monitoring based on bioimpedance spectra. 


Keywords: impedance spectroscopy, bioimpedance, fuzzy logic, feature extraction, 
data mining 


1 Introduction 


Meat quality cannot be measured directly, like a classical measurement quantity. It 
includes different aspects related to microbial attacks, the existence of certain sub- 
stances and hormones and structural changes of the muscle during maturation pro- 
cesses. Especially for beef, maturation and structural muscle changes are important 
to realize a tender state by consumption. In order to follow muscle strutural changes, 
electrical measurement methods, such as bioimpedance. Bioimpedance spectroscopy 
is fast, non destructive and can be realized at low-cost. 

Several investigations have been carried out in the field of meat and food diagno- 
sis. Different features of the impedance spectra can be adopted for data mining such 
as the use of model parameters, the use of characteristic points, the extraction of data 
ata fixed frequency or the use ofthe whole measurement results. An overview is given 
about different results realized with data mining and classification methods related to 
the topic of biological tissue monitoring. 

Guerrero etal. [1] have used Cole-Cole model parameters to describe the relation- 
ship between sensory properties and electrical parameters of different meat qualities 
of hams for the semimembranosus (SM) and biceps femoris (BF) muscles. The used pa- 
rameters are Rọ, which corresponds to the electrical impedance at lower frequencies, 
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Ring, Which corresponds to the electrical impedance at high frequencies, a, which cor- 
responds to the shape adjustment parameter, and f-, which corresponds to the charac- 
teristic frequency of the region under measurement at which the imaginary part of the 
electrical impedance has the largest absolute value. A fifth parameter is introduced: 
the ratio between Rint and Rg, which is proportional to the ratio of extracellular water 
to total water in meat. Principal component analysis is used in order to group different 
hams and to locate them in the biplot. They are classified in three groups according to 
the main components: hams having pH45 (defined as the pH at 45 min postmortem) 
lower than 5.85 were classified as pale, soft and exudative (PSE), hams having pHu 
(defined as the pH at 24h postmortem) higher than 6.0 as dark, firm and dry (DFD) 
and the rest as normal hams. The electrochemical impedance spectroscopy (EIS) pro- 
totype correctly detected 69.2 % and 56.0 % for SM and BF muscles, respectively, of the 
problem hams in terms of pastiness. We suggest that the electrical parameters evalu- 
ated in green hams by the EIS prototype could be useful for predicting pastiness in 
dry-cured ham. 

Also, for the application of bovine tissue classification, Negri et al. [2] have used 
features for classification composed by the Cole-Cole parameters Ro, Roo, T and a. Re- 
sults show that using the fitted Cole-Cole model parameters instead of the full spec- 
trum as the input of a neural network can enhance its classification rate and signifi- 
cantly reduce its topology when tested with experimental bovine tissue spectra. 

The same method has been used by Filho etal. [3] for classification to assess 
bovine milk quality. They have confirmed the use of the Cole-Cole parameters may 
have a better classification than the use of the raw acquired spectrum. Using a mul- 
tilayer perceptron (MLP), the recognition percentage can reach 94.6% to classify 
adulterated milk samples with the type of impurity added being water (H,O) or hydro- 
gen peroxide (H,0,). The used MLP is composed by a hidden layer with two neurons, 
an output layer with three neurons defining each class (milk, milk with H,O, milk 
with H,0,). 

Meat monitoring is done using the three main parts: the experimental investiga- 
tions, the modeling and features analysis and finally the meat type and freshness di- 
agnosis (Fig. 1). 


Meat type and 
freshness degree 
determination 


Fig. 1: Monitoring process of meat 
quality. 
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2 Methodology 


The most important and necessary information from a meat monitoring device is the 
type of the muscle and the freshness of the meat. To ralize a suitable signal processing, 
it is necessary to remove spikes in the measurements due to errors. Another important 
factor that should be respected is the use of a reduced input in order to reduce the effort 
and increase the stability of implementation. The process of classification is therefore 
composed of two principal steps: 

— Accurate feature selection from the measured impedance spectra; 

- Classification with high reliability. 


2.1 Feature selection 


Before classfication, generally feature selection should be carried out to reduce the 
measured data to specific data, which serve as the basis for classification. Selection 
of features represents therefore a decisive step in order to reach a high classification 
probability during learning and test phases. Measurements have been carried out with 
different tissues and for along period of time in order to monitor the degradation state 
with time. The measurement setup consists of a laboratory impedance analyzer (Agi- 
lent 4294 A) connected to a personal computer for data acquisition. The spectrometer 
generated a signal within a frequency band ranging from 40 Hz to 110 MHz (Fig. 2). 
Cylindrical penetrating multielectrodes are used for measurements ensuring the re- 
duction of the anisotropy effect and good contacts between electrodes and the tis- 
sue [4]. 

There are different possibilities for selecting features from impedance spectra, 
such as the use of model parameters, the use of characteristic points, the extraction 


Fig. 2: Meat measurement. 
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of data at a fixed frequency or the use of the whole measurement results. According to 
literature and to experimental investigations, measurement results have been evalu- 
ated by means of the modified Fricke-Cole-Cole model. Results are reproducible and 
correspond to the expected behavior due to aging [5]. The model parameter evolution 
shows principally the same behavior, but differs from one muscle to another in the 
absolute value. This is due to several factors, especially the biologic characteristics of 
each muscle. 

Both model parameters (MP) and characteristic points (CP) can be principally 
used for data mining and represent an accurate feature gathered from the measured 
impedance. The investigated possibilities are MP, CP and MP + CP (Fig. 3). 


Model parameters Model parameters Characteristic points 
with CPE 


Re(Z) {CP1} 
Im(Z) {CP1} 


Re(Z) {CP2} 
Im(Z) {CP2} 


Re(Z) {CP3} 
Im(Z) {CP3} 


Fig. 3: Selected examples for feature extraction from impedance spectra. 


According to the literature, principal component analysis (PCA) has been often used to 
reduce features to only the necessary inputs [6-8]. PCA belongs to statistical methods 
of multivariate analysis [9]. It is a method of transforming the initial dataset of high 
dimensions into a new dataset represented by a vector of samples. The goal of this 
transformation is to concentrate the information and reduce its dimension by consid- 
ering the data variance. The PCA method is a linear transformation, as described in 
equation (1), which transfers the original possibly corresponding data set X to a new 
uncorrelated dataset Few. The new dataset F pew has the same dimension as the orig- 
inal dataset X; A, is the coefficient transformation matrix. We have 


Faew = Ark: (1) 


The first few components (columns) of the new dataset F pew contribute the majority 
of variances and are considered to be the principal components, which can reflect the 
original data matrix without big information losses. They have smaller dimension and 
are uncorrelated to each other, which makes the analysis of the data much easier. The 
input data set X, is in the form of a matrix, as shown in equation (2), where p is the 
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number of the parameters and n is the number ofthe observations. We have 


Xu An Amp 
Xg=| Xn X2 Xy) (2) 
Xm Xn2 Xnp 


The first two components for the feature composed by MP, for the CP and for the 
MP + CP explain more than 90 % of the variances. In order to verify which feature we 
should use for classification, a simple idea is used consisting of the principle of su- 
perposition of the test data in the corresponding region of input data after k-nearest 
neighbors (KNN) calculation. The algorithmic method of KNN was used to give a com- 
plete 2D presentation in a geometrical plot of the data. 


2.2 Classification method 


Different suitable methods of classification can be used in this application. Fuzzy logic 
is selected as a classification method that provides a successful alternative for easy 
classification, as it is in general a robust method. It has also the advantages to deal 
with uncertainty problems and usually gives good results for various areas. Fuzzy logic 
is widely used in medical applications for diagnosis of tuberculosis, cancer, asthma, 
diabetes, aphasia, malaria, HIV, pulmonary embolism, cortical malformations and 
pancreatic diseases [10]. Fuzzy logic has been implemented for an earlier detection 
of breast cancer [11, 12]. Results are found to be very useful and helpful for oncolo- 
gists, radiologists and doctors. Thus, fuzzy logic, as an intelligent method, will be of 
great help especially when saving patient lives is possible. The process of fuzzy logic 
based classification includes three principal steps [13]: 

-  Fuzzification of input variables and choice of suitable membership functions; 

- Construction of fuzzy rules and making decisions; 

—  Defuzzification. 


An explanatory illustration of the fuzzy inference process is shown in Fig. 4. 

In fact, the model parameters (received data) undergo a feature’s extraction pro- 
cess (information), which is a preparatory step for the classification method. The re- 
sulting data are then fuzzified using membership functions. Afterward, fuzzy values 
are analyzed using fuzzy rules to get the fuzzy outputs. After defuzzification, the new 
data are then compared and classified with similarities to get the desired clusters. 

Fuzzification is the conversion of crisp input data into fuzzy values or linguistic 
variables using membership functions. These fuzzy values reflect the human perspec- 
tive of the given system. These fuzzifiers are membership functions determined de- 
pending on a subjective view of the problem and an individual’s perception of the sit- 
uation. Depending on the desired class and performances, simple functions are used 
to build membership functions, such as triangular, trapezoidal and Gaussian ones. 
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Fig. 4: Block structure of a Fuzzy system. 


Fuzzy rules, which have the general form of “IF X is A THEN Y is B”, are destined to 
map the fuzzy input X to the fuzzy output Y. In a logical context, X is referred to be a lin- 
guistic variable and A is a linguistic value. The ‘IF’ part of the rule is called antecedent; 
the ‘THEN’ part is called consequent. Connection between the rules is achieved using 
logical operators AND, OR and NOT. AND or min takes always the minimum of the 
functions value. OR (max) takes the maximum of the presented functions. 

Defuzzification is the conversion of a fuzzy set quantity to a precise and deter- 
ministic set quantity. This step is necessary because in several applications, where 
the implementation is required, machines could only work with crisp or binary values 
and not with linguistic variables. 

After defuzzification, the new data are compared and classified to get desired clus- 
ters [14]. Two important types of fuzzy logic systems exist and are much used for diag- 
nosis: Mamdani inference [15, 16] and Takagi-Sugeno (TS) inference [17]. 

The Mamdani method is suitable for nonlinear systems, characterized by a fuzzy 
set consequent in the fuzzy rule that has the following form: 

R : IF L(y is A),....,8, is Ay) THEN (y; is B4,....,Yn is By), 

where: 

—  X,.., X, are the inputs of membership functions; 

- Aj,.., A, are the fuzzy sets of the antecedent part; 

— Yis- Yn are the outputs of the membership functions; 

— B,,...,B, are the fuzzy sets of the consequent part and L is the logical function 
connecting the propositions. 


For the TS method, the output membership function or more precisely the rule conse- 
quent is only linear or constant. If the rule consequent is constant we call the model 
a zero-order Sugeno model, which has the following form: 

Ro: IF L(x is A4, .. ., Sk is Ay.) THEN y = k, 

where k is a crispy constant defined by the operator. 


Meat quality monitoring by classification based on impedance spectra —— 183 


If the rule consequent is a linear function the Sugeno model is called a first-order 
Sugeno model with the following expression: 
R, : IF LQ is Aj,...,8; is Aj) THEN y = f(X,..Xn)- 


3 Results and discussion 


3.1 PCA computing 


A database composed of six muscles is used in order to use the PCA method well. This 
database will serve to determine the suitable feature. Three different muscles of two 
animals of different ages are used. The test data should be unknown. For that, two 
spectra will be used for learning and two others for testing in every case. 

Singular-value decomposition is chosen following the study done in [18]. The de- 
composition of the normalized matrix “X” is done according to 


X" = USV’, (3) 


where S is the n x p diagonal matrix with singular values g; of the matrix X in diagonal 
form. We also have 


Si =0; (4) 
and singular values are already sorted in a decreasing order, i. e., 


01 > 0, > Oy >> Op. (5) 


The relationship between the singular value g; and the eigenvalue A is given by 


o; = Vii (6) 


where U is the n xn left singular-valued vector matrix of X, V is the p xp right singular- 
valued matrix of X and both are orthogonal matrices. The principal components F can 
be calculated by equation (7), i. e., 


F=XV=USV'V=US. (7) 


Before dealing with PCA computing, the inputs and test data of the three cases, MP, CP 
and MP + CP are defined. The test data have exactly the same dimension as the input 
data. 

The result of the superposition of model parameter test data on the KNN partition 
corresponds to an accuracy equal to 100%, as shown in Fig. 5. The accuracy for the 
characteristic points is 66 % with four erroneous outputs. For the feature composed 
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Fig. 5: Superposition of model parameter test inputs to the KNN graph. 


by the characteristics points and model parameters we get 12 correct outputs that cor- 
respond to an accuracy of 100 %. 

Feature extraction based on the use of model parameters and characteristic points 
shows good results but with the consideration of the use of reduced inputs in the net- 
work we choose to work with the model parameters only for the classification that also 
gives an accuracy of 100 %. 


3.2 Zero-order Sugeno model 


The classification is carried out using the TS method, which is similar to the Mamdani 
model in fuzzifying the inputs and applying the fuzzy operator. The major reason of 
this choice is that the TS method output membership can be either constant (zero- 
order) or linear (first-order), corresponding exactly to the application specifications. 
The zero-order TS method with constant output memberships is exactly the solution 
that permits to get an accurate diagnosis and that corresponds to the problematic. For 
each muscle we appoint a corresponding constant number: 

— 1 for the muscle LD beef (M1); 

— 2 for the muscle RA beef (M2); 

— 3 for the muscle SM veal (M3); 

— 4 for the muscle LD veal (M4). 


And for each muscle we appoint a corresponding constant number as output: 
— 1 for fresh meat (F); 

— 2for edible meat (E); 

— 3 for critical meat (C); 

— 4 for dangerous meat (D). 
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Fig. 6: Fuzzy logic process. 


The fuzzy logic process (Fig. 6) is composed of two fuzzy logic systems that have been 

used for meat diagnosis; the first for the classification of the muscle type, the second 

for finding the freshness, to classify the meat in four classes that correspond to four 

different periods of time: fresh meat, edible meat, critical meat and dangerous meat. 
The output membership function or more precisely the rule consequent is only 

constant. We have 

Ro : IF L(x is A,,...,5K 15 Ay) THEN y = k. 


3.2.1 Membership functions 


The first step of the muscle type classification consists of the choice of the input and 
output variables and selecting carefully their intervals’ ranges. By this step, we try to 
collect the data which generate the same output in one interval range. This process is 
carried out for all the inputs. In one side, we define the three model parameters R,, 
R; and C, extracted from bioimpedance measurements as the inputs [5]; the output 
corresponds to the muscles LDpeer, RApeefo SMyeaı and LDyeq). After several trials for the 
muscle type classification, we select the Gaussian membership functions (MF) defined 
by a central value c and a standard deviation ø > 0. We have 


=(x-0)? 


Ha) =f@%0,0 =e. (8) 


For each of the three inputs R,, R; and C, we attribute the Gaussian membership func- 
tions associated with the intervals’ ranges extracted from the database. It is necessary 
to use different zones for each input due to the performance of the learning. We at- 
tribute five membership degrees to the inputs Re and C and three to the input R;. For 
each input the corresponding number of MF that leads to a high accuracy and that 
avoids the overlapping for the hole database represents the type of the meat. 

For the fuzzy logic of freshness classification, the major difference with the fuzzy 
logic of muscle classification is that we do the learning for each muscle separately 
in order to determine the aging that corresponds to its freshness. The three model 
parameters R,, R; and C are defined as inputs and affected the trapezoidal MF after 
several trials. 
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The fuzzy logic classification results show an accuracy of 100 %. The result of each 
muscle is in the close interval of the corresponding value; for example, for the veal 
SM, the accepted result that corresponds to this muscle should be in the range from 
2.5 to 3.5. 

The fuzzy logic classification for noisy inputs reaches the same results as the used 
inputs. In the diagnosis process, the choice of confidence areas is decisive in order to 
support the learning process. We select and define the period of time such as: 

- fresh, between day 2 and day 3; 

— edible, between day 4 and day 6; 
- critical, between day 7 and day 10; 
- dangerous, more than day 11. 


According to the defined periods, the accuracy is 84.62 % for the beef LD, 92.31% for 
the beef RA, 100 % for the veal SM and 61.54 % for the veal LD. 


4 Conclusion 


Feature extraction using model parameters leads to a high accuracy. Suitable data 
for meat diagnosis are able to correctly classify muscle types and states. Fuzzy logic 
allows for an easy method of classification. It correctly classifies the inputs to the cor- 
responding classes for the first step of meat type with a recognition index equal to 
100%. 

For the freshness classification, fuzzy logic reaches a recognition index equal to 
84.62 % for the beef LD, 92.31 % for the beef RA, 100 % for the veal SM and 61.54 % for 
the veal LD. 

Fuzzy logic has many advantages which can be beneficial in the monitoring pro- 
cess of meat quality. These advantages are especially the consideration of noisy in- 
puts, the tolerance that leads to a compact program according to the conception pro- 
cedure and finally the simplicity. 
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